Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,771,979 G→T intergenic (‑235/‑302) yibI ← / → mtlA DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,771,9790GT92.9% 125.3 / 5.0 42intergenic (‑235/‑302)yibI/mtlADUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (21/18);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

TAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCAC  >  NC_000913/3771736‑3772219
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
tAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTcaac                                                                                                                                                                                                                                           >  2:131915/1‑251 (MQ=255)
 aGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTcaac                                                                                                                                                                                                                                           <  1:124312/250‑1 (MQ=255)
    aaCATATCTGCCTCGGATTCCCGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGCAGGGCGTCAGGCTGATATTTTGACAATCATCCGGGTTCGCGATTCTCGCCATAACACCCACGAATACTTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCCTAACATGCTGTAGATCACATCAGGTGAACACCGTAAGACAATATCTTGTGATTCTGCTCACACATATTCAACAAAc                                                                                                                                                                                                                                       >  2:63634/1‑251 (MQ=255)
              ccgcggactcccGTTTAACAGTGTTTTTTTTGGGCGGGCAGCCAGAAGGGAGACAGGCTGATATGTTGACAATAAGCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTTATTCAGATCACAAATATTCAACAAGCCATcaaaaca                                                                                                                                                                                                                             <  2:32513/240‑1 (MQ=255)
                       cACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAAtgtg                                                                                                                                                                                                                     >  1:383599/1‑250 (MQ=255)
                                       gttTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACACTAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGTTTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGCACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAAAAAACATGTGACACTACTCACATTTa                                                                                                                                                                                                     >  2:475053/1‑250 (MQ=255)
                                                  ggCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTcaaca                                                                                                                                                                                                                                          <  2:368654/202‑1 (MQ=255)
                                                  ggCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTcaaca                                                                                                                                                                                                                                          >  1:368654/1‑202 (MQ=255)
                                                                                   cAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTacac                                                                                                                                                        >  1:470875/1‑251 (MQ=255)
                                                                                       aatCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGc                                                                                                                                                    >  2:97977/1‑251 (MQ=255)
                                                                                                tCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTtgtg                                                                                                                                            <  1:123345/250‑1 (MQ=255)
                                                                                                       tctcGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacg                                                                                                                                     <  2:70925/250‑1 (MQ=255)
                                                                                                          cGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAacgtcac                                                                                                                                 <  1:558144/251‑1 (MQ=255)
                                                                                                                     cAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACct                                                                                                                      <  2:92251/251‑1 (MQ=255)
                                                                                                                      aaaGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctc                                                                                                                     <  1:295255/251‑1 (MQ=255)
                                                                                                                       aaGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctct                                                                                                                    >  1:406716/1‑251 (MQ=255)
                                                                                                                                    agagGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTatat                                                                                                       >  2:132288/1‑251 (MQ=255)
                                                                                                                                       ggTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATAttt                                                                                                     >  2:233041/1‑250 (MQ=255)
                                                                                                                                          gatgatTTTCTTAGTTTCATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTAtgt                                                                                                 <  1:13649/251‑1 (MQ=255)
                                                                                                                                                  gCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGtgattgat                                                                                          >  1:289683/1‑250 (MQ=255)
                                                                                                                                                         tACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCATATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgttg                                                                                                                                                                   >  2:88673/1‑170 (MQ=255)
                                                                                                                                                         tACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCATATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgttg                                                                                                                                                                   <  1:88673/170‑1 (MQ=255)
                                                                                                                                                               aCGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGc                                                                            >  2:225934/1‑251 (MQ=255)
                                                                                                                                                                    tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTc                                                                        <  2:269438/250‑1 (MQ=255)
                                                                                                                                                                    tGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCg                                                                       >  1:283993/1‑251 (MQ=255)
                                                                                                                                                                       aTGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACt                                                                    <  2:180246/251‑1 (MQ=255)
                                                                                                                                                                       aTGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATAGCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATGTATGTGATTGATATCACACAAAAGGCCGTCGACt                                                                    <  2:383599/251‑1 (MQ=255)
                                                                                                                                                                          aCGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                           <  2:116678/209‑1 (MQ=255)
                                                                                                                                                                          aCGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                           >  1:116678/1‑209 (MQ=255)
                                                                                                                                                                          aCGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTgg                                                                  >  1:436707/1‑250 (MQ=255)
                                                                                                                                                                            gAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAATGCCGTCGACTGGAc                                                                >  2:554337/1‑250 (MQ=255)
                                                                                                                                                                              aGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACactact                                                                                                                                                                                                             <  2:20110/107‑1 (MQ=255)
                                                                                                                                                                              aGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACactact                                                                                                                                                                                                             >  1:20110/1‑107 (MQ=255)
                                                                                                                                                                                   tAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAAcc                                                        >  1:38104/1‑251 (MQ=255)
                                                                                                                                                                                    aaCATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAAcc                                                        <  1:131915/250‑1 (MQ=255)
                                                                                                                                                                                       aTGCTGTAGAACACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGAGCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTGTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTATCCGAtt                                                    <  2:208929/251‑1 (MQ=255)
                                                                                                                                                                                                  cacaTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGcc                                             >  1:204850/1‑247 (MQ=255)
                                                                                                                                                                                                  cacaTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGAt                                         <  1:255753/251‑1 (MQ=255)
                                                                                                                                                                                                         ggTGAACGCCGTAAGAAAATATCTTGTGATTCAGAGCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATGTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAg                                  <  1:57569/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGCCAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACa                          >  2:147277/1‑250 (MQ=255)
                                                                                                                                                                                                                    tAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAgg                       >  1:324643/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGACCGGCTACCTCTGCCGCcac  >  1:39867/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                
TAGAAACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTGCTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCAC  >  NC_000913/3771736‑3772219

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: