Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,058,653 T→C 6.8% intergenic (‑243/+13) serA ← / ← rpiA D‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,058,6530TC 6.8% 173.3 / 3.1 59intergenic (‑243/+13)serA/rpiAD‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (26/29);  new base C (2/2);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

CGCCTTCTACCAGCAGAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACTTCTACTGGCAGCGGGAATTTACCC  >  NC_000913/3058356‑3058948
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                       
cGCCTTCTACCAGCAGAAACTTAATCTTGTCTTTCTCCGGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGAGTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTcc                                                                                                                                                                                                                                                                                                      <  1:507214/301‑1 (MQ=255)
      cTACCTGCATAAACTTAATTTTGTCTTTCTCTGGCGATTCCTTTTCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGtt                                                                                                                                                                                                                                                                                                <  1:81656/301‑2 (MQ=255)
       tACCAGCAGAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAACGGTAATTCTTGGATGAATATGCTGGTTTAGGGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACAAGTCACCAAATTTAATGAAGAAAATTTTTTTAACGGGGGAGGTTCCCCCGTCa                                                                                                                                                                                                                                                                                               >  1:182315/1‑301 (MQ=255)
        actagcagAAAATTAATCTTGTCTTTTTCCAGCGATTCCTGTGCCAGTTACCCAATCCTGTCTGTTGAAATGTGGTGTGCGGAGTTGCATCCGCCTGTCAACATATCAAAAAATAATATTGCGGCAAGATGAACGTTTGCGTCGCGATGTTGAGGAAATAGCACGCAAAGGTAATTCTTGGATGAAGATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAAGTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg                                                                                                                                                                                                                                                                                              <  1:34748/298‑1 (MQ=255)
           agcagAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAAGCCTGGCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATc                                                                                                                                                                                                                                                                                           <  2:113310/301‑1 (MQ=255)
              agAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGAAGGTTCCCCAGTCAGATCAtt                                                                                                                                                                                                                                                                                        >  2:444322/1‑301 (MQ=255)
               gAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCAttt                                                                                                                                                                                                                                                                                       >  1:377990/1‑301 (MQ=255)
                ataCTTAATCTTTTCTTTCTCCAGCGATACCTTTGCCATTTACCCCATCCTGTCTTTTGAAATGTTGTTTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAACCTCCCCCGttaaaaaaattc                                                                                                                                                                                                                                                                                      <  2:368483/299‑12 (MQ=255)
                aaaCTTAATCTTGTCTTTCTCCAGCTATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTc                                                                                                                                                                                                                                                                                      <  2:66864/301‑1 (MQ=255)
                   ttttaTCTTGTCTTTCTCCATCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGAGTTGCATCCGCCTTTCAACATATCAAAAAATAAGTTTGCGGGAATAGGAACGTTTGCGTCGGGATGTTGAGGAAATATCACGCAAAGGTAAGTCTTGGATGAATATGCGGGTTTAGAGGAGTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGGAGAGAATTTTTGTAACGGGGGAGGTTCCCCCGTCAGATCGTTTcaca                                                                                                                                                                                                                                                                                   <  1:457944/297‑1 (MQ=255)
                     tAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACaaa                                                                                                                                                                                                                                                                                 >  2:58043/1‑300 (MQ=255)
                      aaTCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATg                                                                                                                                                                                                                                                                                >  2:237542/1‑301 (MQ=255)
                                            aCCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGc                                                                                                                                                                                                                                                          >  2:205400/1‑301 (MQ=255)
                                                tttccgatctcccAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACAcctgacggtgtcaaaac                                                                                                                                                                                                                                                      <  2:241158/291‑16 (MQ=255)
                                                  gCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCa                                                                                                                                                                                                                                                    <  1:57975/301‑1 (MQ=255)
                                                          cccAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACAcctgacggtgtcaaaacagatcggaaa                                                                                                                                                                                                                                            >  1:240891/1‑276 (MQ=255)
                                                                          gAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAg                                                                                                                                                                                                                            >  2:452832/1‑301 (MQ=255)
                                                                           aaaTGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGc                                                                                                                                                                                                                           <  2:378486/301‑1 (MQ=255)
                                                                                tgaaaggcggatgcaaaTCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACa                                                                                                                                                                                                                      >  2:405111/17‑301 (MQ=255)
                                                                                 tgtgtgCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACaa                                                                                                                                                                                                                     <  2:182512/301‑1 (MQ=255)
                                                                                  ttgtgCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAg                                                                                                                                                                                                                    <  2:511378/300‑1 (MQ=255)
                                                                                  gtgtgCGTATTTGCATCCGCCTTTTAACATATCAAAAAATAATATTTCGGCAATATGAACGTTTGCGTCGCGATGTTTAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAg                                                                                                                                                                                                                    <  1:205157/301‑1 (MQ=255)
                                                                                         gATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGTTGTGCCAATCAGCGCAACGTCCGCCCCACGGTTAGAAAACAAGCCAACAg                                                                                                                                                                                                             >  1:311211/1‑301 (MQ=255)
                                                                                            ttGCATCCGTCTTTCAACATATCAAAAAATAATATTGCGGGAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCAAGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTca                                                                                                                                                                                                          <  2:411448/301‑1 (MQ=255)
                                                                                                            aCATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGACAACAGACACCACGCCCAGAAACGc                                                                                                                                                                                          >  1:576843/1‑301 (MQ=255)
                                                                                                                 tCAAAAAATAATATTGTGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTa                                                                                                                                                                                     <  2:81882/301‑1 (MQ=255)
                                                                                                                     aaaaTAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAACCAACAGTCACCACGCCAAGAATCGCATTTATGGc                                                                                                                                                                                 >  2:385566/1‑301 (MQ=255)
                                                                                                                           attttGCGGCAATATGAACGTGTGGGTGGTGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATGTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTc                                                                                                                                                                           <  1:237275/298‑1 (MQ=255)
                                                                                                                                 gggCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAATCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGGCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGc                                                                                                                                                                     <  1:445292/300‑1 (MQ=255)
                                                                                                                                     aaTATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATc                                                                                                                                                                 <  2:19249/301‑1 (MQ=255)
                                                                                                                                             cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCAATCGGGTCAAg                                                                                                                                                         >  1:92860/1‑301 (MQ=255)
                                                                                                                                                 tGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGAAtt                                                                                                                                                     >  2:465232/1‑301 (MQ=255)
                                                                                                                                                    gTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGTGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATAGGGTCAAGGATTTcc                                                                                                                                                  >  2:637814/1‑301 (MQ=255)
                                                                                                                                                          aTGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCg                                                                                                                                            >  2:630612/1‑301 (MQ=255)
                                                                                                                                                               ccgatcgttgcACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGAc                                                                                                                                       <  1:591323/291‑1 (MQ=255)
                                                                                                                                                                     atctCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAg                                                                                                                                 <  2:85614/298‑1 (MQ=255)
                                                                                                                                                                      tatCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGAcagatcgg                                                                                                                                >  2:592084/1‑294 (MQ=255)
                                                                                                                                                                         cACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAgagat                                                                                                                             >  1:85490/1‑297 (MQ=255)
                                                                                                                                                                          aCGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCa                                                                                                                            <  1:629735/301‑1 (MQ=255)
                                                                                                                                                                              aaaTGTAATTCTTGGTGTAATATGCTGGTTGTGAGGGTTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGATTTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGGTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGtt                                                                                                                        <  1:452268/301‑1 (MQ=255)
                                                                                                                                                                                aGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGc                                                                                                                      <  2:92984/301‑1 (MQ=255)
                                                                                                                                                                                        tttGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCg                                                                                                              <  2:311615/300‑1 (MQ=255)
                                                                                                                                                                                              ataaTATGTTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAAGCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTcacc                                                                                                        <  1:26190/299‑1 (MQ=255)
                                                                                                                                                                                                        ggTTTTGAGGTTTTATTTAGAGTAATTGACTATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAGTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTGTTGACACCGTCATGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGAc                                                                                              <  1:391150/301‑1 (MQ=255)
                                                                                                                                                                                                                ggATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTAttac                                                                                      >  2:398333/1‑299 (MQ=255)
                                                                                                                                                                                                                     tatttaGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGgacg                                                                                 <  2:159953/301‑1 (MQ=255)
                                                                                                                                                                                                                        ttaGAGCAATCGACAATTGACTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGAcc                                                                              <  2:577558/301‑1 (MQ=255)
                                                                                                                                                                                                                             gCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCa                                                                         >  2:349051/1‑301 (MQ=255)
                                                                                                                                                                                                                                  tctcAATTGCTTGTTAAAAGCGTGACACATGTCACCAAATTTAATGAAGTGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTc                                                                    <  1:385127/298‑1 (MQ=255)
                                                                                                                                                                                                                                     cAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTACaga                                                                 >  2:385623/1‑299 (MQ=255)
                                                                                                                                                                                                                                                 aaaaGCGTGAACCATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGcc                                                     <  1:404548/301‑1 (MQ=255)
                                                                                                                                                                                                                                                   tttggTGACATGTGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCACTCAGCGCAACGTCAGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGTTCACCACGCCCTGACGGTATTCCGGACCACCGCCCAGACTCACCCGCTGACACCCCAc                                                   >  1:422390/5‑301 (MQ=255)
                                                                                                                                                                                                                                                                            cccttcccgacgctcttccgatctggAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACt                          <  1:87190/277‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             tgaCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGACCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCACCGGGAAAACTTCAACCGGCACCGGGaa         >  1:154329/3‑301 (MQ=255)
                                                                                                                                                                                                                                                                                              aaCGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACCCGCCACTGCACTACCTGCCCTCGGGAACACGTACACTCGCACCGGGAAt        >  1:534214/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               aCGGGGGAACCTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACTTCTACTGGCAGCGGGAAtt       >  2:1807/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                cGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTTACGCGCCACTGCACTACGTGCCATCGGGATAACTTCCACTGGCAGCGGGAAttt      >  1:594066/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  ggggAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACTTCTACTGGCAGCGGGAATTTaa    >  2:14512/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                    ggAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACtagatcggaagagcaaacgtctgaa  >  2:87314/1‑277 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                       
CGCCTTCTACCAGCAGAAACTTAATCTTGTCTTTCTCCAGCGATACCTTTGCCATTTACCCAATCCTGTCTTTTGAAATGTTGTGTGCGGATTTGCATCCGCCTTTCAACATATCAAAAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAGATCATTTCACAATGGTTTTGACACCGTCAGGTGTGCCAATCAGCGCAACGTCCGCGCCACGGTTAGCAAACAAGCCAACAGTCACCACGCCAGGAATCGCATTTATGGCGTTTTCCATCGCTATCGGGTCAAGGATTTCCATGCCGTGGACGTCGAGGATCACGTTGCCATTATCGGTCACCACGCCCTGACGGTATTCCGGACGACCGCCCAGTTTCACCAGCTGACGCGCCACTGCACTACGTGCCATCGGGATAACTTCTACTGGCAGCGGGAATTTACCC  >  NC_000913/3058356‑3058948

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