Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,8950AG29.6% 37.7 / 16.5 29intergenic (+189/‑357)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (16/3);  new base G (5/3);  total (23/6)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCC  >  NC_002947/5988762‑5989037
                                                                                                                                     |                                                                                                                                              
gcgcgGCGCTTCGCGGGTAAAACCCCAACAACCTTTTTTTAGGGCCAGGTATTACGGTGGGACGCCCCCCCCCGCCCCCAGGCGCACGCCACCAAAACGCGTCGTACAAACAACGCCGGCGCGCGCGCCCGGCCCAGGAGGGATTgg                                                                                                                                   >  8:73156/1‑147 (MQ=255)
               ggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCCCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTa                                                                                                                    >  3:42330/1‑147 (MQ=255)
                   aaCCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTGGTACCAACACGGCGGTCGCGCGCGCCTGGTGGAGGAGGGACTGGCC‑GAAACCAACGTAGGAGCg                                                                                                              >  7:49306/1‑148 (MQ=16)
                     cccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGGCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                                                                             <  1:8363/148‑1 (MQ=255)
                          tCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTa                                                                                                        >  7:155130/1‑148 (MQ=255)
                          tCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCTCGACCGCCTTGGGGCATGGCTCGATATCGCGTCGTACTAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAc                                                                                                       <  8:134915/148‑1 (MQ=255)
                           ccTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACACGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAACCGTAGGCGCGGGTTTAcc                                                                                                      >  5:47347/1‑148 (MQ=255)
                             tACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCcg                                                                                                    >  6:154602/1‑148 (MQ=11)
                                gtttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCGTCGATATCGCGGCGTACAAACAAGGCGGTCGCGCGCGCCTGGCCCGGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCcgc                                                                                                   >  7:118854/1‑147 (MQ=255)
                                  ttgttgCGGGCCAGTTCTGCCTGTGGGACGCGCACGCGAGCCCCCGGGCGCACGCCTCGATATCGCGGCGTACAAACAAGGCGGTCGCGCGCGCCGGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGaa                                                                                                <  7:73156/148‑1 (MQ=25)
                                            ccAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCcgcgcg                                                                                      >  2:121788/1‑148 (MQ=255)
                                                ttATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCgcgggcg                                                                                  >  4:121741/1‑148 (MQ=255)
                                                                     gcgACCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGACGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcgct                                                             >  1:46520/1‑146 (MQ=14)
                                                                        accgccttGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGgcg                                                               <  6:7271/135‑1 (MQ=14)
                                                                          cgccttGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGaaaat                                                       >  1:164776/6‑147 (MQ=14)
                                                                            cccTGGCGAACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCtaaaa                                                       >  1:183481/1‑147 (MQ=255)
                                                                             ccTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACt                                                     >  1:22290/1‑148 (MQ=255)
                                                                              ttGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGaaaatcctc                                                   <  2:164776/147‑6 (MQ=14)
                                                                                ggCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATc                                                  >  6:6131/1‑148 (MQ=255)
                                                                                         cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCAc                                         >  2:120176/1‑148 (MQ=255)
                                                                                              tatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGAAAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTg                                      >  2:46253/1‑146 (MQ=255)
                                                                                               atCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGGACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGaa                                   <  1:121788/148‑1 (MQ=255)
                                                                                                   cgtcgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCa                               >  1:112829/1‑148 (MQ=255)
                                                                                                    gtcgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCAt                              >  2:144790/1‑148 (MQ=255)
                                                                                                        tACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCt                           >  5:30481/1‑147 (MQ=255)
                                                                                                                          gcgcgcCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCAc        >  7:125298/1‑148 (MQ=255)
                                                                                                                          gcgcgcCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCAc        >  2:48324/1‑148 (MQ=255)
                                                                                                                              gcCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcg    >  4:106575/1‑148 (MQ=255)
                                                                                                                                cTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcgcc  >  2:114865/1‑148 (MQ=255)
                                                                                                                                     |                                                                                                                                              
GCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCC  >  NC_002947/5988762‑5989037

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