Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,6640CG36.7% 30.1 / 18.0 51pseudogene (940/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base G (6/12);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862428‑3862902
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCgg                                                                                                                                                                                                                                                    <  1:1524975/239‑1 (MQ=255)
  tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGCGGGa                                                                                                                                                                                                                                                  >  1:48787/1‑239 (MQ=255)
  tGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTAGGTAATTCCCACCAGCATGGCGGGTTAAGGTGCgggg                                                                                                                                                                                                                                                  >  1:926614/1‑238 (MQ=255)
  tGGAGATTTGTGGCAGCTTCATATCCGCGTGGTTGCTGAACATCGGGTTCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCGCCGCCCAGGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGa                                                                                                                                                                                                                                                  <  2:1259094/239‑1 (MQ=255)
       aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGCGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGCTCGGTAATTCCCACCAGCATGGCGGCTAAAGTTGCggggaacc                                                                                                                                                                                                                                              >  1:719433/1‑233 (MQ=255)
                            gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAACAGCAAGCCCGCCACCTTTACCCg                                                                                                                                                                                                                      >  1:1805086/1‑239 (MQ=255)
                            gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGTTGCGGGAAAAAGCAAGACCCCCACCTTTAtccc                                                                                                                                                                                                                      >  2:2506911/1‑238 (MQ=255)
                             cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTTGGGGGAACAGGCAAACCCGCCCCTTTTTccga                                                                                                                                                                                                                      >  1:2765707/1‑237 (MQ=255)
                                 gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGTTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGCTAAAGTTTCGGGAAACACGAAGCCCCCCACCTTTTAccgcatt                                                                                                                                                                                                                   >  1:89577/1‑233 (MQ=255)
                                            tCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTTGGGGGATCAACAAACCCCCCCCCTTTTCCCCGTTTTTCCGGGGGGGc                                                                                                                                                                                                       >  1:692300/1‑239 (MQ=255)
                                             cggtgtCCAGTTTTGCAGTAAAACAGAGCCAGACAGACCTCCGCCCATCGTCCCCCCCACCCAAAAGAAATACTTACGGGTCGAAATTGGGGCCGCCAGCAGGGCGTGTACCCAAAAAAGAAACGGGGAAAAGAACAGAAAGGAAAATCCCGGCGTTTCGTTAATCCCCACCAGCAGGGGGTTTAGGTTTGGGGAAATCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                      <  2:1308690/233‑1 (MQ=255)
                                                     aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGAACAGCAAAACCCCCCCTTTTTCCCGCTTTTTCCGGGCGGGAGGGaacaaa                                                                                                                                                                                             >  1:2434753/1‑235 (MQ=255)
                                                      gTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                      <  1:1168283/183‑1 (MQ=255)
                                                        tttGCAGTAAAACCCTGTTCAACAAGCCCCCCGCCATGGTGCCCCACCCCACAAAAGAGATAAATTCGGGCGACCATGGAGGCCGCAAGCAGGGGGTCTACCGAAAACAGAAAGGGAAAAATGAACAGGAAGGTAATTCCCGGGGTTGCGGTAACCCCCCCCACCAGGGGGGTAAGGGTGGCGGAATTCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                             <  2:1600021/236‑1 (MQ=255)
                                                        ggtGCGGTAAAACCTACACAATCAGCCCACGGCCCGGGCACCCCACCACACAGAAGGCTTCCACTCGGGCCGCCATGGGGCCCGCAAGCAGGGCGTGTCCCACAACCGGCAAGGGTAAAATGAAAAGGAAGGTAAATCCCAGGGTTTGGGTAATCCCCCCCAGCCTGGCGTTTAGGGTGGGGGAAATCAGCAAACCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa                                                                                                                                                                                             <  2:692300/233‑1 (MQ=255)
                                                          tGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGAATCAACAAACCCCCCACTTTTTCCCGATTTTTCCGTGCGGCAGTggaaaaaaac                                                                                                                                                                                          >  1:340525/1‑228 (MQ=255)
                                                                aaaaCCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGGTTCGGGGaacaaac                                                                                                                                                                                                                                            >  2:1168283/1‑179 (MQ=255)
                                                                              tCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGGATCAGCAAGCCCCCCCCCTTTTTCCGCGTTTTCGGTGCGGCCGGGAAGGAAAATCGCCAAAGAAATAccccggg                                                                                                                                                                      >  2:135909/1‑234 (MQ=255)
                                                                               gcagccTCCGCCCATGTCGCCCCCCACGCACAAGAGATAGACCACGGTCAGAATTGGGCCCCACGACCCGGGGGGCCCCAAAAACGCGAACGGAAAAAAGAACAGGAAGATAATCCCCAGCGGTACGGAAATCCCCACCAGCAGGCCGGTTGAGGTTGCGGAAATCGGCATGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACgg                                                                                                                                                                      <  2:2315545/232‑1 (MQ=255)
                                                                                                       ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATACGCAAACCCCCCACTTTTACACCACTTTTCGGTGCCGCAGTGAAGTACAACTCTCAAAAAAAAGCCACGGCAGCCAAAACTTTTGTGTTatcatgt                                                                                                                                           >  2:1829381/1‑233 (MQ=255)
                                                                                                                                tCGACGTTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGGTTGGGGGAAAAGCAAACCCCCCCCCTTTTACCCCATTTTCCGGGGGGCCGGGAAGGAAAACCGCTAAGAAAAGCCCCCGGCACCCAAAGATTTTTAGGTTACGGGCCAGGGCAAAACCCCATCCCGGGaa                                                                                                                   >  2:1725538/1‑239 (MQ=255)
                                                                                                                                   aCATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCCGCAAGCCCGCCCCTTTTACCCCGTTTTTCGGTGCGGGAGTGGAAGACATCGCTAAAGAAAAGCCCCCGGCCCCCAAGAATTTTTAGTTTACCTGCACGGGAAAAACCCCCTCCCGGGaacac                                                                                                                >  2:463228/1‑238 (MQ=255)
                                                                                                                                        gcagcccgccAGCGGGGCGGGTACCGCAAAAACAACGTGTAAGATAACAAAGATGATACCTCCCGGCGGTTGGTCAATCCCCACCAGCAGGGCGGTTAAGGTTGCGGGAAACAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTGAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaac                                                                                                                  <  2:2765707/227‑1 (MQ=255)
                                                                                                                                                       ggCTTGAACCGAAAACAAAAACGGTAAAAGAACAGGGATGGAAAATGCCAGCGGTTCAGTAATCCCCACCAGCATGGGGGTTAGGGTGGGGGAAAACAAAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGgcgc                                                                                           <  2:463885/239‑1 (MQ=255)
                                                                                                                                                         cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTGAAGGTTGGGGGGATCAACAAACCCGCCCCTTTTTTCCGAATTTCCGGGGCGGCAGTGGAAAACATCCCTAAAGAAAAGCACACGGCACCCAAGAATTTTTAGTTAACCGGGAGGGCGAAACAGCCCTACGCGGACAACGATATTCACGGCTCCgcgcgc                                                                                         >  2:2500554/1‑237 (MQ=255)
                                                                                                                                                                cGCAAACAGAAACGGTAAAAGAAAGAGGAAGTAAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGGGGTTAGGGTGGCGGAAATCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCa                                                                                                                                                                                                      <  2:544234/127‑1 (MQ=255)
                                                                                                                                                                cGCAAACAACAACGGTAAAAAGAACAGGAAGGAAATTCCCAGCGGTTCGGAAATCCCCCCCAGCAGGGCGGTTAGGGTTGCGGAAATCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGa                                                                                  <  1:463228/239‑1 (MQ=255)
                                                                                                                                                                 aaaaaCAGCACGGGAAAAAGGAACGGGAAGTAAAATCCCGGGGTTTCGGAAATTCCCACCAGCTTGGGGTTGAAGGTGGGGGGAATCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                 <  1:2500554/237‑1 (MQ=255)
                                                                                                                                                                   aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTGGGGGGATCACGAAGCCCCCCACTTTTACCCCATTTTTCGGTGCGGGAGTGAAGGACATCGCTAAAGAAAAGCCCACGCGACCCAAGAATTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCCGGAACAAAGATTTCTGCGCGGCGGCGCTCACACAGAAc                                                                                >  2:525475/1‑238 (MQ=255)
                                                                                                                                                                     aCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAACAGCAAGCCCGCCCCCTTTACCCGATTTTTCCGGGCGGCGGTGGAGGACATAGCTAAAAAAATACCCACGGCCCCCAAGATTTTTGAGTTACCGTGCAGGGCGAAAACGCCTTCCCGGAACAAACATTTCCTCGGGTCGGCGCTCACACAGacaac                                                                              >  2:2502053/1‑234 (MQ=255)
                                                                                                                                                                          aaCGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGTTGCGGGGATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGGAGTGAAGTACa                                                                                                                                                                                           >  1:544234/1‑127 (MQ=255)
                                                                                                                                                                          aaCGGAGAAAGGAACAGGAAGGAAAATCCCGGGGTTTCGATAATCCCCACCAGCATGGCGGTTAAGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  1:1090756/239‑1 (MQ=255)
                                                                                                                                                                            cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                      >  1:128769/1‑69 (MQ=255)
                                                                                                                                                                            cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGc                                                                                                                                                                                                                                                      <  2:128769/69‑1 (MQ=255)
                                                                                                                                                                               ggaaaaagaaaagGAAGAAAAATTCCGGCGGTTGGTTAATCCCCACCAGCAGGGGGGTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                                        <  1:27383/232‑1 (MQ=255)
                                                                                                                                                                                 aaaaTGAACGGGAAGGGAAATCCCAGCGGTTGGGTAATTCCCCCCAGCGGGGCGGTAAAGGTGGCGGGAATCAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                                   <  1:124248/237‑1 (MQ=255)
                                                                                                                                                                                     tGACCAGAAAGGAAAATCCCAGCGGTGCGGTAATCCCCACCAGCTTGGCGGTTAGGGTAGCGGAAAGCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:1061147/239‑1 (MQ=255)
                                                                                                                                                                                     tGAAGAGGAAGAAAATTCCCGGCGTTTCGGTAATCCCCACCAGCAGGGCGGTTAGGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2242422/239‑1 (MQ=255)
                                                                                                                                                                                       aaaaGGAAGTAAAATTCCGGCGTTTCGGTAATCCCCCCCAGCAGGGCGGTTAGGGTTGGGGGAACCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2040202/234‑1 (MQ=255)
                                                                                                                                                                                             gagggTAAATTCCAGGGGTACGGTAATTCCCACCAGCTGGGGGGTTAGGGTGGGGGGAACAAGAAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:2738873/229‑1 (MQ=255)
                                                                                                                                                                                              gaGGTAAAATCCCGGGGTTTCGTTAATCCCCCCCAGCGGGGCGGTTAGGGTGGGGGAAACCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGaa                                                                                 <  1:1626033/212‑1 (MQ=255)
                                                                                                                                                                                              aaGGTAATTCCCAGCGGTTCGGTAATCCCCCCCAGCGTGGCGGTTAGGGTGGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                              <  1:1834308/231‑1 (MQ=255)
                                                                                                                                                                                              aaGGTAAATCCCGGGGTTTCGGTAACTCCCCCCAGCATGGGGTTGAAGGTGGGGGAAACCAGCACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:2021630/234‑1 (MQ=255)
                                                                                                                                                                                                  tAAATCCCAGCGGTCCGTTAATTCCCCCCAGCAGGGCGTTTAAGGTGGGGGAAATAAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                                           <  1:2649972/230‑1 (MQ=255)
                                                                                                                                                                                                    aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTTTGGGAATCAAAAAAATCCCCAATTTTTACCGCTTTTCCg                                                                                                                                                                                                              >  1:920188/1‑85 (MQ=255)
                                                                                                                                                                                                       tCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAGCCCGCCCCTTTTACCCGGTTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGcc                                                                                                                                                                           >  1:584809/1‑116 (MQ=255)
                                                                                                                                                                                                           gCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCCCTTTTACCCGGTTTTTCGGTGCGGCAGGGAAGTACATCGCTAAAAAAAAGCCCACGGCACCAAAGATTTTTGa                                                                                                                                                     >  1:895088/1‑133 (MQ=255)
                                                                                                                                                                                                                               ccgtaaTGGGGGTTAAGTTTGCGGAAATCGAAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  1:1866628/234‑1 (MQ=255)
                                                                                                                                                                                                                               ccaGGAGGGCGGTTAAGGTTGCGGAAACCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  2:1742183/239‑1 (MQ=255)
                                                                                                                                                                                                                               ccaGCAGGGCGGTTAAGGTGGGGGGAACCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa                       <  1:135909/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ggtggTTAAGGTGGCGGAAAGCAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATa                 <  1:2080204/235‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTTAATGGTGGCGGAAAGCAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa               <  1:2614236/238‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaGGTGGCGAAAATCACAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTc          <  1:2095503/238‑1 (MQ=255)
                                                                                                                                                                                                                                                      gCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa                                                                                                                                                                                >  1:367584/1‑69 (MQ=255)
                                                                                                                                                                                                                                                       cGGGTATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGCACACCGCTAAAGAAATGCCCAAGGCACGAAAGATTTTTGAGTTAACGTGTAGGGCCAAACCGCCTACCGGGAACAACGATTTCAGCGGCTCGGCGCTCCAACTGAAATCTTGCAGATGCTGCCTCCAGGACACTTGACTGCCGCCTTCCACAGCTGGCGGACCAAAGCTAACCTGTCCGtt      >  2:984561/1‑234 (MQ=255)
                                                                                                                                                                                                                                                       cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                               <  1:848666/162‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAAct                                                                               >  2:848666/1‑162 (MQ=255)
                                                                                                                                                                                                                                                       cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:2096498/1‑239 (MQ=255)
                                                                                                                                                                                                                                                       cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATg  >  1:784119/1‑239 (MQ=255)
                                                                                                                                                                                                                                                       cGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACACCGCTAAAGAAATACCCACGGCACCAAAGATTTTTGAGGTACCGCGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACAACCGAAAGCCGCCTTCAAAAGCTCACAGGCACAAGCAAAACTGTACGTAGaat  >  2:709623/1‑237 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAGATG  >  NC_000913/3862428‑3862902

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