Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 29.9% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC29.9% 295.3 / 149.2 179intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (58/67);  new base C (46/7);  total (105/74)
Fisher's exact test for biased strand distribution p-value = 3.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-01

Only 100 of 179 total aligned reads displayed.

GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCA  >  NC_000913/3705890‑3706079
                                                                                               |                                                                                                  
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcagat                                                                                             >  1:321559/1‑96 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             >  1:660360/1‑100 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  2:1714430/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:2414257/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:2031458/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                            <  1:1642073/100‑2 (MQ=37)
     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                           <  1:1561210/100‑3 (MQ=37)
        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTata                                                                                        <  1:1866528/100‑3 (MQ=255)
        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                       >  1:1003293/1‑100 (MQ=255)
          gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                     >  2:2176277/1‑100 (MQ=255)
            tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                   <  2:6324/100‑1 (MQ=255)
                cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                               >  2:2527939/1‑100 (MQ=255)
                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                              <  2:250606/100‑1 (MQ=255)
                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTG‑TTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                            >  1:1979869/1‑100 (MQ=255)
                     agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                          <  2:2110829/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                        <  2:1756539/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                        <  1:327576/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                        <  1:2665575/100‑1 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                       >  2:666547/1‑100 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                       >  1:1308065/1‑100 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAACTTc                                                                        >  2:200813/1‑99 (MQ=255)
                            gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                   <  1:1941499/100‑1 (MQ=255)
                              gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                 >  1:354841/1‑100 (MQ=255)
                                cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                               <  1:1694706/100‑1 (MQ=255)
                                cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                               >  2:2631388/1‑100 (MQ=255)
                                 aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                              <  1:1133869/100‑1 (MQ=255)
                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑g                                                           <  2:768951/100‑1 (MQ=255)
                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑ttt                                                   >  2:206616/1‑100 (MQ=255)
                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑ttt                                                   >  2:241697/1‑100 (MQ=255)
                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                                 <  1:729138/100‑1 (MQ=255)
                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCg                                                 <  1:1605600/100‑1 (MQ=255)
                                             gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTa                                               >  2:46913/1‑100 (MQ=255)
                                               ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAgg                                             <  2:1345977/100‑1 (MQ=255)
                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTCC‑AAAGTCGTGCAAGTTCAATACATTGTA‑GA‑Ag‑ttt                                                   >  2:2412171/7‑94 (MQ=14)
                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGc                                            <  1:699052/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGcc                                           <  2:710041/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGcc                                           <  2:2429716/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGc                                            >  1:1061909/1‑99 (MQ=255)
                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCt                                          <  1:554829/100‑1 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTg                                         <  1:300445/100‑1 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTg                                         >  1:1220658/1‑100 (MQ=255)
                                                    gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGa                                        >  1:407172/1‑100 (MQ=255)
                                                    gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGa                                        >  2:942258/1‑100 (MQ=255)
                                                     aCAAAGTGTTCTTTGTTGATGCCGGAGGCGGCGTTAACGCCTTATCCGGCCTAC‑ATAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGAt                                       <  2:606249/100‑1 (MQ=18)
                                                     aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGAt                                       <  1:2648097/100‑1 (MQ=255)
                                                         aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGc                                   <  1:2450090/100‑1 (MQ=255)
                                                          gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCg                                  <  1:2591667/100‑1 (MQ=255)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                 <  2:1711901/100‑1 (MQ=255)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCgg                                 <  1:501753/100‑1 (MQ=255)
                                                             cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAg                               >  2:2447668/1‑100 (MQ=255)
                                                               cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATaagcgaagcg                             >  2:308633/1‑100 (MQ=255)
                                                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAgtg                             <  1:1842417/98‑3 (MQ=255)
                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                         >  1:247865/1‑100 (MQ=255)
                                                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAgg                      >  1:534579/4‑100 (MQ=17)
                                                                      tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                      >  1:2365747/1‑100 (MQ=255)
                                                                       aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:1392348/3‑100 (MQ=18)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                     >  2:101518/3‑100 (MQ=21)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGTTAATTCAATATATTGCA‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGc                     >  2:6433/3‑100 (MQ=18)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:2052756/3‑100 (MQ=21)
                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                     >  1:1507151/1‑100 (MQ=255)
                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                     >  1:1507962/1‑100 (MQ=255)
                                                                        ggctggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATTGTGCAAAATCAATGGATTGCA‑CA‑GC‑TAGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  2:840190/11‑100 (MQ=21)
                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                   >  1:483198/2‑100 (MQ=18)
                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                     >  1:287713/2‑99 (MQ=17)
                                                                        ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCACTATATTGCA‑GA‑GA‑TTGCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                    >  2:1641264/2‑100 (MQ=11)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:2158119/2‑100 (MQ=18)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                    >  1:2152905/2‑100 (MQ=21)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCa                    >  1:2281063/2‑100 (MQ=18)
                                                                        agCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                   >  1:225696/2‑100 (MQ=18)
                                                                         gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑GAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                   >  1:1538163/10‑99 (MQ=21)
                                                                         gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                   >  1:594792/1‑100 (MQ=255)
                                                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                   >  1:1011896/1‑99 (MQ=18)
                                                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGTTAAGCGTAGCGCATCAGGCag                   >  2:537208/1‑99 (MQ=17)
                                                                          ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GA‑GA‑TCATGTAGGCCTGATAAGCGTAGCGCATAAGGCaat                  >  2:2315923/9‑98 (MQ=14)
                                                                          ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                  >  2:87525/9‑98 (MQ=18)
                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                  >  2:551645/1‑100 (MQ=255)
                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGCagc                  >  2:2251430/1‑98 (MQ=18)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑CCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                 >  2:1117463/8‑97 (MQ=18)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                >  1:192492/8‑99 (MQ=21)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:1768415/7‑100 (MQ=18)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:85771/7‑100 (MQ=17)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagct                >  2:2446629/7‑96 (MQ=21)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CTACGTAGGCCTGATAAGCGTAGCGCATCCGGCAgttt                >  2:1073053/7‑100 (MQ=17)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                >  2:2731441/7‑96 (MQ=25)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑TTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:2690385/7‑100 (MQ=17)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GA‑GA‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                >  1:1709520/7‑100 (MQ=17)
                                                                              aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg              <  2:144561/100‑1 (MQ=255)
                                                                               taaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc             >  2:2112853/4‑97 (MQ=17)
                                                                               tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt               >  2:1544002/1‑98 (MQ=255)
                                                                               tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc             <  1:1879496/100‑1 (MQ=255)
                                                                                gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc            <  1:1086667/100‑1 (MQ=255)
                                                                                  ggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgtt          <  1:1107760/100‑1 (MQ=17)
                                                                                  ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTACAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt          >  2:1162677/1‑100 (MQ=255)
                                                                                    cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg        >  1:312692/1‑100 (MQ=255)
                                                                                     gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc       <  2:1507955/100‑1 (MQ=255)
                                                                                      ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc      <  1:2527957/100‑1 (MQ=255)
                                                                                      ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc      <  2:2221099/100‑1 (MQ=255)
                                                                                       tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa     >  2:1454594/1‑100 (MQ=255)
                                                                                       tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa     <  2:748485/100‑1 (MQ=255)
                                                                                         aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc   <  2:1199266/100‑1 (MQ=255)
                                                                                               |                                                                                                  
GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCA  >  NC_000913/3705890‑3706079

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑