Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,990,8830TC34.4% 25.8 / 29.0 32intergenic (+199/‑347)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (6/15);  new base C (1/10);  total (7/25)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCC  >  NC_002947/5990753‑5991015
                                                                                                                                  |                                                                                                                                    
cgGGTAAACCCTCTCCTTCGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGGCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGcc                                                                                                                                 <  2:556253/136‑1 (MQ=255)
 gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCa                                                                                                                                <  1:567489/136‑1 (MQ=255)
                tACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGa                                                                                                                 >  1:669292/1‑136 (MQ=255)
                                           gtgGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcc                                                                                      <  1:670782‑M2/136‑5 (MQ=255)
                                               gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcg                                                                                  <  1:191472‑M2/136‑9 (MQ=255)
                                               gACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcg                                                                                  <  1:337024‑M2/136‑9 (MQ=255)
                                                aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGACGCGCCGcgcgggcgg                                                                                 <  1:331361‑M2/136‑10 (MQ=255)
                                                aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcgg                                                                                 <  1:780339‑M2/136‑10 (MQ=255)
                                                aCGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggccg                                                                                 <  1:173231‑M2/136‑10 (MQ=255)
                                                         cgAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgcccaatc                                                                        <  2:337440‑M2/136‑19 (MQ=255)
                                                         cgAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgcccaatc                                                                        <  1:756319‑M2/136‑19 (MQ=255)
                                                             cccccTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCCCAATCTCac                                                                    <  1:337440/136‑2 (MQ=255)
                                                                            cTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAAca                                                     >  2:742510/1‑135 (MQ=255)
                                                                             tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGCTCTCATGGGCGACCAAAAACTa                                                    >  2:565873/1‑136 (MQ=255)
                                                                             tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaaatcctc                                                   <  2:684192‑M2/136‑40 (MQ=255)
                                                                             tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgcgcgatctcacaggcgctgaaaatcctc                                                   <  1:71899‑M2/136‑40 (MQ=255)
                                                                                         cgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCCCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATAGCCAGGCAcc                                        >  2:197423/1‑136 (MQ=255)
                                                                                               cccaaggCGGTCGCGCGCAAATCCCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaactcttccggcgaacacctcccag                                  <  1:309701‑M2/134‑57 (MQ=255)
                                                                                               cccaaggCGGTCGCGCGCAAATCCCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaactctcccggcgaacacctcccag                                  <  2:112026‑M2/134‑57 (MQ=255)
                                                                                               cccaaggCGGTCGCGCGCAAATCCCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaactctcccg                                                 <  2:773694‑M2/119‑42 (MQ=255)
                                                                                                    ggCGGTCGCGCGCAAATCCCGCAGGAATAATTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaactctcccggcgaacacctcccaggcgct                             >  1:88072‑M2/1‑75 (MQ=255)
                                                                                                        gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGc                         <  2:670782/136‑1 (MQ=255)
                                                                                                            cgcgcgaCTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcacaggcgctgaaaat                                                       <  2:484939‑M2/101‑36 (MQ=255)
                                                                                                            cgcgcgaCTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcacaggcgctgaaaat                                                       >  1:484939‑M2/1‑66 (MQ=255)
                                                                                                                              gggATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcgc   >  2:34689/1‑136 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgctctaccaggcgctgaaacactaccggccaacacctaccaggcgctgaaaataccccggcgaacacctaccac  <  2:47190‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaactctcccggcgaacacctcccaggcgatgaaaatcacccgccgaacaccacccag  <  1:504925‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccaggcgctgaaacactcccggcgaacacctccccggcgcagaaaatcccccggcgaacacctccccg  <  2:203160‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcacaggcgctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccag  <  1:200765‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctagatctcacacgcgctgaaactctcccggcgaacacctcccaggcgctgaaaataccccgccgaacaccaaccac  <  2:623797‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcccgggcggcgctcgatcacacaggcgctgaaactcacccggcgaacacctcccaggcgcggaaaatcccccggccaacacctccccc  <  1:376659‑M2/136‑89 (MQ=255)
                                                                                                                              gggACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcccgcgcggcgctcgatctcccaggccctgaaactctcccggcgaacacctcccaggcgctgaaaatcccccggcgaacacctcccac  <  2:195944‑M2/136‑89 (MQ=255)
                                                                                                                                aatTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccagg                                                                 >  2:255152‑M2/2‑47 (MQ=255)
                                                                                                                                aatTGGCCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGcgcgggcggcgctcgatctcccagg                                                                 <  1:255152‑M2/71‑26 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGAGCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCC  >  NC_002947/5990753‑5991015

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: