Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,7810GC20.7% 68.9 / 11.1 29intergenic (‑112/+45)maeA/sramalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme/stationary‑phase‑induced ribosome‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base C (4/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCA  >  NC_000913/1555659‑1555926
                                                                                                                          |                                                                                                                                                 
tttgGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCAtt                                                                                                                          <  5:33601/148‑1 (MQ=255)
     ttCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGAc                                                                                                                    >  1:11282/1‑149 (MQ=255)
                  actcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCaa                                                                                                                                            >  3:74510/1‑112 (MQ=255)
                  actcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCaa                                                                                                                                            <  4:74510/112‑1 (MQ=255)
                  actcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgt                                                                                                                                        <  6:10810/116‑1 (MQ=255)
                  actcTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGgt                                                                                                                                        >  5:10810/1‑116 (MQ=255)
                   ctctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGttttt                                                                                                                                  <  3:15469/121‑3 (MQ=39)
                   ctctTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCCTGCGCAATGGTGttttt                                                                                                                                  >  4:15469/1‑119 (MQ=39)
                             aaTATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAg                                                                                            >  7:157515/1‑149 (MQ=255)
                                   catccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGc                                                                                      <  8:157515/149‑1 (MQ=255)
                                    atccCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCg                                                                                     >  8:134897/1‑149 (MQ=255)
                                             gggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGcc                                                                            <  1:123848/149‑1 (MQ=255)
                                             gggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGcc                                                                            >  2:123848/1‑149 (MQ=255)
                                              ggCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCg                                                                           <  2:11282/149‑1 (MQ=255)
                                                          cTTTGC‑TTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg                                                                                                                       <  5:9898/92‑1 (MQ=255)
                                                          cTTTGC‑TTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAg                                                                                                                       >  6:9898/1‑92 (MQ=255)
                                                           tttGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                              >  7:28127/1‑149 (MQ=255)
                                                           tttGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTta                                                              >  7:20255/1‑149 (MQ=255)
                                                              gCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGa                                                                                                                 <  1:142444/95‑1 (MQ=255)
                                                              gCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGa                                                                                                                 >  2:142444/1‑95 (MQ=255)
                                                              gCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTtact                                                            <  8:20255/148‑1 (MQ=255)
                                                                                          tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                      <  6:190252/142‑1 (MQ=255)
                                                                                          tataAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCgg                                      >  5:190252/1‑142 (MQ=255)
                                                                                               aGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTggg                                                                     <  2:136459/106‑1 (MQ=255)
                                                                                                                acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCg                                       >  6:119216/1‑119 (MQ=255)
                                                                                                                acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCg                                       <  5:119216/119‑1 (MQ=255)
                                                                                                                acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa         >  7:10691/1‑149 (MQ=255)
                                                                                                                acacCATTGCCTGCGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTa         >  5:75068/1‑149 (MQ=255)
                                                                                                                   ccATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATgg      <  7:134897/149‑1 (MQ=255)
                                                                                                                        gcgcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGgtcca  <  8:28127/148‑1 (MQ=255)
                                                                                                                          gcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt                                                                        <  3:77441/76‑1 (MQ=255)
                                                                                                                          gcAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGtt                                                                        >  4:77441/1‑76 (MQ=255)
                                                                                                                          |                                                                                                                                                 
TTTGGTTCCATGTCACTCACTCTTTTTTGAATATCCATCCCTGGGGGGCTTTTATCGTCTTTGCTTTACCGCCAGGGCGTCGGCCTCAAGTATAAAGCAGATAAAAACAAAAACACCATTGCGCAGGCAATGGTGTTTAATCGTCATTGAGGACTGATGGTTATGAATTACTTTTCAGCGGGGCGTTTTCTGCCGGTTGGGTTATTTACTACGCTGGATTTGTCACCTTCGGTAACTATTTTGCGCTGGTTAGAAATTTTATGGTCCA  >  NC_000913/1555659‑1555926

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: ‑