Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,906,6090GC20.0% 69.9 / 14.5 30intergenic (+26/+34)ygcG/enoTPM domain protein, putative phosphatase/enolase
Reads supporting (aligned to +/- strand):  ref base G (13/11);  new base C (3/3);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCA  >  NC_000913/2906475‑2906752
                                                                                                                                      |                                                                                                                                               
tAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCAtttt                                                                                                                                    >  7:42793/1‑148 (MQ=255)
   aaGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAttttttat                                                                                                                                <  2:154508/149‑4 (MQ=255)
     gTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATc                                                                                                                              <  2:61837/149‑1 (MQ=255)
        aTCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAAt                                                                                                                               >  1:23761/1‑145 (MQ=255)
                     caGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATa                                                                                                                         >  4:87199/1‑138 (MQ=255)
                     caGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATa                                                                                                                         <  3:87199/138‑1 (MQ=255)
                              cGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGc                                                                                                     <  8:42793/149‑1 (MQ=255)
                                            aTAGAGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                       <  4:20105/149‑1 (MQ=255)
                                                  gcggATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCg                                                                                 <  2:23761/149‑1 (MQ=255)
                                                                   ggCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTc                                                                <  3:143290/149‑1 (MQ=255)
                                                                       gTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAtttttta                                                                                                                                 >  1:83437/1‑78 (MQ=25)
                                                                       gTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCTCCGGGCTGGCAtttttta                                                                                                                                 <  2:83437/80‑3 (MQ=25)
                                                                         tcttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTAc                                                                                      <  6:109182/121‑1 (MQ=255)
                                                                         tcttctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTAc                                                                                      >  5:109182/1‑121 (MQ=255)
                                                                             ctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                       >  3:147218/1‑116 (MQ=255)
                                                                             ctGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTa                                                                                       <  4:147218/116‑1 (MQ=255)
                                                                                    tGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGa                                               >  2:169442/1‑149 (MQ=255)
                                                                                               ccGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                             <  2:113258/124‑1 (MQ=255)
                                                                                               ccGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGcc                                                             >  1:113258/1‑124 (MQ=255)
                                                                                                              tCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGc                      >  5:132862/1‑148 (MQ=255)
                                                                                                                 aCCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAaca                  >  3:75657/1‑149 (MQ=255)
                                                                                                                  gatttaaaaaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGCTACGAATCAGCTGGTTGTATTTAGCAacac                 >  3:131828/6‑149 (MQ=39)
                                                                                                                   cATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACg                >  7:75682/1‑149 (MQ=255)
                                                                                                                    aTAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGc                                                                      <  4:89684/94‑1 (MQ=255)
                                                                                                                    aTAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGc                                                                      >  3:89684/1‑94 (MQ=255)
                                                                                                                        aaaaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGt              >  5:146411/1‑146 (MQ=255)
                                                                                                                          aaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                         <  2:29511/101‑1 (MQ=39)
                                                                                                                          aaTGCCAGCCTCCGGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCaga                                                         >  1:29511/1‑101 (MQ=39)
                                                                                                                          aaTGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCCGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGaa         <  8:75682/149‑1 (MQ=255)
                                                                                                                                 gCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCa  >  1:72553/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                               
TAGAAGTTATCATTCAGACAACAGCGGTTCCGCTGGAGGCTCAGATAGCGGCGGATTTAGCGGCGGGGGCGGTTCTTCTGGAGGTGGCGGCGCATCCGGCCGCTGGTAACTCAACCATAAAAAATGCCAGCCCGGAGGCTGGCATTTTTAAATCAGATAAAGTCAGTCTTATGCCTGGCCTTTGATCTCTTTACGACCGTTGTACGGTGCTTTTTCGCCCAGAGCTTCTTCGATACGAATCAGCTGGTTGTATTTAGCAACACGGTCAGAACGGCTCA  >  NC_000913/2906475‑2906752

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