Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 587,266 G→C 6.8% A239G (GCC→GGC)  ybcH ← PRK09936 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913587,2660GC 6.8% 170.2 / 2.6 53A239G (GCC→GGC) ybcHPRK09936 family protein
Reads supporting (aligned to +/- strand):  ref base G (26/23);  new base C (2/2);  total (28/25)
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

CCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAATTTCCGCGC  >  NC_000913/586981‑587557
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                   
ccGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACttt                                                                                                                                                                                                                                                                                      <  2:559578/301‑1 (MQ=255)
  tgtcttgATTTTTTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCCCTGCCAgcggcgttgtttatg                                                                                                                                                                                                                                                                                    <  2:637635/295‑10 (MQ=255)
       cGATCTTCTTTCTTTAGAACGCTTTTTAGCAATTAATCTTGATGGAATTCTGTTGAGGGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAGTGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATGGAGAGTGCCTTTACCGCAGGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTCCGCTGTAAAGGTTTTGCCTTTGCCGGCGAAAAAAAGTTCATAAACAACCCCCCTGCCAgcggcgttgttattgagtat                                                                                                                                                                                                                                                                               >  1:636758/1‑287 (MQ=255)
                 tctttTGCATGCTTTTTTGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGcc                                                                                                                                                                                                                                                                     <  1:220553/301‑1 (MQ=255)
                 tctttAGCACGCTTTTTAGCAATTAATCTTGTTGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAAGGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGcc                                                                                                                                                                                                                                                                     <  1:416237/301‑1 (MQ=255)
                             tttttAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCTAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTAGCAGGGGCAGGACTTTGGCAATCGGCACTGTCCTGTAAATCACGt                                                                                                                                                                                                                                                         >  2:300247/1‑301 (MQ=255)
                                                          gatgaGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCCAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACCACGCCGCTGGCAGGGGCCGGAATTTGGCCATCAGACCTGGCCAGGAAAAAACGGTCACACGTTTCATCGGCCCGGTTAttccc                                                                                                                                                                                                                            >  1:337098/1‑297 (MQ=255)
                                                                            aaGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTAc                                                                                                                                                                                                          <  2:424959/301‑1 (MQ=255)
                                                                              gCCAGGTCGTACCCGACTTACCTGGCGGAGATTTCATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGCACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTAccc                                                                                                                                                                                                        >  2:610433/1‑301 (MQ=255)
                                                                                cAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCaa                                                                                                                                                                                                      >  2:300569/1‑300 (MQ=255)
                                                                                       tACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTaa                                                                                                                                                                                               <  2:337503/301‑1 (MQ=255)
                                                                                        aCCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACACAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAAc                                                                                                                                                                                              >  2:113875/1‑301 (MQ=255)
                                                                                                tactgggaggaGATTTAATACTGTGGAGTGCGGTGCGCTCCGGGGCACTAGCGTAGAGGGTAGTTGGGAGTTTCTTTACCTCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGTCTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTg                                                                                                                                                                                      <  1:593490/296‑1 (MQ=255)
                                                                                                tACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTNTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAANNNTCATAAACAACGCCACTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGCAAATAACGTTCACGCTGTTAAGCGGTCCGTTTGTCCACGCCCCCGCCTCCCTGTCCCCACAAATCAACACCCGttt                                                                                                                                                                                      >  1:468197/1‑300 (MQ=255)
                                                                                                        ttcttccgatctcTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACg                                                                                                                                                                              <  2:352499/289‑1 (MQ=255)
                                                                                                                 ttACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCcagc                                                                                                                                                                     <  2:331629/300‑1 (MQ=255)
                                                                                                                    cTCGAGAATGCCGTGCGCGACGGGCGAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGATGGCCTGTGAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCAAGACATAAAAGCCGGATGCTTAAAAgttatacagaagt                                                                                                                                                                  >  2:393401/1‑290 (MQ=255)
                                                                                                                    cTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTTCTTTAACgaaaccggacgag                                                                                                                                                                  >  1:352036/1‑289 (MQ=255)
                                                                                                                              ccGTGCGCGACGGGTAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGc                                                                                                                                                        <  2:522162/301‑1 (MQ=255)
                                                                                                                                     ggACGGGGAAGTAGGGGAGGGAGAAATGGATGGTGTGTTTACCGGAAGAGTAAGGTTTCGCTAATGAGGAGGTAGTTTCTGGGTCCGGTTTCGGTTTAGGTGTAAATGTTTTGCCATTGCCGGAGACAAAAAGTTCATAAACAACGACGCTGGAAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGGTCACGCTGTTCAGCGGTCAGTTGATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTGGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCggg                                                                                                                                                 <  1:609575/300‑1 (MQ=255)
                                                                                                                                     ggACGGGCAAATGTGGCAGAGGGAATTAGGGAGTGTCTTTACCGCAAGGGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCggg                                                                                                                                                 <  1:300165/300‑1 (MQ=255)
                                                                                                                                     cgACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAAAGTGTTCCAGCAGTTGGAGATAACCATAggg                                                                                                                                                 >  2:353797/1‑301 (MQ=255)
                                                                                                                                        cGGGCAATTAGCTCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCTATCGCCaactgctgg                                                                                                                                              <  2:294181/301‑9 (MQ=255)
                                                                                                                                        cGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATCACGTTCACGCTGTTCAGCGGTCAGTTTATCCCCGCCGCTGCCATCCTGTACCCAGACCTTAACGCCGGTTGCTTTAACGTGTTCCAGCACTTGGCGAAAGCCATCGGGCgc                                                                                                                                              >  1:480357/1‑300 (MQ=255)
                                                                                                                                                 ggctgagagagAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTcc                                                                                                                                     <  1:245220/296‑1 (MQ=255)
                                                                                                                                                  gcgcAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCCCTGCCATCCTCTACCCAGACATTAACGCCGGTTGCTTTAACGTGCTCCACCACTTCGCCACAGCCCTCGGGCGCCATTTTTccc                                                                                                                                    >  1:13227/1‑300 (MQ=255)
                                                                                                                                                   cgcAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCgg                                                                                                                                   <  2:480958/301‑1 (MQ=255)
                                                                                                                                                      agagagAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGa                                                                                                                                >  1:282091/1‑301 (MQ=255)
                                                                                                                                                       gagagaAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGaa                                                                                                                               >  2:614318/1‑301 (MQ=255)
                                                                                                                                                                          tttACCGCAAGAGGAACTTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAAcc                                                                                                            <  2:468761/301‑1 (MQ=255)
                                                                                                                                                                              ccGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCAATCGGGCTACAAGTTTCCGACAAAAAAAACACTGATACAAACCAGtt                                                                                                        >  1:626100/1‑301 (MQ=255)
                                                                                                                                                                                   agagGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTTCt                                                                                                   >  1:242068/1‑301 (MQ=255)
                                                                                                                                                                                                                       gTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAACCTACTGATATAAACCGGTTTTGCTGCAACATCGCAAATCAGCCGCTACGCGCTGATTCAc                                                               >  1:384497/1‑301 (MQ=255)
                                                                                                                                                                                                                             tttcggtttcgTTTTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTTAATAACGTTCACGCTGTGCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGtt                                                         <  1:313443/301‑1 (MQ=255)
                                                                                                                                                                                                                             tttcggtttcgCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGAGTTGTTTAACCATGtt                                                         >  2:301293/1‑301 (MQ=255)
                                                                                                                                                                                                                                  gtttcgCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCaa                                                    >  2:358055/1‑301 (MQ=255)
                                                                                                                                                                                                                                      ttcTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGc                                                <  2:275190/299‑1 (MQ=255)
                                                                                                                                                                                                                                      gggtGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGGCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGGGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGc                                                <  1:141274/298‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAACAAAGGcag                                              >  1:274887/1‑299 (MQ=255)
                                                                                                                                                                                                                                           tAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAACCTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTAACg                                           >  1:242073/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         ttgccGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTgc                             <  2:242340/301‑1 (MQ=255)
                                                                                                                                                                                                                                                          tgccGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTgcg                            <  2:282671/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aCAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGTGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCa                    <  2:384993/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAg                   >  2:33009/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    cttccgttcgtaagtcccgccccTGCCAGCGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAgg                  <  1:464568/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaaGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGgt                 >  2:438073/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       aaGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAACCCGGTTTTGCTGAAACATCGCTCATCAGCCGCTGCCCGTTGTTTAACCATGTTACCAAAGGCTGCCGGGCGGCTCCGCCGCCCCAGTTCCggccc               >  1:1863/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                             aTAAACAACGCCGCTGGCAGGGGCGGGACTATGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCCGTTTATCCACGCCGCTGCCATCCTGTAACCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTCGCGATAGCCATCGGGCGACATGTTTCCGGCGACACAACTACTGATATAAACCGGTTTTGCTGACACATCGCTAATCAGCCGCTGCGCGTTGTTTCCCCATGATAGCACAGGCTTACGGGCGGCTCAGCTCCGCCAGTTACGGTAGTCCAtt         >  1:473623/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              tAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAAttt        >  2:155074/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                              aaagtcccgccccTGCCAGCGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAggaaatcggaaa        >  2:465132/8‑291 (MQ=255)
                                                                                                                                                                                                                                                                                   aacGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAATTTCcgcg   <  2:282394/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aacGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAATTTCcgcg   <  2:378095/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    acGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAATTTCCgcgc  >  2:163120/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                   
CCGATAACGATCTTCTTTCTTTAGCACGCTTTTTAGCAATTAATCTTGATGGAATTCTGATGAGAGCGAAAGAGGTAAGCCAGGTCGTACCCGACTTACCTGGAGGAGATTTAATACTCGAGAATGCCGTGCGCGACGGGCAAATAGCGCAGAGAGAAATAGAGAGTGTCTTTACCGCAAGAGGAACGTTTCGCTAACAGCGAGGCAATTTCTGCGTCCGGTTTCGGTTTCGCTGTAAAGGTTTTGCCTTTGCCGGCGACAAAAAGTTCATAAACAACGCCGCTGGCAGGGGCGGGACTTTGGCAATCGGCGCTGGCCTGTAAATAACGTTCACGCTGTTCAGCGGTCAGTTTATCCACGCCGCTGCCATCCTGTACCCAGACATTAACGCCGGTTGCTTTAACGTGTTCCAGCAGTTGGCGATAGCCATCGGGCGACATGTTTCCGGCGAAAAAACTACTGATATAAACCGGTTTTGCTGAAACATCGCTAATCAGCCGCTGCGCGTTGTTTAACCATGTTAGCAAAGGCTGACGGGCGGCTTCGCTGCGCCAGTTCAGGTCGTCAATTTCCGCGC  >  NC_000913/586981‑587557

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑