Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,7830AG25.0% 41.6 / 13.8 28intergenic (‑114/+43)maeA/sramalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (3/4);  total (14/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 = 4.24e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTC  >  NC_000913/1555649‑1555929
                                                                                                                                      |                                                                                                                                                  
tttttttGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtttt                                                                                                                                      <  3:31591/149‑2 (MQ=255)
         ttttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                             >  2:37707/1‑149 (MQ=255)
                 ccATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACt                                                                                                                      >  6:154107/1‑148 (MQ=255)
                   aTGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGAt                                                                                                                   >  5:92521/1‑149 (MQ=255)
                      tcactcacTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATgg                                                                                                                 <  3:57080/148‑1 (MQ=255)
                         ctcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtt                                                                                                                                        >  1:123057/1‑122 (MQ=255)
                         ctcactcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGtt                                                                                                                                        <  2:123057/122‑1 (MQ=255)
                             ctctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAAt                                                                                                         >  4:106706/1‑149 (MQ=255)
                                                     gggggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTc                                                                                   <  1:37707/147‑1 (MQ=255)
                                                              tATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTgg                                                                         <  5:154107/148‑1 (MQ=255)
                                                                  gTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTAt                                                                    <  3:106706/149‑1 (MQ=255)
                                                                                     ggCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCa                                                                                                <  6:47569/102‑1 (MQ=255)
                                                                                     ggCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCa                                                                                                >  5:47569/1‑102 (MQ=255)
                                                                                     ggCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGt                                                  >  7:18731/1‑148 (MQ=255)
                                                                                              cTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGcc                                                                               >  3:141593/1‑110 (MQ=255)
                                                                                              cTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGcc                                                                               <  4:141593/110‑1 (MQ=255)
                                                                                                 aaGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTaa                                      <  8:18731/148‑1 (MQ=255)
                                                                                                   gTAAAAAGAAGATAAAAACAAAAACACCAATACGCAGGCAATAATGTGGAAAACGCATAGAGTACTGATGTTTATGAATTACATTACAGCGTGGCGTTGTCTGCCGTTTGGGTTAATTACTAAGCTGGATTTGTCACCTTAGGTAACTa                                   >  4:96276/1‑149 (MQ=255)
                                                                                                    tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTc                                                                                   >  1:152919/1‑100 (MQ=255)
                                                                                                    tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTc                                                                                   <  2:152919/100‑1 (MQ=255)
                                                                                                                   aacaaaaacaCCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAg                    >  2:85285/1‑148 (MQ=255)
                                                                                                                        aaacaCCATTGCGCAGGCAATGTTGTTTAATCGTTATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAAttt              <  3:96276/149‑1 (MQ=255)
                                                                                                                         aacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCggg                                                                                           >  7:11384/1‑71 (MQ=37)
                                                                                                                         aacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCggg                                                                                           <  8:11384/71‑1 (MQ=37)
                                                                                                                         aacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGtt                                                                      <  2:51807/92‑1 (MQ=38)
                                                                                                                         aacaCCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGtt                                                                      >  1:51807/1‑92 (MQ=38)
                                                                                                                          acacCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa            <  1:85285/149‑1 (MQ=255)
                                                                                                                                     cAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtccagtc  >  5:106208/1‑148 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
TTTTTTTGTTTTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCAGTC  >  NC_000913/1555649‑1555929

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