Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,388,041 T→G T50P (ACG→CCG)  aaeB ← aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,0410TG98.2% 163.6 / ‑6.0 55T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (29/25);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGT  >  NC_000913/3387797‑3388283
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
aTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggc                                                                                                                                                                                                                                               >  1:89396/1‑251 (MQ=255)
aTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggc                                                                                                                                                                                                                                               <  1:107050/251‑1 (MQ=255)
    tCAGCGCGGTATAACAGGCTAGCACCCACGCATCCGAGTTTTCTATTAGTAACCGAGAGGATATACAGGTACAAACACAGGCCCAGACACAGCACACCAGAATCCTCAATAATGGTGCGCGGAACATCGAAATGATGATCAACAGTCCGGCAATACAGCCAATAAATGTGCCGATGACGCGCAAACAGCCACGAACGCGAAAAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCggaggc                                                                                                                                                                                                                                            >  2:18215/1‑250 (MQ=255)
       gcgcGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTCCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAACAAAGGCCGGACCGGCGGcaaca                                                                                                                                                                                                                                        <  1:296160/251‑1 (MQ=255)
                aaCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCTCGATAGCGAATGGCGGCAGAATACGGTTCGCCTCCCGCAGCGAAGGCCGGACCGGCGGCAACAATCGCCGCt                                                                                                                                                                                                                               <  2:256404/251‑1 (MQ=255)
                         gCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAcc                                                                                                                                                                                                                      <  2:426915/251‑1 (MQ=255)
                            cccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGcc                                                                                                                                                                                                                   >  1:391764/1‑251 (MQ=255)
                                                 ttCGTACAAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTgg                                                                                                                                                                                              >  1:323576/1‑251 (MQ=255)
                                                      aCCAGAGAGGATATCCAGGAACAAAAACCGGCCCACATCCCGCCCACCAGCATCATCAATAATGGTGCGCGGATCATCACAATGATGCTCACCAGTCCGGCACTACCGCCAATACATCTGCCGATGATGCGCAAAAAGCCACGCTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGCTTCCCGCTGGCAGTg                                                                                                                                                                                         >  2:54694/1‑251 (MQ=255)
                                                        cAGCGAGGATACCCAGGTACAAAAACCGGCCAAGATACAGCCCACCCGAATCATCAATCATGGTGCGCGGCTCATCGCCATGATGATCACCAGTCCGGCAATACAGCCAATAACTGTGCCGCTGATGCGCCAAAAGCCACGCTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCACAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGACGTGa                                                                                                                                                                                        >  2:61198/1‑250 (MQ=255)
                                                            gAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGAGGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                                   <  2:398948/251‑1 (MQ=255)
                                                                     cAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCaacaaac                                                                                                                                                                            >  1:76055/1‑249 (MQ=255)
                                                                        gTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcgc                                                                                                                                                                       >  1:395747/1‑251 (MQ=255)
                                                                            aaaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGt                                                                                                                                                                   <  2:337949/251‑1 (MQ=255)
                                                                             aaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGt                                                                                                                                                                   >  2:193305/1‑250 (MQ=255)
                                                                              aaaacGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCCGTACCGCCCAGCGCGCCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTAc                                                                                                                                                                 >  2:229631/5‑251 (MQ=255)
                                                                                 ccGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                              >  2:313851/1‑251 (MQ=255)
                                                                                  cGGCCCAGATACAGCACACCAGAATCATCAATAATGTTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAAGACGGTTCACCTCCCGCAGCAAAGGCCGGAGCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTGTTCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATg                                                                                                                                                             <  1:341275/251‑1 (MQ=255)
                                                                                   ggCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                            <  2:90826/251‑1 (MQ=255)
                                                                                          aTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACGGCCATTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACATACCGCAGCAAAGTCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCGGCGTGGCGTTTACGGCTGGGAGAGAAAGCCAACAAACGGTGCCAGTACAGTGGCggaagg                                                                                                                                                     >  1:218938/1‑247 (MQ=255)
                                                                                          aTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCAGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTAACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAACCGACAAAAAGCGCCAGTACAATGGCGGgcgc                                                                                                                                                     >  1:296747/1‑251 (MQ=255)
                                                                                                  cacCAGAATCATCAAGAAGGGTGCGCGGATCATCGCAATTATGATCACCAGGCCGGCAATACAGCCAATAAATGTGCCGAGGATGCGCAAAAAGCCACGGTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTTCCGCCCAGCGTGGCGTGTCCATCTCGAAGTGAAAGCCAACAAACAGCGCCATTACAATGGCGTTCGCCAGTTTTa                                                                                                                                             <  2:168387/251‑1 (MQ=255)
                                                                                                          tcatcaATAATAGCGCGCGGATCATCGCAATGATGATCACCAGCCCGGAAATACAGCCAATAACTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATCGCGCCAGAATACGGTTCACCTCCCGCAGCAACGGCCGGACCGGCGGCAAAAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGTAGTGCAAGCCAACAAACAGCACCAGCACAATGGCgggcccccgctttacccacaa                                                                                                                                      >  2:25403/1‑231 (MQ=255)
                                                                                                           catcaATAATGGTGCGCGGATCATCGCAATGAAGATCACCAGTCCGGCAATACAGCCAATAAATGAGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCCGGCAGTGAACGCCAACAAACAGCGCCAATACCATGGCGGACGcccagttaaacgcaaaac                                                                                                                                    >  1:42978/1‑235 (MQ=255)
                                                                                                              caATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTACCGCaaaacgg                                                                                                                                  >  1:318625/1‑246 (MQ=255)
                                                                                                                       ggcgcGGATAATCGCAATGAGGAGCACCAGTCCAGCAATACAGCCAATAAAGGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCTCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGCATATGtttgt                                                                                                                        <  2:411015/250‑4 (MQ=255)
                                                                                                                                catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCTGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCTCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATgg                                                                                                               <  1:177743/251‑1 (MQ=255)
                                                                                                                                catcGCAATGACGATCACCTGTCCGGCAATACAGCCAATAAATGTGCCGACGATGAGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCCGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAACGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACACCGCCAGTACAATGGCGGTACCCAATTTTAACGCAAACCGAATACGTTGGTTAGCaaagg                                                                                                               >  2:179671/1‑248 (MQ=255)
                                                                                                                                    gCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGAcc                                                                                                                                                                                                                                                  <  2:285638/116‑1 (MQ=255)
                                                                                                                                    gCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGAcc                                                                                                                                                                                                                                                  >  1:285638/1‑116 (MQ=255)
                                                                                                                                          atgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAcc                                                                                                     <  1:124740/251‑1 (MQ=255)
                                                                                                                                          aggatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAcc                                                                                                     <  2:217164/249‑1 (MQ=255)
                                                                                                                                                  cAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGAGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTACAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACACTGGCGGTCGCCCGTTTTTCCCCCACACGCATATGTTGGTTAGGCATGGGGAAAAATCCCATcctgg                                                                                             >  2:186950/1‑247 (MQ=255)
                                                                                                                                                   aGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGAt                                                                                            >  1:127660/1‑251 (MQ=255)
                                                                                                                                                     cccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTa                                                                                          <  2:339951/250‑1 (MQ=255)
                                                                                                                                                        ggCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTAACGCAACACGAATATGTTGGGTAGGCATTGAGGAAATACCCCTCGTGATTaaac                                                                                       >  2:174419/1‑249 (MQ=255)
                                                                                                                                                            agaCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTGCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGGCGCCAGTTTTACCGCAAAGCGAATATGTTGGTGAGCAATGGAGAAAATACCCATCGTGATTGACCaatc                                                                                   <  2:145846/249‑3 (MQ=255)
                                                                                                                                                              aCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCACAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAtctc                                                                                 >  1:358656/1‑251 (MQ=255)
                                                                                                                                                                           gtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTGAGCAATGGAGAAAATTCCCATCGTGATTAACCAATCTCACGCAGGCGATgg                                                                    <  1:214430/251‑1 (MQ=255)
                                                                                                                                                                               cGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCa                                                               <  2:360658/251‑1 (MQ=255)
                                                                                                                                                                                   gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCAtttt                                                           <  1:220236/251‑1 (MQ=255)
                                                                                                                                                                                         cAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGgaa                                                     <  2:358656/251‑1 (MQ=255)
                                                                                                                                                                                         cAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGga                                                      >  1:301746/1‑250 (MQ=255)
                                                                                                                                                                                            aaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGgaagaa                                                  <  2:323416/251‑1 (MQ=255)
                                                                                                                                                                                                       aGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTTCCGCAAAGCGAATATGTTGGGTAGCAAAGGAGAAAATAACCATAGTGATTAACCAATCTCACACAGGCGGTGGGCCATTTTACGGAAGAAAGAATCCTGGc                                       >  2:297952/1‑251 (MQ=255)
                                                                                                                                                                                                          gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCCAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAACGCAAAGCGGATATGTTGGTTAGCAATGGGGAAAATACCCAACGGGATTTACCAATCTAAAGCCAGAGAAGGGCCATTTTACGGAAGAACGAATCCTGGCttt                                    >  1:296413/1‑251 (MQ=255)
                                                                                                                                                                                                          gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGGTGGTTAGCAATGGAGAAAATACCCATCGGGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCttt                                    <  1:37284/251‑1 (MQ=255)
                                                                                                                                                                                                           aaTAGCCCCAGAATTCGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAGCTCACGCAGGCGATGGGCCATTTTACTGAAGAACGAAGCCGGGCTTtc                                   <  2:42978/251‑1 (MQ=255)
                                                                                                                                                                                                              aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                                <  2:222823/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTcgc                              <  2:117116/251‑1 (MQ=255)
                                                                                                                                                                                                                cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTACAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTTCCGCAAAGCGAAAATGTTGGGTTGCAATGGTGAAAATACCCATCGTGATTAAACAATCTCACCCAGGCGGTGGGCCCTTTTTCGGAAGCACGAATCCTGGCTTTTGTcgc                              >  1:265273/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGGCAGGTTTACCCCAAAGCGAAAATGTTGGTTTGCCATGGAGAAAATACCCATCGTGATTAAACAATCTCACGCAGGGGATGGGGCCTTTTTCGGGAAGACGAATCCCGGCTTTCGTCGCGCTATTGTTTGTCAGTg             >  2:141204/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTTGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGAGGGGCCATTTTACGGAAAAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGGGac           >  1:294995/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACt    <  1:113259/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaaGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGGCCATTTTACGGAAGAAAGAATCCTGGCTTTCGTGTCGATCTTGTTTGcccgggaacacccccgt  >  1:413992/1‑236 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
ATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGT  >  NC_000913/3387797‑3388283

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