Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,828,5690CA24.0% 48.8 / 12.7 25intergenic (+53/‑52)srlD/srlMsorbitol‑6‑phosphate dehydrogenase/sorbitol=responsive srl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (9/10);  new base A (3/3);  total (12/13)
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.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCC  >  NC_000913/2828436‑2828710
                                                                                                                                     |                                                                                                                                             
ctgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAg                                                                                                                                 <  8:19504/148‑1 (MQ=255)
 tgctgTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGAGCGTAGCGCATCAgg                                                                                                                                >  4:117084/1‑148 (MQ=255)
           cGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGAGCGTAGCGCATCAGGCAttctat                                                                                                                        <  3:117084/146‑6 (MQ=255)
                     aaGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTa                                                                                                           >  7:94164/1‑149 (MQ=255)
                                  gCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                                               <  1:68511/116‑1 (MQ=255)
                                  gCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGc                                                                                                                               >  2:68511/1‑116 (MQ=255)
                                         aCAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATgg                                                                                        >  6:93273/1‑148 (MQ=255)
                                                          ccGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGAGCGTAGCGCATCAGGCAtt                                                                                                                            >  8:39387/1‑94 (MQ=38)
                                                          ccGGCGGTCAGGGGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGAGCGTAGCGCATCAGGCAtt                                                                                                                            <  7:39387/95‑2 (MQ=37)
                                                                             cTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgt                                                    <  5:93273/148‑1 (MQ=255)
                                                                                        cGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                                                       >  2:140146/1‑134 (MQ=255)
                                                                                        cGGAGATCCTTAATGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCt                                                       <  1:140146/134‑1 (MQ=255)
                                                                                                            ccGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgt                                                    <  3:128338/117‑1 (MQ=255)
                                                                                                            ccGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGgt                                                    >  4:128338/1‑117 (MQ=255)
                                                                                                                       taGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATc                                                           >  6:40353/1‑99 (MQ=255)
                                                                                                                       taGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATc                                                           <  5:40353/99‑1 (MQ=255)
                                                                                                                       taGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGtttt           <  8:94164/147‑1 (MQ=255)
                                                                                                                            tGCCTGATGAGCGTAGCGCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTg    >  1:174548/1‑149 (MQ=255)
                                                                                                                              ccTGATGAGCGTAGCGCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGcc  <  2:174548/149‑1 (MQ=255)
                                                                                                                               cTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggc                                <  2:24859/118‑1 (MQ=255)
                                                                                                                               cTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggc                                <  2:148190/118‑1 (MQ=255)
                                                                                                                               cTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggc                                >  1:24859/1‑118 (MQ=255)
                                                                                                                               cTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAggc                                >  1:148190/1‑118 (MQ=255)
                                                                                                                                tGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTg                           <  6:2325/122‑1 (MQ=255)
                                                                                                                                tGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTg                           >  5:2325/1‑122 (MQ=255)
                                                                                                                                     |                                                                                                                                             
CTGCTGTTCTACGCCAGTCCTAAGGCGTCGTACTGCACCGGACAGTCGATCAATGTCACCGGCGGTCAGGTGATGTTCTGATCAACAGCGGAGATCCATTAAGGATCTCCGTGAGACTATAGAATGCCTGATGCGCTACGCTCATCAGGCATACAGGACTTCCGCCACTACATTAAGGAAAAGTTATGGTATCCGCACTCATCACCGTCGCCGTTATCGCCTGGTGTGCGCAACTGGCCTTAGGCGGCTGGCAAATTTCTCGTTTTAACCGTGCC  >  NC_000913/2828436‑2828710

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