Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,9050TC28.0% 70.3 / 33.3 50intergenic (+199/‑347)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (18/18);  new base C (5/9);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGC  >  NC_002947/5988761‑5989047
                                                                                                                                                        |                                                                                                                                              
cgcgcgGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTgg                                                                                                                                             <  3:67584/148‑1 (MQ=255)
     ggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGa                                                                                                                                        <  3:282196/148‑1 (MQ=255)
            cgcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaa                                                                                                                                   <  1:22378/146‑1 (MQ=255)
               gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaccaacag                                                                                                                                >  4:171093/1‑143 (MQ=255)
                ggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGCGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAgg                                                                                                                            <  7:214881/148‑1 (MQ=255)
                  tAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                                                                                             >  8:57541/1‑146 (MQ=255)
                       ccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCAAGCGA‑‑‑‑GCCACCTGGA‑‑GCACGCCTGGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCAAGAAACAAACGTAGGAgagcg                                                                                                                      >  1:71492/1‑144 (MQ=25)
                        cGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGCGTGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCgg                                                                                                                       >  1:259294/1‑148 (MQ=255)
                        cGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCggg                                                                                                                      >  5:11852/1‑147 (MQ=255)
                           tCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑‑CCGCCTTGGG‑‑GCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGGAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAc                                                                                                                 >  7:219508/1‑148 (MQ=11)
                                cGTTTGTTTCGGGCCATTTATTCCTGTGGGACGCGCAAGCGATCCTACCCACCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAACGCGGTCGCGCGCTCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTa                                                                                                                  <  1:111269/148‑1 (MQ=255)
                                  tttgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcga                                                                                                           <  7:269750‑M2/148‑8 (MQ=255)
                                    tgtttCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCGGCG‑‑GCAGGCCTCTATATCGCTTCTTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGTATTGGCCAGAAACAAACGTAGTAGCGGGTTTAcccgcga                                                                                                           >  6:263051‑M2/1‑139 (MQ=255)
                                                     tCCTGCGGGATTTGCGCGCGA‑‑‑‑CCGCCTTGGG‑‑GCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                       >  7:150316/1‑148 (MQ=14)
                                                         gtgGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctc                                                                                      <  1:262938‑M2/147‑29 (MQ=255)
                                                         gtgGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcga                                                                                    <  1:50628‑M2/148‑31 (MQ=255)
                                                          tgGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGACACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcga                                                                                    <  5:138533‑M2/147‑31 (MQ=255)
                                                             gACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatc                                                                                  <  5:7107‑M2/146‑33 (MQ=255)
                                                                     cgcgA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGgcgc                                                                        >  4:38987/1‑148 (MQ=255)
                                                                     cgcgA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTAGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGgcg                                                                         >  6:325573/1‑147 (MQ=255)
                                                                             accgccttggg‑‑GCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaa                                                                    <  2:15927‑M2/143‑47 (MQ=255)
                                                                               cccccTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCaaaa                                                                  <  2:262938/148‑1 (MQ=255)
                                                                               cccccTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCaaaa                                                                  <  2:50628/148‑1 (MQ=255)
                                                                                gccttGGG‑‑GCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAAAGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggctctgaaaat                                                                 <  6:264173‑M2/146‑50 (MQ=255)
                                                                                    cTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTa                                                              >  7:207485/1‑148 (MQ=255)
                                                                                               gtctcGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaaca                                                    <  8:233809‑M2/146‑63 (MQ=255)
                                                                                                    gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTg                                                >  3:113060/1‑148 (MQ=255)
                                                                                                     atatCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACACACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTCTCGCCAGGCACCt                                                 >  6:234807/1‑146 (MQ=255)
                                                                                                         cgcgTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGc                                           >  3:285632/1‑148 (MQ=255)
                                                                                                          gcgTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCGGCAAGc                                           >  2:296574/1‑147 (MQ=255)
                                                                                                           cgtcgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCa                                         <  7:57541/148‑1 (MQ=255)
                                                                                                                 acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaa                                             >  8:219413‑M2/1‑68 (MQ=255)
                                                                                                                 acaaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaa                                             <  7:219413‑M2/137‑70 (MQ=255)
                                                                                                                  caaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccg                                 <  8:171791‑M2/148‑82 (MQ=255)
                                                                                                                     acaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAAtt                               >  2:329265/1‑148 (MQ=255)
                                                                                                                        aGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAACCTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTcc                            >  7:240475/1‑148 (MQ=255)
                                                                                                                        aGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTcc                            >  1:188661/1‑148 (MQ=255)
                                                                                                                          gCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTg                          <  4:285632/148‑1 (MQ=255)
                                                                                                                           cGGTCGCGCGCGCCTGGCACAGGAGGGATTAGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTg                          <  5:325573/147‑1 (MQ=255)
                                                                                                                                 gcgcgcgaCTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaatt                    >  7:233809‑M2/1‑53 (MQ=255)
                                                                                                                                          tGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTcgcc            <  6:11852/147‑1 (MQ=255)
                                                                                                                                            gCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCgccgcc         <  8:207485/148‑1 (MQ=255)
                                                                                                                                            gCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCgccgcc         <  2:259294/148‑1 (MQ=255)
                                                                                                                                                aggaggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctggatctcacaggcgctgaaaatcctccggcaaacccctccaaattccaccggcgtcgcctccaaattccaccgccgcccccgc    <  2:308248‑M2/148‑111 (MQ=255)
                                                                                                                                                  gaggGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATAACAGCGCCGCCATGAc    >  6:170648/1‑147 (MQ=255)
                                                                                                                                                  gaggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctca   <  3:242638‑M2/147‑112 (MQ=255)
                                                                                                                                                   aggGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGAcgc  >  4:91739/1‑148 (MQ=255)
                                                                                                                                                   aggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctca   <  6:25842‑M2/146‑112 (MQ=255)
                                                                                                                                                   aggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgct                                                                                       >  6:80596‑M2/1‑35 (MQ=255)
                                                                                                                                                   aggGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgct                                                                                       <  5:80596‑M2/62‑28 (MQ=255)
                                                                                                                                                        |                                                                                                                                              
CGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑‑GCCCCCTGGC‑‑GCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCACCTCGCCGCCATGACGC  >  NC_002947/5988761‑5989047

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: