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.4% 197.8 / NA 61T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (26/34);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAG  >  NC_000913/3387800‑3388285
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
aTCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                             <  2:85316/251‑1 (MQ=255)
aTCAGCGCGGTAGAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCggcggc                                                                                                                                                                                                                                             <  2:261248/251‑1 (MQ=255)
    gcgcGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCAGcaaca                                                                                                                                                                                                                                         <  1:151062/251‑1 (MQ=255)
       cgGTATACCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATc                                                                                                                                                                                                                                      >  1:402761/1‑251 (MQ=255)
            aaaCCGGCCAGCCCCCACGCAGACGAGTTATCTATTCGAACCAGCGAGGATAGCCAGGTACAAAAACCGGCCCAGATACAGCACTCCAGAATCAGCAATAATGGTGCGCGGAGCATCGCAATGATGATCACCTGTCCGGCAAGAGAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAAGGCGGTTCA‑CCTCGCGCAGGAAAGGCCGGACCGGCGGCAACAATcgccgc                                                                                                                                                                                                                                 <  1:240399/250‑1 (MQ=255)
               ccGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCAGCAATAATTGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                                              <  1:451248/251‑1 (MQ=255)
                 ggCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                             >  1:87061/1‑250 (MQ=255)
                  gCCAGCCCCCACGCATACGAGTTTGCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGACCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGTCAACAATCGCCGCTGTCAg                                                                                                                                                                                                                           <  2:137123/251‑1 (MQ=255)
                           cACGCATTCGAGTTTACTATTCGTACCAGCTAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCTCGGATCATCGCAATGATGATCACCAGTCCGGCGATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATCCGGTTCT‑CTTCCCTCAGCAAATGCCGGACCGGCGTCAACAATCGCCGCTGTCAGTACAGCCCa                                                                                                                                                                                                                  <  2:173178/251‑1 (MQ=255)
                             cGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGc                                                                                                                                                                                                                >  1:5242/1‑251 (MQ=255)
                                  aCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTgg                                                                                                                                                                                                            <  1:187004/250‑1 (MQ=255)
                                                 gTACCAGCGAGTATATCCAGGTACAAAACCCGGCCTAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATTATCACCAGTCCGGCAATACAGCCCATAAATGTGCCTATGATGCTCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGTTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCGAGCGTGGCGTTTCCAGCTGGAAg                                                                                                                                                                                            <  1:481037/251‑1 (MQ=255)
                                                          aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCCCACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCca                                                                                                                                                                                   <  2:102831/251‑1 (MQ=255)
                                                          aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                                    >  1:430367/1‑250 (MQ=255)
                                                          aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCca                                                                                                                                                                                   <  1:53863/251‑1 (MQ=255)
                                                               atCCAGGTACAAAACCCGGCCCAGATACAGCACACCAGAATCCTCAATAATGGCGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCCATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaacaa                                                                                                                                                                               >  1:429483/1‑250 (MQ=255)
                                                                  cAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAg                                                                                                                                                                           <  1:210163/251‑1 (MQ=255)
                                                                           aaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTAc                                                                                                                                                                  <  1:16088/251‑1 (MQ=255)
                                                                                ggCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGCTCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGCTAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATgg                                                                                                                                                             >  2:91119/1‑251 (MQ=255)
                                                                                         aCAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGccc                                                                                                                                                    >  1:360185/1‑250 (MQ=255)
                                                                                          cAGCCCACCAGAAGCATCAATAATTGTGCGCGGCTCATCGCAATGAGGAGCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACTGGCGGCAACAAGCGCCGCGGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACATCGCCAGTACAATGGTGGTCGCCAg                                                                                                                                                   <  1:309295/251‑1 (MQ=255)
                                                                                                        catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGc                                                                                                                                     <  2:459393/251‑1 (MQ=255)
                                                                                                            aataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCggat                                                                                                                                 >  1:377981/1‑248 (MQ=255)
                                                                                                               aatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATg                                                                                                                              >  1:83626/1‑251 (MQ=255)
                                                                                                                             catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATgg                                                                                                                <  2:149929/251‑1 (MQ=255)
                                                                                                                             catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAAAAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCCAAGCGAATATGTTTGGTAGCACTgg                                                                                                                >  2:420958/1‑251 (MQ=255)
                                                                                                                                 gCAATGATGATCACCAGTCCGGCAATACAGCCCATAAATGTGACGATGATACGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCAGCTGTCAGTACCGCCCAGAGCGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTAGCCAGGTTTACCGCAAAGCGAATAAGTTGGTTAGCAAAGgaga                                                                                                             >  1:317409/1‑250 (MQ=255)
                                                                                                                                      gatgatAACCAGACCGGCAATACAGCCAATAAATGTGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCACGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTAAGTACCGCCCAGCGTGGCGTTTCCAGCTGGCAGTGAAAGCCAACAAACAGCGCAAGTACAATGGCGGTCTCCAGTTTTTCCGCAAACCGAATATATTGGGTATCAATGGAGAAAATa                                                                                                        >  2:468905/1‑250 (MQ=255)
                                                                                                                                       atgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGCTTTTTCGCAAAGCGAAAATGTTGGTTTGGAAAGGAGAAAAAAcc                                                                                                      >  2:129460/1‑251 (MQ=255)
                                                                                                                                          atCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAt                                                                                                   <  1:531002/251‑1 (MQ=255)
                                                                                                                                            cACCAGTCCGGCATTCCAGCCAATAAATGGGCCGATGATGCGCACAAAGCCACGAGAGCGAATAGCGCCAGAATGCGTTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGAGTCAACATTCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTTGAAGTGAAAGCCAACAAACAGCGCCAGTACAATTGCGGTCGCCAGGTTTACCGCAAAGCGAATATGTTGTTTAGCAATTGATAGAATACCCATc                                                                                                  <  2:328820/250‑1 (MQ=255)
                                                                                                                                            cACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCg                                                                                                 <  2:360185/251‑1 (MQ=255)
                                                                                                                                            cACCAGGCCGGCAATACAGCCAATAAATGTGCCGGTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGGAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAAtttccagc                                                                                                  <  2:339002/250‑8 (MQ=255)
                                                                                                                                              ccAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTg                                                                                               <  1:129460/251‑1 (MQ=255)
                                                                                                                                                      gCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGAAACAATCGCCGCTGTCAGTAACGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTAACGCAAAGCGAATATGTTGGTTAGCCATGGGGAAAATACCCATCgggaataaac                                                                                        >  1:100432/1‑241 (MQ=255)
                                                                                                                                                              gCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAAAACACATCGGGATTAAACAAACTCACg                                                                               >  1:339260/1‑251 (MQ=255)
                                                                                                                                                                  aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAgg                                                                           <  1:524777/251‑1 (MQ=255)
                                                                                                                                                                      atgcgcctatgatgCGCAAAAAGCCACGATAGCGAAGAGCGACATAAGACGTTTCA‑CCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGACAGTCCCGCCCAGCGTCGCGGTTCCAGCCGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCTAATATGTGGGGAAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGAt                                                                       <  2:259357/243‑1 (MQ=255)
                                                                                                                                                                       tgtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCCTTTTAACGCAAAGCGAACATGTTGGTTAGCAAAGGAGAAAATACCCACCTTGATTAACCAAAATCACCCCGGCGATg                                                                      >  2:470093/1‑251 (MQ=255)
                                                                                                                                                                        gtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCCATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATgg                                                                     >  1:420100/1‑251 (MQ=255)
                                                                                                                                                                          gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAAAATGTTGGGTAGAAATGGAGAAAATAACCATCGGGATTAAACCAACACACGCAGGCGATgggg                                                                   >  2:386067/1‑250 (MQ=255)
                                                                                                                                                                               tgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGGcc                                                                  <  1:271859/247‑1 (MQ=255)
                                                                                                                                                                               tgatgCGCAAAAAGCCACGATAGAGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCCGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGTTTTTCCCCCACACGCAAATGGTGGGTTGCAATGGAGGAAATACCCATAGGGGTTAACCCAACTCACGCAGGCGGGGGGCCaga                                                               >  2:271859/1‑248 (MQ=255)
                                                                                                                                                                                      cAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaa                                                       <  2:377981/251‑1 (MQ=255)
                                                                                                                                                                                          aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACGCACGCAGGCGATGGGCCATTTTACGGAAGAAc                                                   <  1:179624/251‑1 (MQ=255)
                                                                                                                                                                                                 gATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCCTTTTACGGAAGAACGAATCCt                                            >  1:198253/1‑251 (MQ=255)
                                                                                                                                                                                                    aGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGc                                         >  1:482942/1‑251 (MQ=255)
                                                                                                                                                                                                      cGAATATCGCCAGAATACGGTTCAACCTCC‑GCAGCAAAGGTCGGGCCGGCGGCAACAATCGCAGCTGTACGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCTAACAAACAGCGCCAGTACAAGTGCGGTCGCCAGCTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGTGATAATTCCCAGCGTGATTAACCAAACTCACGCAGGCGATGCGCCATGTTTCGGAAGTACGTATCCTgggtt                                       <  2:85530/251‑4 (MQ=255)
                                                                                                                                                                                                      cGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGGATATGTTGGTTAGCAATGGAGGAAAAAACCATCGTGATTAACCCAACTCCCGCAGGCGATGGGCCATTTTACGGAAGAACGAATccctgctt                                       >  2:145344/1‑245 (MQ=255)
                                                                                                                                                                                                      cGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCtt                                       <  2:90433/251‑1 (MQ=255)
                                                                                                                                                                                                       gAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATTTTGGTTAACAAAGGAGAAAATACCCACCGGGATTAACCAAACAAACGCCGGCGAGGGGGCATTTTACGGAAGAACGAAACCCGGCttt                                      >  2:209521/1‑251 (MQ=255)
                                                                                                                                                                                                        aaTAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                     >  1:422496/1‑251 (MQ=255)
                                                                                                                                                                                                        aaTAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAGCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATGCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                     <  1:334372/251‑1 (MQ=255)
                                                                                                                                                                                                              cccAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTcgcg                               <  1:468905/250‑1 (MQ=255)
                                                                                                                                                                                                                   aaTACGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCtt                          <  2:402761/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTTCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                      <  2:429200/251‑1 (MQ=255)
                                                                                                                                                                                                                          ttCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAGCCATGGGGAAAAAACCCATCGGGATTAAACAAACCCACCCCGGCGATGGGCCATTTTACGGAAGAAAGAATCCTGGCTTtcgtccccgtcttttttgcc                   >  2:181779/1‑236 (MQ=255)
                                                                                                                                                                                                                           tCA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGcc                   <  2:87061/250‑1 (MQ=255)
                                                                                                                                                                                                                            cA‑CCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTTTTTGcc                   >  1:267845/1‑249 (MQ=255)
                                                                                                                                                                                                                                       agcaAATGCCGGACCGTCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCTTGGCGTTTCCAGCTGGAAGTGAAAGCCAGCGAACAGCGCCAGTTCAATGGCGGTCGCCAGTGTGACCGCAAAGCGAATATGTTGGTTAGCAATGGTGAAAATGCCCATCGTGATTAACCAAACTCGCTCAGGCGAGGGGCCATTTTACGGAAGAACGAATCCTTGCTGTCTTAGCGATCTTGGTTGCCTGTTGCcaccac       <  2:300821/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAACCTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCGGTGACCACCACTGTAg  <  2:46256/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    
ATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCA‑CCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAG  >  NC_000913/3387800‑3388285

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