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.6% 180.9 / NA 69T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (23/45);  total (23/46)
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

ATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG  >  NC_000913/3387797‑3388282
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggc                                                                                                                                                                                                                                              >  2:376513/1‑251 (MQ=255)
          cgGTAAAACCGGCCAGCCCCAAAGAATACGAGTGTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGAGCACCAGTCCGGCAATACAGCCAATAAATGTGCCGAGGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGACGGACCGGCGGcaacaa                                                                                                                                                                                                                                      <  2:397018/249‑1 (MQ=255)
                aaCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCt                                                                                                                                                                                                                              >  1:275224/1‑251 (MQ=255)
                aaCCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACAAGTCCGGCAATACAGCCAATAAATGTGCCGATGCTGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGCCCGGCGGCAACAATCGCCGCt                                                                                                                                                                                                                              >  2:86056/1‑251 (MQ=255)
                     gCCAGCCCCCACGCATCCGAGTTTTCTATTCGTACCAGCGAGGATATCCCGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAAGGAGGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAg                                                                                                                                                                                                                         <  2:433968/251‑1 (MQ=255)
                          cccccACGCATACGCGTTTTCCATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCCCACCAGAATCATCAATAATGGTGCGCGGGTCCTCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAACTGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCATcgccacttgcaataac                                                                                                                                                                                                                     >  2:382656/1‑238 (MQ=255)
                            cccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGcc                                                                                                                                                                                                                  <  1:516065/251‑1 (MQ=255)
                            cccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGAACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACACACGCAGCGGTCAGTACCGcc                                                                                                                                                                                                                  >  2:225914/1‑251 (MQ=255)
                             ccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCGCACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGccc                                                                                                                                                                                                                 <  1:518744/251‑1 (MQ=255)
                             ccACGCATACGAGGTTTCTATTCGTACCAGCGAGGATATCCAGGCACAAAAACCGGCCCAGATACAGCACACCCGAATCATCAATAATGGTGCGCGGATCATCGCAATTATGCTCACCAGTCCGGCATTACAGCACATAAATGTGACGATGATGCGCAAAAAGCCACGCTCGCGACTAGCGCCCGACACCGGTTCACCACCCGCAGAAAAGGCCGGACCGGCGGCAACAATCCCCGCTTTCAGTACCGccc                                                                                                                                                                                                                 >  2:411889/1‑251 (MQ=255)
                                  cATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCCCGATAGCGAAAAGCGCCAGCATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATAGCCGCCGTCCGTACCGCCCAGCgg                                                                                                                                                                                                            >  2:51914/1‑250 (MQ=255)
                                  cATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCACAAATCGCACCTGTCAATACCGCCCcgcgg                                                                                                                                                                                                            >  2:10996/1‑250 (MQ=255)
                                               tATTCGTCCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGCCCGGCGGCAACAATCGCCGCTGTCCGTACCGCCCAGCGTGGCGTTTCCCGCt                                                                                                                                                                                               >  2:196083/1‑251 (MQ=255)
                                               tATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTACAGCt                                                                                                                                                                                               >  1:443459/1‑251 (MQ=255)
                                                  tcggaacatcGAGGACATCCAGGATCAAAAACCGGCCCAGATTCAGCTCACCAGAAACATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGcggg                                                                                                                                                                                             <  2:468491/241‑4 (MQ=255)
                                                    gTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGa                                                                                                                                                                                            >  1:334555/1‑249 (MQ=255)
                                                      aCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATTATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCATCAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGt                                                                                                                                                                                         <  1:472659/250‑1 (MQ=255)
                                                       ccagctgggattgccAGGTACAAAAACCGGCCCAGATACAGCCAACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAGCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAgcga                                                                                                                                                                                       <  2:38643/238‑4 (MQ=255)
                                                         aGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGaaa                                                                                                                                                                                     <  1:30615/251‑1 (MQ=255)
                                                             aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                                  >  2:171817/1‑250 (MQ=255)
                                                                      aGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgc                                                                                                                                                                        >  1:117258/1‑251 (MQ=255)
                                                                                 ccGGCCCATATACATCCCACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAATCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCTGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAt                                                                                                                                                             <  1:411889/251‑1 (MQ=255)
                                                                                    gCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGAGCACCAGTCCGGCAATACAGCCAATAAAGGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGc                                                                                                                                                          <  2:492693/251‑1 (MQ=255)
                                                                                    gCCCAGATACAACACACCAGAATCATCAATAATGGTGCGCGGATCATCACAATGATAATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTCCCATCTGGTAGTGAACGACAACAAAACGCGCAAGTACACTGGc                                                                                                                                                          >  1:233622/1‑251 (MQ=255)
                                                                                            acatcacaCCAGAATCATCAATAATGGTGCGCGGATCATCGAAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCa                                                                                                                                                 <  2:307877/247‑1 (MQ=255)
                                                                                             cAGCACTCCAGAATCAGCAATAATGGAGCGCGGATCATAGAAATGATGATCACCAGCCCGGCAATACAGCAAATAAATGTGCAGATGATTCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGCTCACATCACGCAGCAAAGGCCGGACCGGAGGCAACAACCGACACTGTAACTACAGACCCGCGTGGCGTTTCAAGCTGGCAGGGAAAACAAACAAACAGCGACAGAACaaagggagacgcccg                                                                                                                                                 >  2:516352/1‑238 (MQ=255)
                                                                                                cacaCCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCTCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGttt                                                                                                                                             >  1:31695/1‑251 (MQ=255)
                                                                                                  acccAGAATCATCAATAATGGGGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTa                                                                                                                                           <  2:388363/249‑1 (MQ=255)
                                                                                                    ccAGAATCAACAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGACCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCAGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAcc                                                                                                                                         <  1:29646/251‑1 (MQ=255)
                                                                                                     cAGAATCATCAATAAGGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCg                                                                                                                                        <  2:479829/251‑1 (MQ=255)
                                                                                                     cAGAATCAACAATAATGGTGCGCGGATAACCGCAATGTTGATCACAAGCCCGGCATTACCGCCAATAAACGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCTCCTGCCGCAGCAAAGGCCGGCCCGGCGTCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTaacg                                                                                                                                        <  2:471507/251‑4 (MQ=255)
                                                                                                             ccaATAATGGGGCGCTGAGCATCTCAATGATGACCCACAGTCCGGCAATACAGCCAATAAATGCGCCGATGATGCGCAAAAAGCCACGATAGCGATTAGGGCCAGAATACGGTTCACCTCCCGCAGCAAATGCCGGTCAGGCGGCAACAATCGCCGCTGTCAGTACCGCCGAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCTGTCGCCAGTTTTACCGCAAAGCGa                                                                                                                                <  2:537211/250‑1 (MQ=255)
                                                                                                              caATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAACGa                                                                                                                                >  1:39592/1‑250 (MQ=255)
                                                                                                                ataatGGTGCGCGAATCATCGCAAAGATGATAACCAGTCCGGCAATACATCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAata                                                                                                                             <  2:275224/251‑1 (MQ=255)
                                                                                                                    tGGTGCGCGGATCCTCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCACCAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAAAAAACAGCGCCAGTACAATGGCGGTCGACCAGTTTACCCCAAAACGAATATGtt                                                                                                                         >  2:271592/1‑251 (MQ=255)
                                                                                                                    tGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCGATTACAGCTGGAAGAGAAAGCCAACAAACAACGCCAGTACAATGGCGGGCGCCAGTTTTTCCACAAAACGAataatgt                                                                                                                         >  2:333890/1‑247 (MQ=255)
                                                                                                                          gcgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAg                                                                                                                   <  1:376513/251‑1 (MQ=255)
                                                                                                                           cgGATCATCGCAATGATGATCACCAGTCCGGCAAGACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGc                                                                                                                  <  1:419868/251‑1 (MQ=255)
                                                                                                                            gtatcaccgaaaggatgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAACAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGGTGCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGTTCGCCAGGGTTACCTCAAAGGGAATATGTTGGGTAGCa                                                                                                                 <  2:230210/238‑1 (MQ=255)
                                                                                                                             gATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGc                                                                                                                                       <  1:347826/228‑1 (MQ=255)
                                                                                                                             gATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCaa                                                                                                                <  1:237227/251‑1 (MQ=255)
                                                                                                                                catcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATgg                                                                                                             <  1:59778/251‑1 (MQ=255)
                                                                                                                                catcGCAATGATGATCACCAGGCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATg                                                                                                              <  1:451160/250‑1 (MQ=255)
                                                                                                                                     aaaTGAAGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGaaa                                                                                                        <  2:419956/250‑1 (MQ=255)
                                                                                                                                       atgatgATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAAt                                                                                                      <  2:233017/251‑1 (MQ=255)
                                                                                                                                          atgatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAcc                                                                                                   <  1:362992/251‑1 (MQ=255)
                                                                                                                                               cACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCg                                                                                              <  1:24135/251‑1 (MQ=255)
                                                                                                                                                 ccAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTg                                                                                            <  1:182053/251‑1 (MQ=255)
                                                                                                                                                     gacGGCAAGACAGCCACGAAATGCGCCGATGAGTCGCAAACTGCCACGATAGCGAATAGCGCCAGAATCCGGTTCACCCCCCGAAGGAAAGGCCGGACCGGCTGCAACAAGCGACGCTGGTAGTCGCGCCCAGCGTGGCGTTGCCAGCTGGAAGTGAAAGACAACAACCAGCGTCAGGACAATGGCGGTCGACCGTGTTAACGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAAGACCCATCGTGATTa                                                                                        <  1:123514/249‑1 (MQ=255)
                                                                                                                                                               cAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCTGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCa                                                                              <  1:235611/251‑1 (MQ=255)
                                                                                                                                                                aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCTCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCAc                                                                             <  2:330222/251‑1 (MQ=255)
                                                                                                                                                                 gCCAATAAAAGTGACGAGGATTCGCAAAAAGCCACGATAGCGAATGGCGCCATATTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCACCGCAGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACATCGCCAGTACAATGGCGGTCGCCAGTGGTACCGCAAAGCGAATATCGTGGGTAGCAATGGAGAAAATACCCATCGTGATTAACCACACTCACg                                                                            <  2:220419/251‑1 (MQ=255)
                                                                                                                                                                   cAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTAAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCa                                                                          <  1:39041/251‑1 (MQ=255)
                                                                                                                                                                           gtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAATTTTACCGCAAAGAGGATATGTTGGGTAACCAAGGAGGAAATAACCATCGTGATTAACCAAACTCACGCAGGCGaggg                                                                  >  1:436038/1‑248 (MQ=255)
                                                                                                                                                                                        gcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATGCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAATTTACGga                                                     <  1:519707/251‑1 (MQ=255)
                                                                                                                                                                                                 cACGACAGCGAATAGCGACAGAATACGGTGCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCTTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAAGATGTGGGTTAGCAAGGGAGAAAATACCCATCGTGATTAACCTAACTCACGCAGGCTATTGGCCATTTTACGGAATAACGAAt                                            <  2:214372/251‑1 (MQ=255)
                                                                                                                                                                                                    gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCt                                         <  2:311414/251‑1 (MQ=255)
                                                                                                                                                                                                     aTAGAGAATAGCGCCATAATACGGTTCACCTCCCGAAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTGAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAGTTTACGGAAGAACGAATCCTg                                        <  2:423137/251‑1 (MQ=255)
                                                                                                                                                                                                            agaGCGCCATAATACGTGTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACTCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcg                                 <  2:281316/249‑1 (MQ=255)
                                                                                                                                                                                                              cgcgcCAGAATACGGTTCACCCCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCGGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgt                                <  2:460784/249‑1 (MQ=255)
                                                                                                                                                                                                              aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                               <  2:180941/251‑1 (MQ=255)
                                                                                                                                                                                                               gcgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGCCCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtcgtcg                              <  2:185005/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGAt                          <  2:334555/251‑1 (MQ=255)
                                                                                                                                                                                                                      aaTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCCCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGCAGTGAAAGCCAACAAACAGCGCCAGTACCATGGCGGAAGCAACTTTTACCTCAAAGCGAATATATTGGTACGCAAAGGGGGAAATACCCATCGTGATTAACCAAACACACGCAGGCGATGGGCCATATTACGGAAGAACGAAACCGGGCCATCGTCCCGATcct                       >  1:238656/1‑249 (MQ=255)
                                                                                                                                                                                                                       aTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAAAATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGCTGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCttt                      >  1:405158/1‑250 (MQ=255)
                                                                                                                                                                                                                           ggTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTg                  <  1:77435/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGACAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAACGCCAAACGCATATGTTTGGTAGCCATGGGGAAAATACCCAATGTGATTAACCAAACTCACCCAGGCGATGGGCCATTTTACGGGAGAACGGATACCGGCTTTCGTCGGGCTCTTTTTTGc                 >  1:142524/1‑251 (MQ=255)
                                                                                                                                                                                                                              tCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGcc                <  1:368387/250‑1 (MQ=255)
                                                                                                                                                                                                                                           caAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTGCGTCGCGATCTTGTTTGCCAGTGACCACCACTg  <  2:405158/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTG  >  NC_000913/3387797‑3388282

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