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

CCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCT  >  NC_002947/5988760‑5989033
                                                                                                                                       |                                                                                                                                          
ccgcGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTg                                                                                                                                <  1:176090/148‑1 (MQ=255)
  gcgcgGCGCTTCGCGGGTAAACCCTCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGc                                                                                                                              <  4:42101/148‑1 (MQ=255)
           ttCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTACTGTGGGACGCGCCCGCGAGCCCCCAGGCGCACGCCGCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCCCAGGAGGGCTTGGCCCGaaaaaa                                                                                                                     >  2:37767/1‑148 (MQ=25)
              gcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaa                                                                                                                    <  5:98276/146‑1 (MQ=255)
                    aaaCCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAg                                                                                                            <  7:103507/148‑1 (MQ=255)
                       cccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                                                                         <  4:122903/148‑1 (MQ=255)
                       cccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCGCGATATCGCGTCGGACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                                                                         >  5:1332/1‑148 (MQ=255)
                        ccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCggg                                                                                                        >  8:154139/1‑148 (MQ=255)
                            tCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTCTCGATATCGGGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGTGCAGGAGGGACTGGCC‑GCAACAAACGTAGGAGCGGGTTTAc                                                                                                   >  5:133280/1‑148 (MQ=11)
                              cTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGGCCGCCTTGGGGGATGTGTCGATATCGCGTCGTACAAACAAGGCGGTGGCGCGCGTCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGGTTAccc                                                                                                 >  3:4615/1‑148 (MQ=255)
                              cTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGACCGCCTTGGGGCATGTGTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGGCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAccc                                                                                                 >  3:1947/1‑148 (MQ=11)
                               tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAccc                                                                                                 <  2:194727‑M2/148‑4 (MQ=255)
                               tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAccc                                                                                                 >  2:165144‑M2/1‑145 (MQ=255)
                                 cGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccg                                                                                                <  4:31943‑M2/147‑5 (MQ=255)
                                       tttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcg                                                                                         <  3:103141‑M2/148‑12 (MQ=255)
                                       tttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcg                                                                                         <  3:174469‑M2/148‑12 (MQ=255)
                                                 gTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGTGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggc                                                                               <  2:40026‑M2/148‑22 (MQ=255)
                                                               aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctca                                                                 <  7:182045‑M2/148‑36 (MQ=255)
                                                                cgcgCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCTCGAAGCGCCGCGCTGGCGGCGCTCGATCTCAt                                                                <  4:174469/148‑1 (MQ=255)
                                                                 gcgcACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCGCGATCTCATg                                                               >  5:132024/1‑148 (MQ=255)
                                                                     aCGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACagga                                                                                                                                       >  2:165761/1‑72 (MQ=255)
                                                                      cgcgAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                                                                      <  1:165761/72‑1 (MQ=255)
                                                                      cgcgAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGgcgc                                                          >  5:158056/1‑148 (MQ=255)
                                                                         gaccgccttggggCATGTCTCGATATCGCGTCGTACAACCAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAACCGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctga                                                       <  3:2197‑M2/142‑46 (MQ=255)
                                                                                    cgcACGCCTCGATATCGCGTCGTACAAACACGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATc                                              <  6:1332/146‑1 (MQ=255)
                                                                                     gcACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCATGGCACACGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCAGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGcc                                           >  7:59364/1‑148 (MQ=255)
                                                                                       atgtctcgaTATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaa                                         <  1:186674‑M2/143‑60 (MQ=255)
                                                                                          tctCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcacgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcct                                                <  2:156474‑M2/136‑53 (MQ=255)
                                                                                           cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctc                                               >  1:156474‑M2/1‑84 (MQ=255)
                                                                                                  tCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAg                              >  3:101216/1‑148 (MQ=255)
                                                                                                    gcgTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGcc                            >  2:164526/1‑148 (MQ=255)
                                                                                                      gtcgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCAt                          >  7:135102/1‑148 (MQ=255)
                                                                                                        cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATcc                        <  3:31943/148‑1 (MQ=255)
                                                                                                         gtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCt                       >  7:18726/1‑148 (MQ=255)
                                                                                                          tACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCt                       >  4:23351/1‑147 (MQ=255)
                                                                                                             aaacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCCCCTGGAAGCCATCCTGCaa                   >  1:21945/1‑148 (MQ=255)
                                                                                                               acaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtcc                <  2:14444‑M2/148‑85 (MQ=255)
                                                                                                                caaGGCGGTCGCGCGCGCCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacccctcaaaattccaacgccgtcca               <  5:62376‑M2/148‑86 (MQ=255)
                                                                                                                            gcgcgcCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCAc    <  1:194727/148‑1 (MQ=255)
                                                                                                                             cgcgaCTGGTGCAGGAGGGACTGGCC‑GAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaa  >  6:171927‑M2/1‑50 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCT  >  NC_002947/5988760‑5989033

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