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% 159.5 / NA 56T50P (ACG→CCG) aaeBaromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (26/28);  total (27/29)
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

ACAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGT  >  NC_000913/3387794‑3388283
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
acaATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGACTCATCAATAATGGTGCGAGACTCATCGCAATGATAATCACCAGTCCAGCAATACAGCCACTAAATGTGCCAATGATGCGCAAAAAGCCCCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGACCGCAcc                                                                                                                                                                                                                                                  >  2:499397/1‑251 (MQ=255)
     gATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGCTCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggcgg                                                                                                                                                                                                                                             >  1:489915/1‑251 (MQ=255)
       tCAGCGCGGTATAACCGGACAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGAGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAAAAGCGCCAGAATAAGGTTCACCTCCCGCAGCAAAGTCCGGACCggcggcc                                                                                                                                                                                                                                           >  2:621911/1‑250 (MQ=255)
                 ataACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGAAACATTcgccg                                                                                                                                                                                                                                 >  1:103228/1‑251 (MQ=255)
                 ataACCGGCCAGCCCCCACGCATACGAGTTCTCAATTCGTACCAGCGAGGATATCCAGGGACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAAGTGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgccg                                                                                                                                                                                                                                 <  2:9865/251‑1 (MQ=255)
                    cccGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATAAAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCATCAACGGCCGGACCGGCGGCAACAATCGCCGCTg                                                                                                                                                                                                                              >  1:299449/2‑251 (MQ=255)
                                  aCGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAg                                                                                                                                                                                                                <  1:76960/251‑1 (MQ=255)
                                             ttttCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGGACCGCCCAGCGTGGCGTTTc                                                                                                                                                                                                     >  1:157937/1‑251 (MQ=255)
                                                         aCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTg                                                                                                                                                                                         <  1:628953/251‑1 (MQ=255)
                                                          ccAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTg                                                                                                                                                                                         >  1:237082/1‑250 (MQ=255)
                                                            aGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGaa                                                                                                                                                                                       >  1:627507/1‑250 (MQ=255)
                                                               gAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                                                                   >  1:26653/1‑251 (MQ=255)
                                                                  gagatCCAGGTACAAAATCCTGCCCATATACGGAACACCAGAATCAGCAATAATGGTGCGCGGATCATCGCAATGAAGAACACCAGGCCGGCAATACAGCCAATAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGGCGCTGTCAGTACCGCCCAGTGTGGCGTTTCCAGCTCGAAGTGATAGCcaac                                                                                                                                                                                <  2:502898/248‑1 (MQ=255)
                                                                   atatCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaaca                                                                                                                                                                               <  1:621911/251‑1 (MQ=255)
                                                                   atatCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGGGAAAGCcaac                                                                                                                                                                                >  2:677598/1‑250 (MQ=255)
                                                                      tCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAACGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcaaaaaaa                                                                                                                                                                             >  1:448962/1‑243 (MQ=255)
                                                                        cAGGTACAAAACCCGGCCCAAATATAGCACACCAGAATCAACAATAATGGTGCGCGGATCAACGCAATGATGAACACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTGCCAGCTGGAAGTGAAAGcgaccaaccag                                                                                                                                                                          <  2:154958/251‑11 (MQ=255)
                                                                         aGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgc                                                                                                                                                                         <  1:369648/251‑1 (MQ=255)
                                                                              cAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGcc                                                                                                                                                                      >  1:577948/1‑249 (MQ=255)
                                                                               aaaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGt                                                                                                                                                                   <  1:575391/251‑1 (MQ=255)
                                                                                           cgccttaccgatcaaGAATCAGAAATAATGGTGCGCGGATCATCGCAATGATGACCCCCAGGCCGACAATACAGCCAATAAATTTGCCGAGGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTc                                                                                                                                                       <  2:109928/237‑1 (MQ=255)
                                                                                              tACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGACAGTACAATGGCggggccc                                                                                                                                                    >  1:454140/1‑246 (MQ=255)
                                                                                               ccAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCa                                                                                                                                                   <  1:617098/250‑1 (MQ=255)
                                                                                                 aGCACAACAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAACGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACACACGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCACTAGAAGTGTAAGCCAACAAACCACGCCCGTACCATGGCggatgccagg                                                                                                                                                 >  2:311875/1‑243 (MQ=255)
                                                                                                    acacCAGAAGCAACAATAATGGTGCGGGGAACATCGCAATGATGATCACCATTCCGGCAATACAGCCAATAAATGTGCCGATAATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGATGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCTAGtttt                                                                                                                                              <  2:65500/251‑1 (MQ=255)
                                                                                                         aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATACATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCCAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAGGCCACCAAACTGCGCCAGTACTATGGCGGTcagggcggtggag                                                                                                                                           >  1:109928/1‑237 (MQ=255)
                                                                                                                          ggcgcGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGt                                                                                                                        <  1:330933/250‑1 (MQ=255)
                                                                                                                             gcgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCGCCAAGCGaaaattgtgggttg                                                                                                                     >  1:570210/1‑239 (MQ=255)
                                                                                                                                         aaTGATGATCACCAGTCCGGCAATACAGCCTATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGaaaa                                                                                                         <  1:311875/251‑1 (MQ=255)
                                                                                                                                               gatCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAccc                                                                                                    <  1:411192/250‑1 (MQ=255)
                                                                                                                                                  aaCCAGTCCGGCAATACAGCCAAGAAATGGGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCTAGCTGGAAGTGAAAGCCAACAAACCGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAATGCGAATGTGTTGGTTAGCAATGGAGAAAATCCCCATCg                                                                                                <  1:48752/250‑1 (MQ=255)
                                                                                                                                                    ccAGTCCGGCATAACAGCCAATAAATGTGCCGAGGACGCGCAAAAAGCCACGATAGCGAATAGGGCCAGAATACGGTTCACCGCCCGCATCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTGAGCAATTGAGAAAATGCCCATCGTg                                                                                              <  2:299449/251‑1 (MQ=255)
                                                                                                                                                         ccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTaa                                                                                         <  1:677598/251‑1 (MQ=255)
                                                                                                                                                                   aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTAGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCCCAAAGAGGATAAGTTGGTTTGCAATGGGGAAAATACCCACCGGGGTTAACCAAACttcc                                                                              >  1:206858/1‑248 (MQ=255)
                                                                                                                                                                            aTGTGCCGATGATGCGCAAAAAGCCACGAACGAGAATTGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGAt                                                                     <  1:499397/251‑1 (MQ=255)
                                                                                                                                                                              gtgCCGAGTATGCGCAAAAAGCCACTATATCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGCCCGGCGGCAACAATCGCCGCTGGCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAATCATCGACAGTACAATGGCGGTCGCCATTTTTACCTCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCGCGCAGGCGATgg                                                                   <  2:21319/251‑1 (MQ=255)
                                                                                                                                                                                gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGGAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATgggg                                                                 >  1:193164/1‑250 (MQ=255)
                                                                                                                                                                                gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGAGCCCGATTTACCGCAAAGCGCATATGTTGGGTAGAGATGGAGAAAATACCCATCATGATTAATCAAACTCACCCAGGCGCTGGGc                                                                 >  1:418990/1‑251 (MQ=255)
                                                                                                                                                                                    atgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGGCAATATTACCGCAAAGCGAATATGTTGGGTTGCAATGGAGAGAATAACCCTCGTGATTAACACAAAccacccacgggcgaggccatt                                                              >  2:82208/1‑231 (MQ=255)
                                                                                                                                                                                     tgatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCCAAGCGAAAATGTTGGTTTGCAATGGAGAAAATACCCATCGGGGTTAACCAAACTCACGCCGGCGAAGGGCCAttt                                                            >  1:11728/1‑251 (MQ=255)
                                                                                                                                                                                          cgcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACgg                                                       <  1:535627/251‑1 (MQ=255)
                                                                                                                                                                                                aaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGa                                                                                             <  1:546370/208‑1 (MQ=255)
                                                                                                                                                                                                     aCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAACGCAAAGCGAATATGTTGGTTAGCAATGGAGGAAAAACCCATCGTGATTAACCAAACTCACGAAGGCGATGGGCCATTTTACGGAAGAACGaaac                                            >  1:625472/1‑249 (MQ=255)
                                                                                                                                                                                                       gagatcgcaTAGCTCCAGAATACGGTTCACCTCCCCCAGAAAAGGCCGGTCCGGAGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACATTGGCGGTCGCCTGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAATCAAACTCACGCAGGCGATGGGCCATTTTACGGAATAACGAATCCt                                          <  2:448962/243‑1 (MQ=255)
                                                                                                                                                                                                              aaTAGCGCCAGAACACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                   <  1:571517/251‑1 (MQ=255)
                                                                                                                                                                                                               aTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAAAGAATCCTGGCTTtcg                                  >  1:408519/1‑251 (MQ=255)
                                                                                                                                                                                                                             cGGTTCCCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                    <  2:627507/251‑1 (MQ=255)
                                                                                                                                                                                                                             cGGTTCACCTCCCGCAGCAAAGGCCAGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTAAGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                    >  1:89537/1‑251 (MQ=255)
                                                                                                                                                                                                                              ggTTCCCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGACGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTTACGCAGGCGATGGGCCATTTTACGGAAGATCGAATCCTGGCTTTCGTCGCGATCTTGTTTg                   <  1:429489/251‑1 (MQ=255)
                                                                                                                                                                                                                               gTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAc                                                                                                                                            >  1:8936/1‑130 (MQ=255)
                                                                                                                                                                                                                               gTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAc                                                                                                                                            <  2:8936/130‑1 (MQ=255)
                                                                                                                                                                                                                                       cccGCAGCAAAGGCCGTACCGGCGGCAACAATCGACGCTTTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGTCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAATAACGAAGCCTGGCTTTCTTCGCGATCTCGTTTGCCAGTgccc          <  2:625472/251‑4 (MQ=255)
                                                                                                                                                                                                                                       cccGCAGCAAAGGCCGGACCGGCGGCAACAAGCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCGTTCGTCGCGATCTTGTTTGCCAGTGa            <  2:408519/249‑1 (MQ=255)
                                                                                                                                                                                                                                        ccGCAGCAAAGGCCGGACCGGCGGCAACAAACGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAAGCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGacca         <  2:193164/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCGCAAAGCGGATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGATTAACCAAAATCACGCAGGCGGTGGGCCATTTTAAGGGAGAAAGAATCCTGGCTTTCGTCGCGGTCTTGTTTGCCAGTGaccaccc      >  1:531209/1‑250 (MQ=255)
                                                                                                                                                                                                                                               aaaGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCGCAAAGCGAAAATGTTGGGTAGCAATGGAGAAAAAACCCATCGGGAATAAACAAAATCACGCAGGCGAAGGGCCATTTTTCGGAAAAACGAATCCTGGCTTTCGGTGCGGTATTTTTTTCCAGTGACCCCCACTgg  >  1:98349/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
ACAATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGT  >  NC_000913/3387794‑3388283

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