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

AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387796‑3388279
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
aaTGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATAAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGCATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGAACcgg                                                                                                                                                                                                                                                    >  1:478495/1‑251 (MQ=255)
                               cACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCAGCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCa                                                                                                                                                                                                                     <  1:184368/251‑1 (MQ=255)
                                  gCATAAGAGTTTTCTATTCGTACAAGCGAGGATATCCAGGGATAAAACCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGAGCGCGGTTCATCGCCATGATCATCACCAGTCCGGCAATACAGACAATAAATGAGCGGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCCAAAGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGACCAGCg                                                                                                                                                                                                                  >  2:221612/1‑251 (MQ=255)
                                    agaagagtgatctattaaaaacAGCGAGGATATCCAGGAACAAAAACCGTCCCGGTTACAGCACACCAGAATCAGCCATAATGGGACGCGGAACATCGCAGTGATGATCACCAGGACGGCAATACCGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCCAAGGCCGGACCACCGGCAACCATCGCCGCTGTCCGTACCGCCCCGCGTg                                                                                                                                                                                                                >  2:65/22‑251 (MQ=255)
                                     tACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTgg                                                                                                                                                                                                               >  1:114985/1‑251 (MQ=255)
                                       cGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGCTACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGAGCAAAAAGCCAAGATAGCGAATAGCGCCAGAATACGGTTCACCTACCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCCCTGTCAGTACCGCCCAGCGTGGCg                                                                                                                                                                                                             >  2:269921/1‑251 (MQ=255)
                                            tttCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTc                                                                                                                                                                                                         <  1:436174/250‑1 (MQ=255)
                                             ttCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGGGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACCGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAAGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTccg                                                                                                                                                                                                       >  1:305955/1‑250 (MQ=255)
                                                 aTTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAAGCAGCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAg                                                                                                                                                                                                      <  1:206956/248‑1 (MQ=255)
                                                              aGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACTGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCTAGCTGGAAGTGAAAGCca                                                                                                                                                                                      <  1:132332/251‑1 (MQ=255)
                                                                   atCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCaacaaa                                                                                                                                                                                 <  1:49800/251‑1 (MQ=255)
                                                                   atCCAGGAACAAAACCCGGCCCAGATACAGCACACCAGAAGCATAAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCcaacaa                                                                                                                                                                                  <  2:89950/250‑1 (MQ=255)
                                                                     ccAGGTACAAAAACCGGCCCAGATACAGCCCACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGACGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCcaacaacc                                                                                                                                                                                >  2:501160/1‑248 (MQ=255)
                                                                             caaaaacagCCCAGATACAGCACACCAGACTCATCAATAATGGTGCGCGGATCATCGCAATGAAGATCAACAGTCAGGCAATACAGCCAATAAATGCCCCGATGCTGCGCAAAAAGCCACGATAGCGAATAGCCCCAGACTCCGCTTCACCCCCCGCAGCAAAGACCGGACCGGCGGCAACAAACGCCACTGTCCGTACCGCCCACCCTGGCGTTTCCAGCTGGATGCGAAAGCCaacaaacaacgacaat                                                                                                                                                                       >  2:43117/9‑243 (MQ=255)
                                                                                aaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGAGCATCGCAATGATGATCTCCAGTCCGGCAATACAGCCAATAAATGGGCCTATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCGGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACa                                                                                                                                                                    <  1:76282/251‑1 (MQ=255)
                                                                                 aCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACaa                                                                                                                                                                   <  1:361051/251‑1 (MQ=255)
                                                                                      cccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCg                                                                                                                                                              >  2:171751/1‑251 (MQ=255)
                                                                                        cAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCggg                                                                                                                                                            >  1:170632/1‑250 (MQ=255)
                                                                                          gATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCg                                                                                                                                                         <  1:128794/251‑1 (MQ=255)
                                                                                                      cAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCg                                                                                                                                           <  2:442910/251‑1 (MQ=255)
                                                                                                         aaTCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCACGTTTTAccca                                                                                                                                          >  2:311108/1‑248 (MQ=255)
                                                                                                           acatcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAg                                                                                                                                      <  2:64604/250‑1 (MQ=255)
                                                                                                            catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAg                                                                                                                                      <  1:132083/250‑1 (MQ=255)
                                                                                                                   aatGGTGCGCGGATCATGTAAATGATTATCACCAGTCCTGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCATTTTTACCGCAAAGCGAatat                                                                                                                              <  1:260266/250‑1 (MQ=255)
                                                                                                                    atGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTACCCCGTTTTACCGCAACGCGAATATgg                                                                                                                            >  1:378618/1‑250 (MQ=255)
                                                                                                                      ggTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGACAGTACAATGGCGGGCGCCCGGTTTACCGCCAAGCGAATATGTTg                                                                                                                         >  1:429934/1‑251 (MQ=255)
                                                                                                                           tcgaatcatcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACTCCCGTTTTAACGCCAAGAGCAAATGTTGGttt                                                                                                                     >  2:156630/5‑249 (MQ=255)
                                                                                                                                  atcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGga                                                                                                             <  1:369977/251‑1 (MQ=255)
                                                                                                                                             gatCCCCAGTCCGGCAATACAGACAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAGGCGAATATGTTGGTTAGCAATGGAGAAAATAccc                                                                                                 <  2:90964/250‑1 (MQ=255)
                                                                                                                                             gatCACCAGTCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCAAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAccc                                                                                                 <  2:150706/250‑1 (MQ=255)
                                                                                                                                                      cccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTa                                                                                      <  2:125488/250‑1 (MQ=255)
                                                                                                                                                          gCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATTGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCa                                                                                  <  1:285543/251‑1 (MQ=255)
                                                                                                                                                           cAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                                                 <  1:184823/115‑1 (MQ=255)
                                                                                                                                                           cAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                                                 >  2:184823/1‑115 (MQ=255)
                                                                                                                                                                 aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGGCAGTTTTACCGCAAAGCGAATATGTTGGGTAGCAATGGGGGAAATACCCATCGTGATTAACCAAAATCAc                                                                           >  1:219943/1‑251 (MQ=255)
                                                                                                                                                                 aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCAc                                                                           <  1:357176/251‑1 (MQ=255)
                                                                                                                                                                   ccAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCACAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGc                                                                         >  1:244786/1‑251 (MQ=255)
                                                                                                                                                                   ccAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGc                                                                         <  2:109388/251‑1 (MQ=255)
                                                                                                                                                                    cAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAGACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCa                                                                        <  2:234664/251‑1 (MQ=255)
                                                                                                                                                                      aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAgg                                                                      <  1:210986/251‑1 (MQ=255)
                                                                                                                                                                         aaTGTTCCGATGACGCGAAAAAAGCCCCGATAGCGCATAGCGTCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCACCAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCAAGTACAATGGCGGTTCCCCGGTTTTACACCAAGCGCAATTGGTGGTTTGCAATGGAGAAAATATCCATCGCGATTAAACaaaatccccccggggg                                                                   >  2:71123/1‑238 (MQ=255)
                                                                                                                                                                          atttgCCGATGATTCGCAAAAAGCCACGAAAGCGAATAGCGCCAGTATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCTCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGAt                                                                  <  2:352650/248‑1 (MQ=255)
                                                                                                                                                                               ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGAGAAAGCCAACAAACAGCGCCAGTACAATGGCGGCCGCAAGTTTTAACGCAAAACGAATATGATGGTTTGCAATGGAGAAAAAACCCATCGGGATTTACCAAAATCACGAAGGGGATGAGcc                                                             >  1:201755/1‑251 (MQ=255)
                                                                                                                                                                                    gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCAGAGGCAAGAATCGCCGCTGTCAGTACCGCCCAGGGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGAGCCAGAACAAGGGCGGTCGCCAGTTTTGCAGGAAAACGAATTTGGTGGGTTGGAATGGGGAAAAATACCATCGTGGTTGACCAAACTCACGGAggagggggggcatttt                                                        >  1:496834/1‑237 (MQ=255)
                                                                                                                                                                                      tgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGAGCCAGATTTACACCAAAGAGAATATGGTAGGTTGCCATGGGGAAAATACCCCACGAGATTAAACAAACTTACGCAGGCGAGTGGCCAGTTTa                                                       >  2:388910/1‑250 (MQ=255)
                                                                                                                                                                                           aaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTGTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaag                                                 <  1:343671/251‑1 (MQ=255)
                                                                                                                                                                                            aaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaaga                                                <  2:381351/251‑1 (MQ=255)
                                                                                                                                                                                                     gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGGTGGGCCATTTTACGGAAGAACGAATCCt                                       <  2:446161/251‑1 (MQ=255)
                                                                                                                                                                                                         gCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATAACCATCGTGATTAAACAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCt                                   >  1:16260/1‑251 (MQ=255)
                                                                                                                                                                                                           gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGAGGGGCCATTTTACGGAAGAACGAATCCTGGCttt                                 <  1:109585/251‑1 (MQ=255)
                                                                                                                                                                                                            aaTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGTGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                <  2:252807/251‑1 (MQ=255)
                                                                                                                                                                                                             aTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTTCCGCCAAACGAATATTGTGGGTAGCAATGGGGAAAATACCCCTCGTGATTAAACAAAATCACGCAGGCGGTGGGCCATTTTACAGAAGAACGAATCCTGGCTTtcg                               >  2:414824/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAAGGGAGAAAATACCCATCGTGAGTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTg                   <  1:496411/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAg           <  1:66124/249‑1 (MQ=255)
                                                                                                                                                                                                                                  cctcaagcagcaACGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACAGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGAGGGCGCCAGTATTACCCCAAAGGGAATTTCTTGGGTAGGCATGGGGAGAAAAACCCATGTGGTTAAACAAAATCACACAGGCGATGGGACACTTTTCGGAAGAAAGCAACCTGGCTTTCGTCCCGATTTTGTTTTCCatgg         >  2:329741/7‑248 (MQ=255)
                                                                                                                                                                                                                                  ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGACCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCAGTTTTCCGCAAAGAGAATATGTTTGTTAGAAATGGAAAAAAAACCCCACCTGGTGTAACAAACCCACACCGGCGATGAGCCATTTTTAAGAAGAACGGAACCTGGCTATGGTCGCGATCTTGTTTGTCAGTg         >  2:286592/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACACACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCGCAAAGCGCATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGCATCCTGGCTTTCGTCCCGGTCTTGTTTTCCAGGGac       >  1:57154/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCCGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACACCCGTTTTACCGCCAAGCGAATATGTTTGTTAGCAATGGAGGAAATACCCATCGTGATTAACCAAACTCCCGCCAGCGATGGGCCATTTTACTGAAAAAAGAATCATGGCCTTCGGTGCGATCTTTGTTTCCAGTGaccaccc  >  2:124440/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387796‑3388279

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