Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,048 G→T 28.2% intergenic (‑258/+50) 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,706,0480GT28.2% 601.4 / 228.5 313intergenic (‑258/+50)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (76/141);  new base T (64/18);  total (153/160)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 321 total aligned reads displayed.

TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147
                                                                                                         |                                                                                                   
tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCgg                                                                                                     <  1:3476098/100‑1 (MQ=255)
 gcgcTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    >  2:2524683/1‑100 (MQ=255)
 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGa                                                                                                    <  1:542938/100‑1 (MQ=255)
 gcgcTTGGTGGATGCCGTAGGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GCAT‑‑TTTC‑‑GTAGGTATGATAAGCgg                                                                                                     <  2:3865933/99‑1 (MQ=17)
  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAg                                                                                                   <  2:685038/100‑1 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  >  1:3551034/1‑100 (MQ=255)
   gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                                                                                                  <  1:1844845/100‑1 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:2337018/1‑100 (MQ=255)
    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                                                                                                 >  1:2234460/1‑100 (MQ=255)
     gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAgcgc                                                                                                >  1:28312/10‑100 (MQ=18)
      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                                                               <  1:635910/100‑1 (MQ=255)
      atttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCA‑GAAA‑‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCa                                                                                               >  1:253299/9‑100 (MQ=18)
       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCAt                                                                                              >  1:3894808/8‑100 (MQ=18)
         tttATGCCGGATGCGTCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            <  1:4070090/100‑1 (MQ=255)
         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                                                            >  1:3654130/1‑100 (MQ=255)
         tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                                                            >  2:586815/6‑100 (MQ=18)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGTG‑GGCT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:452644/5‑100 (MQ=21)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AGAT‑‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                                                           >  2:2685580/5‑100 (MQ=17)
          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                                                                                           >  1:2081561/5‑100 (MQ=21)
          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  2:368802/5‑100 (MQ=18)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  1:92434/100‑1 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                                                          <  2:2747110/100‑1 (MQ=255)
           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTAC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                                                                                           >  2:2564800/1‑99 (MQ=255)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                                                         >  2:2992827/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGA‑‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc                                                                                         >  1:3151707/3‑100 (MQ=18)
            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                         >  1:3217352/3‑100 (MQ=18)
             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                         >  1:1037291/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:1277254/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                                                                                        >  1:1462990/2‑100 (MQ=21)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑AAT‑‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  1:3000158/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCA‑GGAT‑‑ACGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                        >  2:4063685/2‑100 (MQ=18)
             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                                                       >  2:2929945/2‑100 (MQ=21)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:2663455/1‑99 (MQ=18)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TGAG‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  2:1063557/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                                                       >  1:1226572/1‑99 (MQ=17)
              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag                                                                                       >  2:383007/1‑99 (MQ=18)
              gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                                                       >  1:2356037/1‑100 (MQ=18)
                tggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                                                                                     <  1:3426685/93‑4 (MQ=25)
                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                                                                                     >  2:3637288/8‑100 (MQ=17)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:2875162/7‑100 (MQ=18)
                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GAGT‑‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  1:2021854/7‑100 (MQ=17)
                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                                                    >  2:875970/7‑100 (MQ=18)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                                                                                    >  2:3104334/7‑96 (MQ=25)
                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                                                   >  2:2104867/7‑98 (MQ=21)
                 ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                                                                                    >  1:878592/1‑100 (MQ=255)
                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                                                                                   >  1:3978452/6‑100 (MQ=18)
                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                                                  <  1:2634783/100‑1 (MQ=255)
                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                                                                <  1:694945/100‑1 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCACTACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATGcttgcct                                                                               >  2:4120221/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:39669/1‑100 (MQ=255)
                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                                                                               >  1:959821/1‑100 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            <  1:3804816/100‑1 (MQ=255)
                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                                                            >  2:1152345/1‑100 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:1509777/100‑1 (MQ=255)
                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc                                                                          <  1:794578/100‑1 (MQ=255)
                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATctg                                                                     >  2:2496787/1‑98 (MQ=255)
                                    tATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt                                                                 >  1:3398643/1‑100 (MQ=255)
                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  2:800362/100‑1 (MQ=255)
                                     aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                                                <  1:2306532/100‑1 (MQ=255)
                                         ggCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                                                            <  2:2202202/100‑1 (MQ=255)
                                           ccTAC‑AAAATCGTGCTAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                                                           <  2:2418874/99‑1 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:1336299/100‑1 (MQ=255)
                                           ccTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGc                                                          <  2:4040525/100‑1 (MQ=255)
                                               c‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                                                      <  1:2601241/100‑1 (MQ=255)
                                                      cGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTTATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                                                <  1:1184777/100‑1 (MQ=255)
                                                      cGTGAAAATGCAATAAATTGCA‑TCAT‑‑CTTC‑‑TTAGCCTTGTTAAGCGGAGCTCATCAGGCATTTTTGCCTTTCCCATCAGTCCTGTATGGCTTTTAatt                                                  <  2:1721596/98‑1 (MQ=255)
                                                       gTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt                                               <  2:2337005/100‑1 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  2:1454734/1‑100 (MQ=255)
                                                        tGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc                                              >  1:1909469/1‑100 (MQ=255)
                                                         gCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCg                                             <  2:3551011/100‑1 (MQ=255)
                                                          cAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa                                            <  2:1530810/100‑1 (MQ=255)
                                                            aaTTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAtn                                          <  1:3333189/100‑2 (MQ=255)
                                                             aTTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag                                         <  1:1107641/100‑1 (MQ=255)
                                                              ttCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt                                           <  2:1251386/97‑1 (MQ=255)
                                                                  aTACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg                                    <  2:1852365/100‑1 (MQ=255)
                                                                     cATTGCA‑GACT‑‑TTTC‑‑GTCGGCCTGATAAGCGGCGCGCATCAGTCATTTTTGCCTTTGCCATCAGTCCTGTACGGCCTTTAATTATTCGATAGAGATGttt                                 >  2:1966019/1‑100 (MQ=255)
                                                                       ttGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTc                               >  1:3267176/1‑100 (MQ=255)
                                                                        tGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCg                              <  2:3223579/100‑1 (MQ=255)
                                                                           a‑GAAT‑‑TTTC‑‑GTAGGCCTGTTAAGCGGAGCGCATCAGGCATTTTTGGCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAg                           <  2:1796508/100‑1 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  2:1766039/1‑100 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          <  1:1939843/100‑1 (MQ=255)
                                                                             gAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGt                          >  1:2952919/1‑100 (MQ=255)
                                                                                t‑‑tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGAt                       <  1:1424498/100‑1 (MQ=255)
                                                                                   tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaagttgcaagattttgtagtccggataaggcgttaacgcc                         <  2:1032666‑M2/97‑65 (MQ=255)
                                                                                   tttC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATg                      <  2:3398623/100‑1 (MQ=255)
                                                                                     tC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGtt                    <  2:2619780/100‑1 (MQ=255)
                                                                                           aGGCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggc                <  2:1208489‑M2/100‑74 (MQ=255)
                                                                                            ggCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggca               <  1:3105598‑M2/100‑75 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCctacgaaaattctgcaatgtattgaatttgcacgattttgtaggccggataaggcgttaacgccgcatccggca               >  2:1433885‑M2/1‑26 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATATAGACGTTTTCGAAGTGATGTTTGCCt               >  2:519829/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               >  2:2270031/1‑100 (MQ=255)
                                                                                            ggCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCt               <  1:4091775/100‑1 (MQ=255)
                                                                                             gCCTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaa              <  1:1383917‑M2/100‑76 (MQ=255)
                                                                                              ccTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgattttgtaggccggataaggcgttcacgccgcatccggcaag             <  2:919671‑M2/100‑77 (MQ=255)
                                                                                              ccTGATAAGCGTAGCGCATCAGGCctacgaaaatcttgcaatgtattgaattttcatgatttggtaggccggataaggcgttcacgccgcatccggcaag             <  2:3292919‑M2/100‑77 (MQ=255)
                                                                                                tGATAAGCGTAGCGCATCAGGCctacatggtcctgcaatatattgaattggcaagatttttgtaggccggataaggcgttcacgccgcatccggcatgaa           <  1:1615250‑M2/100‑79 (MQ=255)
                                                                                                tGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATg           <  2:1626204/100‑1 (MQ=255)
                                                                                                   tAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGAt        <  1:2785223/100‑1 (MQ=255)
                                                                                                     aGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCa      >  1:3369503/1‑100 (MQ=255)
                                                                                                        ggAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAAcc   <  1:1795122/100‑1 (MQ=255)
                                                                                                         |                                                                                                   
TGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTTTCGAAGTGATGTTTGCCTAATGGATCAACCA  >  NC_000913/3705949‑3706147

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