Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,548,571 G→T 15.2% T500T (ACC→ACA)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,5710GT15.2% 140.0 / 10.0 57T500T (ACC→ACA) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base G (33/15);  new base T (2/7);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACA  >  NC_000913/3548423‑3548719
                                                                                                                                                    |                                                                                                                                                    
gCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACcggtacg                                                                                                                                                      >  8:201022/1‑149 (MQ=255)
    gTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCggtacggtac                                                                                                                                                  <  6:384253/149‑1 (MQ=255)
         gtgtCGCCTCAACCTCCATCGACTGCTCCGGAAACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTAGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCTCTGCTGCGCAGCTTACCGACAATCTCTACCggtacggtaca                                                                                                                                                 >  4:301387/1‑144 (MQ=255)
         gtgtCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCATCTTACCGACAATCTCTACCGGTACGGTAACgag                                                                                                                                              >  1:387242/1‑148 (MQ=255)
         gtgtCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGAACTTTGTAAGAGTACACACCGCTGCTGCGCATCTTACCTACAATATCTACCGGTACTGTAcc                                                                                                                                                 >  6:103664/1‑145 (MQ=255)
                        ccATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCAcc                                                                                                                              <  3:217918/149‑1 (MQ=255)
                          aTCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCa                                                                                                                             >  4:247933/1‑148 (MQ=255)
                           tCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATa                                                                                                                           >  8:106517/1‑149 (MQ=255)
                             gACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTAACCACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACAATAc                                                                                                                          >  8:126272/1‑148 (MQ=255)
                                       gATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCAc                                                                                                                               >  6:49257/1‑133 (MQ=255)
                                        aTACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                                                                                              >  1:306144/1‑149 (MQ=255)
                                         tACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGttt                                                                                                             <  7:201022/149‑1 (MQ=255)
                                            gCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACt                                                                                                          >  2:325840/1‑149 (MQ=255)
                                                       aCGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGt                                                                                                   <  3:76262/145‑1 (MQ=255)
                                                       aCGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGt                                                                                                   >  4:76262/1‑145 (MQ=255)
                                                          gAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCACTCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGa                                                                                             >  7:374979/1‑148 (MQ=255)
                                                              gTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTAAACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATATTCAACAATCACCATACAGACATGACGTTTACTTTCTAGTGCCACAATCg                                                                                         >  2:157894/1‑148 (MQ=255)
                                                                cttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCga                                                                                        >  7:147402/1‑147 (MQ=255)
                                                                      gTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTACCCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGagag                                                                                     >  3:39903/1‑144 (MQ=255)
                                                                        ggTCGTTTTCGAAACCGAGCACTTTGTAAGAGTACACACCGCTTCTGCGCAGCTTTTTTTTTTTATTTATAGGGAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatca                                                                              <  4:245421/149‑1 (MQ=255)
                                                                         gTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatcat                                                                             >  3:315993/1‑149 (MQ=255)
                                                                              tttCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTg                                                                                                  <  5:214634/123‑1 (MQ=255)
                                                                              tttCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTg                                                                                                  >  6:214634/1‑123 (MQ=255)
                                                                                  gAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGa                                                                     >  8:88941/1‑148 (MQ=255)
                                                                                   aaaTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCCCCAATCCCCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGa                                                                     >  7:77860/1‑147 (MQ=255)
                                                                                    aaTAGAGCACTTTGTAAGATTACACACCGCTGCTGCGCAGCTTTTTTTTTTTATTTATAGGGAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAg                                                                  <  2:306144/149‑1 (MQ=255)
                                                                                    aaTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTCACGACAATCTATACCGGTACTGTACCGAGATCTTCACCAATCACCATACATCGAAGACGTTTACTTTCGAGTGCCAGAATCGAGCGCAGATCATCCACCGGAAAg                                                                  >  1:285877/1‑149 (MQ=255)
                                                                                     atatAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCTACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGt                                                                 <  1:157894/149‑1 (MQ=255)
                                                                                     aTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCTCAGCTTTTTTTTTTTATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGt                                                                 <  2:188208/149‑1 (MQ=255)
                                                                                     aTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc                                                                         <  1:309079/141‑1 (MQ=255)
                                                                                     aTAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCAAACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc                                                                         >  2:309079/1‑141 (MQ=255)
                                                                                      tAGAGCACTTTGTAAGAGTACACACCGCTGCTTCGCAGTTTTTTTTTTTTATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTg                                                                <  2:387242/149‑1 (MQ=255)
                                                                                      tAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGt                                                                 >  3:44237/1‑148 (MQ=255)
                                                                                        gagCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATAACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACAGGATAGTGaa                                                              >  1:166923/1‑149 (MQ=255)
                                                                                        gagCACTTTGTAAGAGTACACACCGCTGCTGCGCAGATTATTTTTTTTATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGaa                                                              <  8:147402/149‑1 (MQ=255)
                                                                                           cACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCCCCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATACTGACCa                                                            >  5:67451/1‑148 (MQ=255)
                                                                                                        tACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGAACTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGACAGCAGATCATCCACCGGATAGTCAACAAACCCGACCTGAt                                               >  6:275282/1‑148 (MQ=255)
                                                                                                        tACACACCGATGCTGCGCAGCTTACCGACAATCTCTACCGGTACGCTACCGAGATCTTCACCAATCACCATACAGCCATGACGATTACTTTCGAGTGCCAGAATCCAGAGCAGATCAACCACCGGATAGTGAACATACGCTCCCTGAt                                               >  2:358735/1‑148 (MQ=255)
                                                                                                          cacacCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTg                                            >  5:92419/1‑149 (MQ=255)
                                                                                                          cacacCGCTGCTGCGCAGATTATTTTTATTATTTATAGGTAGTGTACCGAGATCTACACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCt                                             <  6:67451/148‑1 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTTTCTTTAATATTTATAGGTAGTGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAc                                                         <  2:377177/133‑1 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCgaga                                                                                      >  8:108620/1‑104 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCgaga                                                                                      <  7:108620/104‑1 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAAc                                                             >  4:73355/1‑129 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAAc                                                             <  3:73355/129‑1 (MQ=255)
                                                                                                             acCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACCTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAc                                                         >  1:377177/1‑133 (MQ=255)
                                                                                                               cGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGtc                                       >  8:286985/1‑149 (MQ=255)
                                                                                                               cGCTGCTGCGCAGCTTACAGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCATAATCGAGAGCAGATCATCCACCGTATAGTGAACATACGCGCCCTGATCTGCCGtc                                       <  5:103664/149‑1 (MQ=255)
                                                                                                                    ctgcGCATCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGa                                                <  1:83037/135‑1 (MQ=255)
                                                                                                                       cgcAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAcgc                                                       <  1:120160/125‑1 (MQ=255)
                                                                                                                       cgcAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATAcgc                                                       >  2:120160/1‑125 (MQ=255)
                                                                                                                        gcAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGAACATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATAc                              >  2:223393/1‑149 (MQ=255)
                                                                                                                             ttACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTAt                         >  4:210762/1‑149 (MQ=255)
                                                                                                                             ttACCGACAATCTCTACCGGTACGAAACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTAt                         >  4:317085/1‑149 (MQ=255)
                                                                                                                               aCCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATc                        >  4:43896/1‑148 (MQ=255)
                                                                                                                                 cGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATccac                     <  3:210762/149‑1 (MQ=255)
                                                                                                                                       tctcTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACaaa               >  6:95414/1‑149 (MQ=255)
                                                                                                                                             ccGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGACCATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACgcagc         <  1:358735/149‑1 (MQ=255)
                                                                                                                                             ccGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACgca           <  5:95414/147‑1 (MQ=255)
                                                                                                                                             ccGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACgca           <  3:247933/147‑1 (MQ=255)
                                                                                                                                               ggtacggtacCGAGAGCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCCTCCACCTGATAGTGAACATACGCGCCCTGATCTGCCGTCTCTCCATACGGTATCCACCACAAACGCAGCAt       <  3:301387/149‑1 (MQ=255)
                                                                                                                                                gtacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATGTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATc      <  5:275282/149‑1 (MQ=255)
                                                                                                                                                gtacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATc      >  2:209120/1‑149 (MQ=255)
                                                                                                                                                gtacggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATc      <  1:223393/149‑1 (MQ=255)
                                                                                                                                                   cggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGa    <  7:106517/148‑1 (MQ=255)
                                                                                                                                                    ggtacCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAAACGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACa  <  8:374979/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGATAGTGAACATACGCGCCCTGATCTGCCGTCTCGCCATACGGTATCCACCACAAACGCAGCATCGACA  >  NC_000913/3548423‑3548719

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑