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% 155.4 / NA 55T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (23/31);  minor base A (1/0);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG  >  NC_000913/3387799‑3388282
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
gATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggcg                                                                                                                                                                                                                                              <  1:17612/250‑1 (MQ=255)
       gcgGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAAt                                                                                                                                                                                                                                      <  2:184036/251‑1 (MQ=255)
          gTATAACCGGCCAGCCCCCCCGCATACGAGTTTTCTATTCGTACCAGCGAGGCTATCCAGCTACAAAAACCGGACCAGATACAGCACAAAAGAATCCTCTATAATGGAGCGCGGATCATCGCAATGATGATCCCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCAACGATAGAGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGACCCGGCGGGAACAATcgc                                                                                                                                                                                                                                   >  2:498744/1‑251 (MQ=255)
                 cGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                            >  1:103340/1‑251 (MQ=255)
                    ccAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGt                                                                                                                                                                                                                         <  2:317847/251‑1 (MQ=255)
                    ccAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGt                                                                                                                                                                                                                         <  2:120567/251‑1 (MQ=255)
                       gCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAcc                                                                                                                                                                                                                      <  1:234228/251‑1 (MQ=255)
                               gCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCg                                                                                                                                                                                                              >  1:107281/1‑251 (MQ=255)
                                              aTTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTg                                                                                                                                                                                               <  2:107281/251‑1 (MQ=255)
                                                                   cAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCaacaaaca                                                                                                                                                                           >  1:62204/1‑250 (MQ=255)
                                                                    aGGTACAAAAACCGGCCCAGATACAGCACCCCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgc                                                                                                                                                                         <  2:54048/251‑1 (MQ=255)
                                                                               ccGGCCCAGAAACAGCACACCAGAACCATCAATAATGGTGCGCGGATCATGGCAATGATGATCACCAGTCCGGCAATACAGCCAAAAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCCCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                              <  2:2831/251‑1 (MQ=255)
                                                                                 ggCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                            >  1:347498/1‑251 (MQ=255)
                                                                                 ggCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                            >  1:204381/1‑251 (MQ=255)
                                                                                          aCAGCACACCAGAATCATCAATAATGGTGCTCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTTCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGcc                                                                                                                                                    <  2:419253/250‑1 (MQ=255)
                                                                                                gacCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAACTGGGCCGTTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGTCCGGATCGGCGGCAACAATCGCCGCTGGCCGTACCGCCCAGCGTGTCGTTTCCAGCTGGAAGTGAACGCCAACAACCAGCGCCAGTGCATTGGCGGGCGCCAGttttt                                                                                                                                             >  2:133869/2‑250 (MQ=255)
                                                                                                 acCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATACGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTTAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGATTTAc                                                                                                                                            >  1:51686/1‑251 (MQ=255)
                                                                                                  ccAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAAGGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAcc                                                                                                                                           <  1:348678/251‑1 (MQ=255)
                                                                                                              gtaatGGGGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATGAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGGCAGAATACGGTTCTCCTCCCGGAGCAAAGGCCGGACCGGCGGCAACAATCGCCTCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAATCCAACATACAGCGCCAGGACAATGGCGGTCGCCAGTTTTACCGCAAAGCGa                                                                                                                                  <  1:300669/247‑1 (MQ=255)
                                                                                                                aatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCCGTACAATTGCGGTCGCCACTTTTACCCCAAAGCGActatt                                                                                                                            >  2:348757/1‑250 (MQ=255)
                                                                                                                aatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAatatt                                                                                                                            >  1:389294/1‑250 (MQ=255)
                                                                                                                aatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATg                                                                                                                            <  1:299311/251‑1 (MQ=255)
                                                                                                                   ggTGCGCGGATCATCTCAATGATGATAACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTg                                                                                                                         <  2:115067/251‑1 (MQ=255)
                                                                                                                    gTGCGCGGATCATCGTAATGATGATCCCCAGTCCGAAAATACAGCCAATAAAGGTGCTGAGGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCGGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGGTTCCAGCTGGAAGTGAAAGCCAAGAAACAGCGCCAGTACTATGGCGGTCGCCGGTGTTAGCGCAAAGTGAATATGTTgg                                                                                                                        <  2:435127/251‑1 (MQ=255)
                                                                                                                         cgGATCATCGCAATGATGATCCCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGCAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTATTACAGCCAAGCCCATATGTTGGTTAGc                                                                                                                   >  2:114978/1‑251 (MQ=255)
                                                                                                                            atcatcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAAt                                                                                                                <  1:414426/251‑1 (MQ=255)
                                                                                                                            atcatcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGATTTACCGCAAAACGAATATGTTGGTTAGCaa                                                                                                                 >  2:187882/1‑250 (MQ=255)
                                                                                                                              catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATg                                                                                                               <  1:361318/250‑1 (MQ=255)
                                                                                                                                    aaTGAGGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGaaaa                                                                                                       <  2:330085/251‑1 (MQ=255)
                                                                                                                                     atgatgATCACCAGTCCGGCCATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGttt                                                                                                                                               >  1:345473/1‑212 (MQ=255)
                                                                                                                                     atgatgATCACCAGTCCGGCAAGACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGttt                                                                                                                                               <  2:345473/212‑1 (MQ=255)
                                                                                                                                                cAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCCAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCCTCGTgg                                                                                           >  1:107251/1‑250 (MQ=255)
                                                                                                                                                     cGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAAc                                                                                      <  2:259616/251‑1 (MQ=255)
                                                                                                                                                                              gatgatGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGCTTTACCGCAAAGCGAATATGTTGGGTAGCAATGGGGAAAATACCCCTCGGGATTAACCAAACTCACGCAGGCgcgggggca                                                              >  1:284027/1‑242 (MQ=255)
                                                                                                                                                                               atgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCt                                                                                                                                                                                                >  1:328601/1‑121 (MQ=255)
                                                                                                                                                                               atgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCt                                                                                                                                                                                                <  2:328601/121‑1 (MQ=255)
                                                                                                                                                                                 gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAttt                                                           <  1:164501/250‑1 (MQ=255)
                                                                                                                                                                                    gcgcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTAc                                                        <  1:56819/250‑1 (MQ=255)
                                                                                                                                                                                              ccACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCACAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTTCCGCCAAACGAAAATGTTGGGTTGCAATGGGGAAAATAACCATCGTGATTAACCAAACTCACCCAAGCGCTGGGGCCTTTTATGggagaaaga                                              >  2:258142/1‑247 (MQ=255)
                                                                                                                                                                                                  gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGTAAGAACGAATCCt                                         <  2:333938/251‑1 (MQ=255)
                                                                                                                                                                                                   aTAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTg                                        <  2:276852/251‑1 (MQ=255)
                                                                                                                                                                                                         aaTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                  <  2:284027/251‑1 (MQ=255)
                                                                                                                                                                                                          aTAGCGACACAACACGGTTCACCGCCCTCAGCACAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTTGCGTTTCCAGCTCGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCCCCAGTTTTCCCGCAAAGCTAATAAGTTAGTTAGCACTGCAGAAAATACCCACCGTGATTAACCAAGCTCACGCACGCGATGGCCCATTTTAAGGAAGAACCAATCCTTGCtgtct                                 <  2:38010/251‑5 (MQ=255)
                                                                                                                                                                                                            aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgt                                <  1:92561/250‑1 (MQ=255)
                                                                                                                                                                                                                ccAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGACAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGTTTTTCCGCAAAGCGAATATGGTGGTTTGCCATGGCGGAAATACCCCTCGGGAGTTACCAAACCCACGCCGACGGTGGGCACTTTTTCGGACAAACGAATACTGGCttttctcgcg                            >  2:187796/1‑243 (MQ=255)
                                                                                                                                                                                                                       aCGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAGCAATTGAGAAAATAAACATCGTGATTAACCCAACTCACGCAGGCGGTGGGCCCTTTTACGGAAGAAAGGAACCCGGGATTCGTCGCGATCTTGtt                    >  2:314239/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATgg                                                                                                              <  1:441616/160‑1 (MQ=255)
                                                                                                                                                                                                                             cACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCa               <  1:325758/250‑1 (MQ=255)
                                                                                                                                                                                                                               ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTg            <  2:209669/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cccGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGacc         <  1:187882/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gcagcaAAGGCCGGACCGGAGGCAACAATCGCCGCTGTCAGTACCGCACAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCTACAAACAGCGCCAGTACAATGGCGGTGGCCAGGTTTTCCGACAAGCGGATATGTTGGGTAGCAATGGAGAACAAAACCCACGTGATTAACCAAACACAAGCAGGCGATAGGGCATTTTTCGGAAGAACGACACCCGGTTTTCGTCGCGATCTTGTTTAGCAGTGaccacc      >  2:35237/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCGCAAAGCGAAAATGTTGGTTTGCAATGGAGAAAATACCCATCGTGGTTAACCAAACTCACGCCGGCGGTGGGCCAATTTAAGGGAGAACGACACCTGGCTTTCGGCGCCATCTTGGTTTCCAGTGaccaccccc   >  2:185461/1‑248 (MQ=255)
                                                                                                                                                                                                                                        gcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGGAATGGAGGAAATACACATCGGGATTAACCAAAATCCAGCCGGCGCTGGGCCATTTTTCGAAAGAACGAATCCCGGGTTTCCTCGCGATCTTTTTTGCCCGGCACCACCACt   >  2:369312/1‑251 (MQ=255)
                                                                                                                                                                                                                                         caAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCCGTACAATAGCAGTCGCCAGATTTACCGCAAACGCCATATGTTGGTTAGCAATGGCGAAAATCCCCATCCTGATTAAACCAACTCACGCAGCCGAAGGGCCATTTTACGGAAGAAAGACTCCTCGCTTTCGGCGCGGTCTTGTTTGCCAATGaccccccctg  >  1:243359/1‑244 (MQ=255)
                                                                                                                                                                                                                                         caAAGGCCGGACCGGCGGCAACAATCGCCGATGTCAGTACCGCCCAGAGCGGCGTTTCCAGCTGGAAGTGCAAGCCAACAAACAGCGCCAGTACAATGGCGATCGCCAGTTTTACCGCCAAGCGAATATGTTGGTTAGAAATGGAGAAAATACCACTCGTGATTAACCAAACTCACGCAGGAGATGGGCCATTTTCCGGAAGAAAGAATCCTTGCTATTGTCGCGATCTTGTTTGCCAGTGAccaccaatg  >  1:209404/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
GATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG  >  NC_000913/3387799‑3388282

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