Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 776,311 G→A intergenic (+34/‑31) tolR → / → tolA Tol‑Pal system protein TolR/Tol‑Pal system protein TolA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913776,3110GA98.1% 179.0 / ‑3.9 52intergenic (+34/‑31)tolR/tolATol‑Pal system protein TolR/Tol‑Pal system protein TolA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (22/29);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGGTAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGTAGTTGAGCAGTACAAACGCATGCAAAGCC  >  NC_000913/776065‑776549
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
cAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAAc                                                                                                                                                                                                                                            <  1:314721/251‑1 (MQ=255)
     cacCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCg                                                                                                                                                                                                                                       <  2:367311/251‑1 (MQ=255)
      acCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCg                                                                                                                                                                                                                                       >  2:283476/1‑250 (MQ=255)
       ccGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGaa                                                                                                                                                                                                                                     <  1:29646/251‑1 (MQ=255)
        cGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAAc                                                                                                                                                                                                                                    >  1:77831/1‑251 (MQ=255)
         gtggtgGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACa                                                                                                                                                                                                                                   <  2:301680/251‑1 (MQ=255)
                      aaGATAGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAAcc                                                                                                                                                                                                                      >  1:134862/1‑251 (MQ=255)
                           cGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGagag                                                                                                                                                                                                                 <  2:44608/251‑1 (MQ=255)
                                 gAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCaa                                                                                                                                                                                                           <  2:113650/251‑1 (MQ=255)
                                     gTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGc                                                                                                                                                                                                       <  2:77831/251‑1 (MQ=255)
                                                        ggtggCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACaa                                                                                                                                                                                     <  1:112841/250‑1 (MQ=255)
                                                           ggCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTc                                                                                                                                                                                 >  1:352325/1‑251 (MQ=255)
                                                               gAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGc                                                                                                                                                                             <  1:337864/251‑1 (MQ=255)
                                                                               tCAAGGCTAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTAttt                                                                                                                                                              <  2:127172/250‑1 (MQ=255)
                                                                               tCAAGGCCAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATAGGTTGGTTTAATGACGAAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGCGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTc                                                                                                                                                             >  1:298038/1‑251 (MQ=255)
                                                                                     ccAACCCGGAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGATGCAGCCGATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAATCCGAGATTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGt                                                                                                                                                       <  2:432757/251‑1 (MQ=255)
                                                                                     ccAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGt                                                                                                                                                       <  1:303573/251‑1 (MQ=255)
                                                                                      cAACCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGtg                                                                                                                                                      <  1:341154/251‑1 (MQ=255)
                                                                                       aaCCCGAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCGTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGtgc                                                                                                                                                     <  1:283476/251‑1 (MQ=255)
                                                                                            gAAAACGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGACATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCAt                                                                                                                                                >  1:258167/1‑251 (MQ=255)
                                                                                              aaaCGATCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGt                                                                                                                                              <  1:129870/251‑1 (MQ=255)
                                                                                                acgatCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCa                                                                                                                                            <  1:440825/247‑1 (MQ=255)
                                                                                                           gATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTAttt                                                                                                                                    <  2:158282/248‑1 (MQ=255)
                                                                                                              cGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTgcggcg                                                                                                                              <  2:211422/251‑1 (MQ=255)
                                                                                                                                 ccTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGttcgttcg                                                                                                           <  2:383168/251‑1 (MQ=255)
                                                                                                                                  cTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTAATCTTATTTGCGGCGCTGATCTGGAGttcgttcgt                                                                                                          >  2:163624/1‑250 (MQ=255)
                                                                                                                                       gATGAAATAATTCAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATg                                                                                                        >  2:433777/1‑248 (MQ=255)
                                                                                                                                          gAAATAAGTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAata                                                                                                  <  2:91857/251‑1 (MQ=255)
                                                                                                                                             ataatTCAAGAACTGAACCTGTTACATAGTGCGGGCGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGATTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGAAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAatat                                                                                                 >  2:179558/1‑249 (MQ=255)
                                                                                                                                                atTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCTGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAACCGACAAGCTCAAGCGGGCGATAATTATTTCAGTTGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGacg                                                                                            <  1:203620/251‑3 (MQ=255)
                                                                                                                                                           gAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGa                                                                                 <  1:291041/251‑1 (MQ=255)
                                                                                                                                                             aCTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGACAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAgg                                                                               >  2:133098/1‑251 (MQ=255)
                                                                                                                                                              cTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAggc                                                                              >  2:394430/1‑251 (MQ=255)
                                                                                                                                                                tGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCATCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAggc                                                                              <  2:357193/249‑1 (MQ=255)
                                                                                                                                                                     cATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCggtg                                                                       >  1:135414/1‑251 (MQ=255)
                                                                                                                                                                      aTAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCggtgg                                                                      >  2:5903/1‑251 (MQ=255)
                                                                                                                                                                        aGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGtt                                                                    >  1:120431/1‑251 (MQ=255)
                                                                                                                                                                        aGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGtt                                                                    >  1:270062/1‑251 (MQ=255)
                                                                                                                                                                         gTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTc                                                                   >  1:254459/1‑251 (MQ=255)
                                                                                                                                                                          tacggtggtgAAATCGGTTGGTTTAATGACGCAGCCTATATAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAGCAAAACGACAAGGTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATATTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCg                                                                  >  2:30369/8‑251 (MQ=255)
                                                                                                                                                                           gCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGt                                                                 <  2:252651/251‑1 (MQ=255)
                                                                                                                                                                                   gAAATCGGTTTGTTTAATGACGCAGGCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAATCGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGAc                                                         <  1:430647/251‑1 (MQ=255)
                                                                                                                                                                                      aTCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAGCAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCTATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCt                                                      <  1:411050/251‑1 (MQ=255)
                                                                                                                                                                                             ggttTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATgg                                               >  1:216918/1‑251 (MQ=255)
                                                                                                                                                                                                 taacgacgtacccTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGAAGAACAGTTTTTGGTAACCGCGAGTGTCACAGGCATCCGACCAAAACGACAAGCTCAAGCGGGCGTTTCTCATTTCAGCCGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAGGCTTCAGCAGGAGGCGGCGTGGGTTTGTCCATCGACGCTGTCATGGTTg                                            >  2:203023/12‑250 (MQ=255)
                                                                                                                                                                                                  aaTGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGa                                           <  2:451478/250‑1 (MQ=255)
                                                                                                                                                                                                        gCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGt                                    <  2:276856/251‑1 (MQ=255)
                                                                                                                                                                                                             cTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGATCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGt                               >  1:152959/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tgcttgAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGTAGTTGAGCAGTACAAACGCATg         <  2:307460/249‑1 (MQ=255)
                                                                                                                                                                                                                                       cttgAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGTAGTTGAGCAGTACAAACGAATGCaa      >  1:168896/1‑250 (MQ=255)
                                                                                                                                                                                                                                          gAAAGAGAGCGGATAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGTAGTTGAGCAGTACAAACGCATGCAAAGcc  >  2:107862/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gCGGGTAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTc                                        >  1:263633/1‑205 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
CAGTACACCGTGGTGGTTGAGAAAGATCGCCTGGAGCGTTTACCACCAGAGCAGGTGGTGGCGGAAGTGTCCAGCCGTTTCAAGGCCAACCCGAAAACGGTCTTTCTGATCGGTGGCGCAAAAGATGTGCCTTACGATGAAATAATTAAAGCACTGAACTTGTTACATAGTGCGGGTGTGAAATCGGTTGGTTTAATGACGCAGCCTATCTAAACATCTGCGTTTCCCTTGCTTGAAAGAGAGCGGGTAACAGGCGAACAGTTTTTGGAAACCGAGAGTGTCAAAGGCAACCGAACAAAACGACAAGCTCAAGCGGGCGATAATTATTTCAGCAGTGCTGCATGTCATCTTATTTGCGGCGCTGATCTGGAGTTCGTTCGATGAGAATATAGAAGCTTCAGCCGGAGGCGGCGGTGGTTCGTCCATCGACGCTGTCATGGTTGATTCAGGTGCGGTAGTTGAGCAGTACAAACGCATGCAAAGCC  >  NC_000913/776065‑776549

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