Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 21.2% intergenic (‑194/+114) 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,9840TG21.2% 651.2 / 180.8 251intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (103/94);  new base G (51/2);  total (155/96)
Fisher's exact test for biased strand distribution p-value = 2.09e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 253 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083
                                                                                                 |                                                                                                         
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  2:1473459/1‑100 (MQ=255)
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                         >  1:943263/1‑100 (MQ=255)
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCGTGGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                        >  2:157112/1‑100 (MQ=255)
   gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                      >  2:1244809/1‑100 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                  <  2:3153067/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                 <  2:2714226/100‑3 (MQ=37)
           gcgcTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                             <  2:2203751/100‑1 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                               >  1:2172834/1‑99 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACGAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                                            >  1:1928874/1‑100 (MQ=255)
             gttGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                           >  2:2338309/3‑100 (MQ=255)
              cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                                          >  2:2641531/1‑100 (MQ=255)
               tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                         >  1:338329/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                      <  1:1138382/100‑1 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGt                                                                                      >  2:1786807/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                     >  2:3349735/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                     <  2:1902549/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                     >  1:457585/1‑100 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                     <  1:38623/100‑1 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                   >  2:2747907/1‑100 (MQ=255)
                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                                  >  1:2655090/1‑99 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                              <  2:383993/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                             <  2:2029878/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                             <  2:1145506/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                           <  1:903231/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                       <  2:2465653/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                       >  1:1250060/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                      <  1:1214455/100‑1 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                      >  2:1359710/1‑100 (MQ=255)
                                   cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                     >  1:2231321/1‑100 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                    <  1:1988304/100‑1 (MQ=255)
                                         gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At                                                             <  2:2238163/100‑1 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt‑‑t                                                       >  2:2394041/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCg                                                     >  2:3120581/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCg                                                     <  1:2920629/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCg                                                     >  2:1981182/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTAGTTTATGCCGGATACGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt‑‑t                                                       >  1:79632/1‑98 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGt                                                    >  1:2443640/1‑100 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGt                                                    >  2:1051815/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGTCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTa                                                   >  1:203407/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGt                                                    <  1:1608798/99‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTa                                                   <  2:1603235/100‑1 (MQ=255)
                                                   cataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGc                                                >  1:10325/6‑100 (MQ=18)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCt                                              >  1:732815/1‑100 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTg                                             <  2:563523/100‑1 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGa                                            >  1:3166985/1‑100 (MQ=255)
                                                       gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TAGTAGGCCTGa                                            >  1:2273254/1‑100 (MQ=255)
                                                          aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATaa                                         <  1:1316700/100‑1 (MQ=255)
                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAg                                        >  1:2542253/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGa                                    >  1:735502/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGa                                    <  1:2338290/100‑1 (MQ=255)
                                                                cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAg                                   <  2:3317055/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAgc                                  >  1:410848/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAgcg                                 <  1:63398/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAgcg                                 <  1:2832996/100‑1 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAgcg                                 >  2:876262/1‑100 (MQ=255)
                                                                  cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑CCGTAGGCCTGATAAGCGGAg                                   >  1:963954/1‑98 (MQ=255)
                                                                      ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑AT‑‑TGGTAGGCCTGATAAGCGCAGCGCATc                             >  2:2006981/7‑100 (MQ=21)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCa                            <  2:2599090/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                           <  1:1474492/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                           >  2:2879289/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑ACGTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:750789/5‑100 (MQ=17)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑ACGTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:249762/5‑100 (MQ=17)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑‑T‑GTAGGCCTGATAAGCGTAGCGCATCAg                           >  2:2085240/5‑100 (MQ=18)
                                                                         taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TC‑‑ACGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:1300739/4‑100 (MQ=17)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:3165432/100‑1 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2526943/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:1027473/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑CA‑‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:3262955/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  2:1264129/3‑100 (MQ=18)
                                                                           tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  1:2687757/1‑100 (MQ=255)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑‑A‑TT‑CAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  2:825341/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CT‑‑GCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:771773/2‑100 (MQ=17)
                                                                            gacggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑AT‑‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  2:344159/10‑99 (MQ=25)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:2185395/1‑100 (MQ=255)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      <  2:2777647/100‑1 (MQ=255)
                                                                              tggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttg                    >  1:977133/8‑99 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  2:649910/7‑98 (MQ=21)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTGAATATATTGCA‑GA‑AA‑CC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:2986309/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑CC‑‑ACGTAGGCCTTATAAGCGTAGCGCAAAAGGCAgttt                    >  1:2962707/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TT‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  2:505524/7‑98 (MQ=25)
                                                                                gataaggcgTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                  >  2:586030/6‑97 (MQ=21)
                                                                                gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATGG‑GC‑CA‑‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAatttt                   >  2:2554182/6‑100 (MQ=18)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑AC‑‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                   >  2:1687849/6‑100 (MQ=18)
                                                                                  taaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CC‑‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                 >  1:828953/4‑100 (MQ=18)
                                                                                  taaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CA‑‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgc                 >  2:1283102/4‑100 (MQ=18)
                                                                                  taaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CT‑‑GCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttggc                 >  2:2897580/4‑97 (MQ=18)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:2255934/100‑1 (MQ=255)
                                                                                  tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  1:2964701/100‑1 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               <  1:1786791/100‑1 (MQ=255)
                                                                                     ggcgTGAACGCCATAGCCGGCCGAC‑AAAATCGTGCAGATTCAGTACAGTGCA‑TA‑AG‑GG‑‑GCGTATGCCTGAGAAGCGGAGCGCATCAGGCATTtttgcctt              >  1:2017791/1‑100 (MQ=12)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  1:3447954/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  2:977525/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         <  2:338821/100‑1 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  2:848318/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGAATCAGGCATTTTTGCCTTTGCCa         >  1:1192970/1‑98 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:44088/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:1359698/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  1:2310902/100‑1 (MQ=255)
                                                                                               ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:568494/1‑100 (MQ=255)
                                                                                                 ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  2:1503866/100‑1 (MQ=255)
                                                                                                 |                                                                                                         
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TT‑‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705887‑3706083

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: ‑