Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,951,7660TG50.0% ‑3.8 / 10.3 28intergenic (+66/+47)PP_2583/oguAacetyltransferase/8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base T (6/8);  new base G (6/8);  total (12/16)
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
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTT  >  NC_002947/2951634‑2951899
                                                                                                                                      |                                                                                                                                     
aTACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTgcg                                                                                                                                    >  2:1740444/1‑136 (MQ=255)
     tGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCCCCCCCAAAGGGGGGCGCagg                                                                                                                                >  2:2574418/1‑134 (MQ=255)
      gAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCCGCCCCAAGGGCGGGGCCTGcc                                                                                                                              >  2:736842/1‑136 (MQ=255)
      gAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGgcggcggcttcc                                                                                                                              >  2:1808772/1‑130 (MQ=255)
          ccgtcgccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGACCAAAGGGCGGCGGATGcccgtg                                                                                                                          >  2:1605604/1‑133 (MQ=255)
               gccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCttcattca                                                                                                                     >  2:2726008/1‑136 (MQ=255)
                 cgtcgCACCCCGGCGGAAGGGCCTGCTTTGAGTTGGCATTGCCCCTTTGTGTGGACATTAATTGTTTGGGGGCGCTACCCCCTGCTTGCCTTGGTGCCAGCGCAGCCCAACGGGCTGCGCATGCCCTTCATTCAcc                                                                                                                   <  2:1631421/136‑1 (MQ=255)
                  gtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCCCCCCAAAAGGGTGGGCCTGGCCGTAATTAAccc                                                                                                                  >  2:2983901/1‑136 (MQ=255)
                        ccAGGGCGGAAGGGCGGGGCATGAGTGGGGCTTGCCCCTTTGTGTGGAAATTACATGGATGTGTGCCCTCCCGGCAGCTTGCCTCGGGGGCAGCGCAGCCAAATGGGCTGCGCATGCCCTTCATTCACCCGGCGAt                                                                                                            <  2:2084239/136‑1 (MQ=255)
                        ccACGGCGGCTTGGCCGGGCTTGTGTTGGTCATGGACGTTTGTGTGGACATGACATGGTTGTGTGCGCTCCCCGCTGCTTGCCTTGGGGGCTGCGCAGCCCACCGGGCGGCCCATGCCCTTCATTCACCCGGCGAt                                                                                                            <  2:2336320/136‑1 (MQ=255)
                                  tgggCCGGGCTTGAGTTGGCCTTGCACTTTTGTGGTGACATCACATCGTTGGGGGCGCACCCCGCAGCTGGCCTTGGGGCCCCCGCAGCCAAATGGGCTGGGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGa                                                                                                  <  2:2417744/134‑1 (MQ=255)
                                             tGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCCGCCCAAAGGGGGGGGCATGCCCGTAATTAACCCGGCGCTGCCCAGCTGATGCCCGGccgc                                                                                       >  2:2756206/1‑134 (MQ=255)
                                                               ttGCGTGTACTTTACATGGTTGTGTGCGCTACCCGCTGCTTGCCTCGGGGGCCGCGCAGCCAAACGGGGGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCt                                                                     <  2:313152/136‑1 (MQ=255)
                                                                      aaatttacatCGTTGGGGGCGCCCCCCGCCGCTGGCCTCGGGGGCCGCCCAGCCAAACGGGGGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                              <  1:141899/130‑1 (MQ=255)
                                                                          tgacatCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAACGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGc                                                          <  2:1639321/136‑1 (MQ=255)
                                                                               tCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCCCCCCCAAGGGCGGCGCCAGCCCATGATTAAACCCGCGCATCCCAGCAGATGCCACGCCACAACTTCCACAAACCAGCAGATGGGCCAGCGGGtt                                                     >  2:3109659/1‑136 (MQ=255)
                                                                               tCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                     <  2:2492916/136‑1 (MQ=255)
                                                                               tCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGAAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                     >  2:98941/1‑136 (MQ=255)
                                                                               gcGTAGTGGCCCCTACCCGCAGCTTGCCTAGGGGGCAGCGCAGCCCAACGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                     <  2:2418435/135‑1 (MQ=255)
                                                                                                    gCTTGCCTAGGTGGCAGCGCCGCCCCAAGGGCGGCGCCTGCCCGTAAATAACCCCGCGCTGCCGCGCTGATGCCAGGGCACAAGTTCCACAGACCAGCGGAAGGGGCCGCAGGGTGTTGCGGTGCAGTAAGCGCgg                                >  2:1617198/1‑136 (MQ=255)
                                                                                                    gCTTGCCTAGGTGGCAGCGCCGCCCAAAAGGGTGCGCCTGCCCGTAATTAAACCCGCGCTGCCGCGCTGATGGCAGGCCGCAAGTGCCACAGACCAGCAGCTTGGGCCGCAGGTTTTGGGGTTGCAGTACGCGCgg                                >  2:798306/1‑136 (MQ=255)
                                                                                                                        ccccccAACGGGCGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGc             <  1:2197029/133‑1 (MQ=255)
                                                                                                                          gtccAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg           <  2:3161262/134‑1 (MQ=255)
                                                                                                                          ccccAACGGGCGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg           <  2:1650888/135‑1 (MQ=255)
                                                                                                                          ccccAACGGGCGGCCCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg           <  1:2161697/135‑1 (MQ=255)
                                                                                                                          ccccAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg           <  1:1808772/135‑1 (MQ=255)
                                                                                                                           cccAAACGGGCGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTg           <  1:736842/136‑1 (MQ=255)
                                                                                                                            ccaacggggCGGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTgg          <  2:1457293/130‑1 (MQ=255)
                                                                                                                                    gCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCttt  >  2:2744856/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                     
ATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTT  >  NC_002947/2951634‑2951899

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