Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,978,083 +A coding (115/351 nt) flhD ← DNA‑binding transcriptional dual regulator FlhD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,978,0831.A95.3% 136.6 / 0.9 43coding (115/351 nt)flhDDNA‑binding transcriptional dual regulator FlhD
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (25/16);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01

TTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TT‑TATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACC  >  NC_000913/1977834‑1978324
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
ttttCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑tt                                                                                                                                                                                                                                                    <  1:441753/250‑1 (MQ=255)
    cACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑ttatat                                                                                                                                                                                                                                                <  1:506907/250‑5 (MQ=255)
     aCTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGc                                                                                                                                                                                                                                              >  1:634207/1‑251 (MQ=255)
         aTGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCg                                                                                                                                                                                                                                            >  1:638745/1‑249 (MQ=255)
           gttCAGGCCCTTTTCTTGCGCATCGCTTCTTCATGCTGATTATCATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGACGTGTCAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGTAAGAGTCAGTGCCGCTAACGTTGTCGCTATTTCTTCA‑TTATATGCCGAGAc                                                                                                                                                                                                                                        <  1:398516/249‑1 (MQ=255)
                                    cttcttCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTg                                                                                                                                                                                                               >  1:216578/1‑251 (MQ=255)
                                         tcAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACa                                                                                                                                                                                                           >  1:584316/1‑250 (MQ=255)
                                           aGGCTGATTAACATCATTCAGCCCACGTGTTGAGAGCATGATGCCGGTATGAATTTGATGGAGATCGACAACGCTGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAACACGGAAGTGACAAACAAGTTGATTGGTTTCTGCCCGCTTAACCATTTGCGGAAGAGTCCGTGCCGCTAACGTTGTCGCCCTTTCTTCA‑TTATATGCCGAGCCGAACCATAGCTGCCGCTTTGTCCTGAACAAtt                                                                                                                                                                                                        >  2:258058/1‑250 (MQ=255)
                                                  ttAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCt                                                                                                                                                                                                 >  1:56690/1‑251 (MQ=255)
                                                   tAACATCATTCAGCAAGCGTGTTGAGAGGATGATGGCGGTATGCATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGGCAACTTAGTAATCGGCTGGTGGCTGTCAAAACGGAAGGGACAAACCAGTTGAGTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTGCAGTA‑TATGCCTATACGAAACAGAGCGGACGTTTTGTCCTGAACAATCAAACGCt                                                                                                                                                                                                 <  2:697374/250‑1 (MQ=255)
                                                         catTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCGGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGGCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAag                                                                                                                                                                                          <  1:671493/251‑1 (MQ=255)
                                                                  aGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATAAAACGCTGTGCAAGAAGTAAat                                                                                                                                                                                  >  2:316558/1‑250 (MQ=255)
                                                                       gtggagaGCATGATGCCGGTATGTATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGAGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGAc                                                                                                                                                                             <  2:386294/247‑1 (MQ=255)
                                                                         tGAGAGCATGATGCCGGTAGGAATTTGCTGGAGATCGTCAACGCTGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACaa                                                                                                                                                                           <  2:476305/250‑1 (MQ=255)
                                                                                atgatgCCGGTATGAATTTGCGGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGATTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACTTTGTCGCCATTTCTTCA‑TTATATGCCGAGAGGAAACATAGCGGACGCTTTGTCCGGAACAATCAATCACTGTGCAAGTAGGAAATGTGACAAGGTGATGt                                                                                                                                                                   <  2:335476/251‑1 (MQ=255)
                                                                                atgatgCCGGTATGAATTTGATGGAGATCGTCAACGCGGGAATCTTGCGTCCACTGAGTAATCGTCTGGTGGCTGTCACAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCC‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTAAACAATCAAACGCTGGGCAAGTAGTAAAAATGCAAAGTTGATGt                                                                                                                                                                   >  2:490531/1‑251 (MQ=255)
                                                                                 tgatgCTGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACGGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGAt                                                                                                                                                                     <  1:660152/248‑1 (MQ=255)
                                                                                 tgatgCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGAAGTAAATATGAAAAGTTGaa                                                                                                                                                                     >  2:660152/1‑247 (MQ=255)
                                                                                     gCCGGTATGAATTTGCTGGAGATCGTCACCGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCAGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGCGACGAAACATAGCGGACGCTTTGTCATGAACAATCAAACGCTGTGCACGTAGTAAATATGACAAGTTGATGTCATa                                                                                                                                                               >  2:664055/1‑250 (MQ=255)
                                                                                      ccGGTATTAAGTTGCTGGAGAGCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGGCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCTGACGCTTTGTGCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATaaa                                                                                                                                                             <  2:158009/251‑1 (MQ=255)
                                                                                               aTTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTc                                                                                                                                                     >  2:174598/1‑250 (MQ=255)
                                                                                                   gCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCaa                                                                                                                                                <  1:109591/251‑1 (MQ=255)
                                                                                                    cTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGATGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTCACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGACACATAGCGGACGCTTTGTCCTGAACAATAAAACGCTGTGCAAGTCGTAAATATGACAATTTACTGACATAAATGTGTTTccacaac                                                                                                                                               >  2:467952/1‑245 (MQ=255)
                                                                                                                aaCGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAagta                                                                                                                                                                                        >  2:48450/1‑198 (MQ=255)
                                                                                                                aaCGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAagta                                                                                                                                                                                        <  1:48450/198‑1 (MQ=255)
                                                                                                                     gggAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCattat                                                                                                                              >  1:402525/1‑251 (MQ=255)
                                                                                                                           cTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGCCGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTcccac                                                                                                                        >  2:164986/1‑251 (MQ=255)
                                                                                                                              gCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTcccaccca                                                                                                                     <  1:566521/251‑1 (MQ=255)
                                                                                                                                gTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCC‑TT‑TATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGa                                                                                                                   >  2:307475/1‑250 (MQ=255)
                                                                                                                                  cAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCAGTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTTCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGaa                                                                                                                  <  1:295512/250‑1 (MQ=255)
                                                                                                                                   aaCTGAGTAATCGTCTGGTGGCTGTCAAAACGGACGTGACAAACCAGTTGATTGGTTTCTGCCAGCTATACCATTTGCGGAAGCGGCAGAGCCGCTAACGTTGTGCCCATTTCTTCA‑TTATATGCCGAGATGAAACATAGCGGTCGCTTTGTCCTGAGCAATCAAACACTGGTCAAGTAGTAACTATGACATGTTGATGTCACAAATGTGTTTCAGCACCTCGGAGGTATGCATTATTCCAGCCCAGAAt                                                                                                                 >  2:196760/1‑250 (MQ=255)
                                                                                                                                                    gtggCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGATTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTAt                                                                                               >  1:703782/1‑251 (MQ=255)
                                                                                                                                                          gTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGtt                                                                                         >  1:80805/1‑251 (MQ=255)
                                                                                                                                                            cAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGCGACGAAACATAGCGGACGCTTTGt                                                                                                                                                                                                                   >  1:642106/1‑127 (MQ=255)
                                                                                                                                                                aCGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCg                                                                                   >  1:668848/1‑251 (MQ=255)
                                                                                                                                                                                 gTTGATTGGTTTCTGCAAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCAATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGATGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCCCCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGAGATGTCGcc                                                                  >  2:572040/1‑251 (MQ=255)
                                                                                                                                                                                           ttCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTACTGAACAATCAAACGATGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccc                                                        >  1:632118/1‑251 (MQ=255)
                                                                                                                                                                                             cggCCAGCTTTACCATTTCCGGAAGAGTCAGTGCCGCTAACGTTCTCGCAATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAAAGCTGTGCAAGGAGTAAATATGACAAGTTGATGTCATAACTGTGTTTCAGCAACTCGGAGGTATGCCTTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTAGCCGCACCCTGTGATGACACCGGGTAGCCCCgg                                                      >  2:16226/3‑251 (MQ=255)
                                                                                                                                                                                               gCCAGCTTAACCAGTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCTCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTTTGCATTATTCCCACCCAGAATAACCAACTGTATTTTTATGCGGTTTCCCCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTa                                                    <  1:480500/251‑1 (MQ=255)
                                                                                                                                                                                                   gCTTAACCATTTGCGGAAGAGGCAGTGCCGCTAACGTTGTCGCCATTTCTTCC‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTaaaaa                                                <  1:258058/251‑1 (MQ=255)
                                                                                                                                                                                                      tAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCGCCCCGTGATGTCGCCGGGAAGCCCCGGTaaaaaa                                               <  1:382260/249‑1 (MQ=255)
                                                                                                                                                                                                                 ggAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCAtt                                    >  1:691103/1‑249 (MQ=255)
                                                                                                                                                                                                                             gCCGCTAACGTTGTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGAt                      <  2:56690/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gTCGCCATTTCTTCA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTTAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTTATGTCATAAATGTGTTTCAGCAACTCGTAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCTGTTTCACCGCACCCCGTGATGTCGCCGTGAATCCCCTGTAAAAAATAATTAGCATTAGAATAGTTGCTATAATCTGCAATa           <  2:319885/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ttcttcA‑TTATATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAacc  >  1:695490/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
TTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAGAGCATGATGCCGGTATGAATTTGCTGGAGATCGTCAACGCGGGAATCTTGCGTCAACTGAGTAATCGTCTGGTGGCTGTCAAAACGGAAGTGACAAACCAGTTGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCA‑TT‑TATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACC  >  NC_000913/1977834‑1978324

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