Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 41.4% 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,9850TC41.4% 234.0 / 304.8 242intergenic (‑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 (70/72);  new base C (89/11);  total (159/83)
Fisher's exact test for biased strand distribution p-value = 3.94e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-02

Only 100 of 246 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                            
gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                             >  1:2116179/1‑100 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                           >  2:519938/1‑100 (MQ=255)
      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                       <  1:2809485/100‑1 (MQ=37)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                      <  1:1699370/100‑2 (MQ=37)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                     <  2:1683930/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                                     <  1:221226/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                     <  2:743460/100‑3 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑a                                                                                                  <  1:213316/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                                  <  1:2517435/100‑1 (MQ=25)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                                 >  1:161693/1‑100 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                                 >  2:3460089/1‑99 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                               <  2:1529982/100‑1 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                                             >  1:2248387/1‑100 (MQ=255)
                gCACGCAGAGGCACAGCGCCATCGGGTAGACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                                             <  2:2831839/100‑11 (MQ=18)
                gCACGCAGAGGAACAGCGCCACCGGGTATACAGATAGTGGACAAAGTGCGCTTTGTTTATGCCGAATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                                              <  2:1717787/99‑10 (MQ=18)
                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                                         <  2:2057814/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                       <  1:1370036/100‑1 (MQ=255)
                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                                       >  2:613041/1‑100 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                      <  1:1974731/100‑1 (MQ=255)
                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                                      <  2:3582598/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                                <  1:3035695/99‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCa                                                                               >  1:3140673/1‑100 (MQ=255)
                            gcaGCGCCATCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                              >  1:3643810/2‑100 (MQ=255)
                            acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                                              >  1:1794866/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                             <  2:671520/100‑1 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                                             <  2:310388/100‑1 (MQ=255)
                                 gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                                         >  1:134892/1‑100 (MQ=255)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAtt                                                                        >  1:950880/1‑100 (MQ=255)
                                          tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t                                                             >  1:2340773/1‑100 (MQ=255)
                                           ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t‑t                                                           >  1:2943577/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                      <  1:1469346/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                      <  2:647631/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑g                                                      <  1:3068839/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTa                                                    >  2:3064913/1‑100 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTa                                                    >  2:3407641/1‑100 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGcc                                                <  1:3311824/100‑1 (MQ=255)
                                                    ataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGcc                                                >  1:2819664/5‑100 (MQ=18)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCt                                               <  2:1056185/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATa                                           >  1:93419/1‑100 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGc                                        <  1:3213095/100‑1 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGAA‑GA‑AT‑T‑TTC‑‑GTGGGACTGATatg                                         >  2:3579623/1‑97 (MQ=17)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                       >  1:3020061/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                       >  1:115850/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCg                                       <  2:1865346/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GG‑AT‑T‑TTC‑‑GTAGGCCTGCTAAGCGGa                                     >  1:3262708/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGa                                     <  1:805363/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAg                                    <  2:1587848/99‑1 (MQ=255)
                                                                  agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑CAC‑‑GTAGGCCTGATAAGCGAAg                                    >  2:754061/11‑100 (MQ=2)
                                                                   tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAgcgc                                 >  2:328768/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                             >  1:1692790/1‑100 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  1:1925295/5‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAATATATTGCA‑GA‑AA‑T‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1822416/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCagatcgga                         >  2:1140318/3‑93 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1847401/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                          >  1:1833242/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc                          >  2:488448/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACTCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑T‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3564957/3‑100 (MQ=11)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1988866/3‑100 (MQ=21)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑T‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:1032579/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:2900952/3‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCC‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:2154372/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2437847/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:145228/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2463804/2‑99 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:13546/2‑100 (MQ=18)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                        >  1:1772547/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCac                        >  1:736990/1‑99 (MQ=17)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2319068/1‑99 (MQ=18)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2222745/1‑100 (MQ=18)
                                                                             ccggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                       >  1:572264/9‑98 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  1:1554762/1‑100 (MQ=255)
                                                                             ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                       >  2:2763196/1‑98 (MQ=21)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTATGCCTGATAAGCGTAGAGCATCAGGCa                         >  1:1318405/8‑97 (MQ=17)
                                                                              cggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCACGCAAACTCAATATATTGCA‑GAGAT‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                      >  1:2044099/8‑97 (MQ=17)
                                                                              cGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                      <  2:2670477/100‑1 (MQ=255)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑AT‑CCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGCCAtttgc                    >  2:694816/7‑98 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1127222/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2485018/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:1005414/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1740579/7‑100 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGAAgttt                     >  1:3546941/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  1:2577577/7‑96 (MQ=25)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCA‑GG‑GA‑C‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                     >  1:459185/7‑96 (MQ=21)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GA‑GA‑T‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  1:3495876/7‑100 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                    >  1:2710035/6‑100 (MQ=17)
                                                                                gataaggcgTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttg                    >  1:123204/6‑99 (MQ=18)
                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                   <  1:1065965/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  2:3140650/100‑1 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:620544/1‑100 (MQ=255)
                                                                                          tAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          <  1:3407666/100‑1 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑CT‑T‑TTC‑‑GTAGGCCTGATAAGAGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  2:1497602/1‑99 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑C‑TTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCATTtttgcctttgtcct         >  2:2211785/1‑96 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:3377901/1‑100 (MQ=255)
                                                                                            aCGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  2:1957047/100‑1 (MQ=255)
                                                                                             cGCCTTATCCGGCCT‑‑AC‑AAGATCGTGAAAATTCAATACATTGCA‑GC‑AT‑T‑TTC‑‑GTAGGCCGGATCAGCGGAACGCATTAGGCATTtttgcctttggta          >  2:3022143/1‑94 (MQ=17)
                                                                                             cGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  2:3341122/1‑99 (MQ=255)
                                                                                              gCCTTATCCGGCGT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg      <  2:1266239/100‑1 (MQ=255)
                                                                                                 ttATCCGGCTT‑‑AC‑AAAATCGTGCAAATTCAATAAATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCAATCAGTCt   <  2:1240889/100‑1 (MQ=255)
                                                                                                 ttATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt   <  1:1135717/100‑1 (MQ=255)
                                                                                                  tATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:2451083/1‑100 (MQ=255)
                                                                                                  |                                                                                                            
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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