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

AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTT‑GTTTGCCAGTGACCACCA  >  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)
                                                                                          gATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑Cg                                                                                                                                                         <  1:128794/251‑1 (MQ=255)
                                                                                                      cAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCg                                                                                                                                           <  2:442910/251‑1 (MQ=255)
                                                                                                         aaTCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGACTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGG‑CGCC‑ACGTTTTAccca                                                                                                                                          >  2:311108/1‑248 (MQ=255)
                                                                                                           acatcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAg                                                                                                                                      <  2:64604/250‑1 (MQ=255)
                                                                                                            catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAg                                                                                                                                      <  1:132083/250‑1 (MQ=255)
                                                                                                                   aatGGTGCGCGGATCATGTAAATGATTATCACCAGTCCTGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑TTTTTACCGCAAAGCG‑Aatat                                                                                                                              <  1:260266/250‑1 (MQ=255)
                                                                                                                    atGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑ACCC‑C‑GTTTTACCGCAACGCG‑AATATgg                                                                                                                            >  1:378618/1‑250 (MQ=255)
                                                                                                                      ggTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGACAGTACAATGGCGGG‑CGCC‑C‑GGTTTACCGCCAAGCG‑AATATG‑TTg                                                                                                                         >  1:429934/1‑251 (MQ=255)
                                                                                                                           tcgaatcatcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGA‑CTCC‑C‑GTTTTAACGCCAAGAG‑CAAATG‑TTGGttt                                                                                                                     >  2:156630/5‑249 (MQ=255)
                                                                                                                                  atcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGga                                                                                                             <  1:369977/251‑1 (MQ=255)
                                                                                                                                             gatCCCCAGTCCGGCAATACAGACAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAGGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TAccc                                                                                                 <  2:90964/250‑1 (MQ=255)
                                                                                                                                             gatCACCAGTCCGGCAATCCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCA‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TAccc                                                                                                 <  2:150706/250‑1 (MQ=255)
                                                                                                                                                      cccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑Ta                                                                                      <  2:125488/250‑1 (MQ=255)
                                                                                                                                                          gCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATTGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCa                                                                                  <  1:285543/251‑1 (MQ=255)
                                                                                                                                                           cAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                                                 <  1:184823/115‑1 (MQ=255)
                                                                                                                                                           cAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                                                 >  2:184823/1‑115 (MQ=255)
                                                                                                                                                                 aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGGC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGGTAGCAATGGGGG‑AAA‑TACCCATCGTGAT‑TAACCAAAATCAc                                                                           >  1:219943/1‑251 (MQ=255)
                                                                                                                                                                 aGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCAc                                                                           <  1:357176/251‑1 (MQ=255)
                                                                                                                                                                   ccAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCACAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGc                                                                         >  1:244786/1‑251 (MQ=255)
                                                                                                                                                                   ccAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGGTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGc                                                                         <  2:109388/251‑1 (MQ=255)
                                                                                                                                                                    cAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAGACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCa                                                                        <  2:234664/251‑1 (MQ=255)
                                                                                                                                                                      aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAgg                                                                      <  1:210986/251‑1 (MQ=255)
                                                                                                                                                                         aaTGTTCCGATGACGCGAAAAAAGCCCCGATAGCGCATAGCGTCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCACCAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCAAGTACAATGGCGGTTCCCC‑G‑GTTTTAC‑ACCAAGCGCAAT‑TG‑GTGGTTTGCAATGGAGA‑AAA‑TATCCATCGCGAT‑TAAACaaaatccccccggggg                                                                   >  2:71123/1‑238 (MQ=255)
                                                                                                                                                                          atttgCCGATGATTCGCAAAAAGCCACGAAAGCGAATAGCGCCAGTATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCTCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCAGCGTGAT‑TAACCAAACTCACGCAGGCGAt                                                                  <  2:352650/248‑1 (MQ=255)
                                                                                                                                                                               ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGAGAAAGCCAACAAACAGCGCCAGTACAATGGCGGC‑CGCA‑A‑GTTTTAACGCAAAACG‑AATATG‑ATGGTTTGCAATGGAGA‑AAA‑AACCCATCGGGAT‑TTACCAAAATCACGAAGGGGATGAGcc                                                             >  1:201755/1‑251 (MQ=255)
                                                                                                                                                                                    gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCAGAGGCAAGAATCGCCGCTGTCAGTACCGCCCAGGGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGAGCCAGAACAAGGGCGGT‑CGCC‑A‑GTTTTGCAGGAAAACG‑AATTTG‑GTGGGTTGGAATGGGGA‑AAAATACC‑ATCGTGGT‑TGACCAAACTCACGGAggagggggggcatttt                                                        >  1:496834/1‑237 (MQ=255)
                                                                                                                                                                                      tgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGGGAAAGCCAACAAACAGCGCCAGTACAATGGCGGG‑AGCC‑A‑GATTTACACCAAAGAG‑AATATGGTAGGTT‑GCCATGGGGA‑AAA‑TACCCCACGAGAT‑TAAACAAACTTACGCAGGCGAGTGGCCAGTTTa                                                       >  2:388910/1‑250 (MQ=255)
                                                                                                                                                                                           aaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTGTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaag                                                 <  1:343671/251‑1 (MQ=255)
                                                                                                                                                                                            aaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaaga                                                <  2:381351/251‑1 (MQ=255)
                                                                                                                                                                                                     gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCAGCGTGAT‑TAACCAAACTCACGCAGGCGGTGGGCCATTTTACGGAAGAACGAATCCt                                       <  2:446161/251‑1 (MQ=255)
                                                                                                                                                                                                         gCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TAACCATCGTGAT‑TAAACAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCt                                   >  1:16260/1‑251 (MQ=255)
                                                                                                                                                                                                           gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGAGGGGCCATTTTACGGAAGAACGAATCCTGGCttt                                 <  1:109585/251‑1 (MQ=255)
                                                                                                                                                                                                            aaTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGTGCCATTTTACGGAAGAACGAATCCTGGCTTtc                                <  2:252807/251‑1 (MQ=255)
                                                                                                                                                                                                             aTAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTTCCGCCAAACG‑AATATT‑GTGGGTAGCAATGGGGA‑AAA‑TACCCCTCGTGAT‑TAAACAAAATCACGCAGGCGGTGGGCCATTTTACAGAAGAACGAATCCTGGCTTtcg                               >  2:414824/1‑251 (MQ=255)
                                                                                                                                                                                                                        aTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAAGGGAGA‑AAA‑TACCCATCGTGAG‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTT‑g                   <  1:496411/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTT‑GTTTGCCAg           <  1:66124/249‑1 (MQ=255)
                                                                                                                                                                                                                                  cctcaagcagcaACGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACAGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGAGGG‑CGCC‑A‑GTATTACCCCAAAGGG‑AATTTC‑TTGGGTAGGCATGGGGAGAAA‑AACCCAT‑GTGGT‑TAAACAAAATCACACAGGCGATGGGACACTTTTCGGAAGAAAGCAACCTGGCTTTCGTCCCGATTTT‑GTTTTCCatgg         >  2:329741/7‑248 (MQ=255)
                                                                                                                                                                                                                                  ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGACCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCCCA‑GTTTT‑CCGCAAAGAG‑AATATG‑TTTGTTAGAAATGGAAA‑AAA‑AACCCCACCTGGTGTAAC‑AAACCCACACCGGCGATGAGCCATTTTTAAGAAGAACGGAACCTGGCTATGGTCGCGATCTT‑GTTTGTCAGTg         >  2:286592/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACACACAGCGCCAGTACAATGGCGGT‑CGCC‑C‑GGTTTACCGCAAAGCG‑CATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGCATCCTGGCTTTCGTCCCGGTCTT‑GTTTTCCAGGGac       >  1:57154/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCCGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGA‑CACC‑C‑GTTTTACCGCCAAGCG‑AATATG‑TTTGTTAGCAATGGAGG‑AAA‑TACCCATCGTGAT‑TAACCAAACTCCCGCCAGCGATGGGCCATTTTACTGAAAAAAGAATCATGGCCTTCGGTGCGATCTTTGTTT‑CCAGTGaccaccc  >  2:124440/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
AATGATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGT‑CGCC‑A‑GTTTTACCGCAAAGCG‑AATATG‑TTGGTTAGCAATGGAGA‑AAA‑TACCCATCGTGAT‑TAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTT‑GTTTGCCAGTGACCACCA  >  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: ‑