Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,4940TA69.2% 14.6 / 16.0 27intergenic (+665/+71)ycdU/serXuncharacterized protein YcdU/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base A (4/14);  total (5/22)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTCGTTCCTGTCTTGCTGGAACGGGCGC‑TAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAA  >  NC_000913/1097259‑1097728
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
gCGTTCTCGGTCTGTCCAACTGGCTGGCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTGTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCgg                                                                                                                                                                                                                                   <  1:483176/250‑1 (MQ=16)
        ggtttgtgcAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTGCCCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCaa                                                                                                                                                                                                                          <  2:371695/243‑1 (MQ=16)
             gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAata                                                                                                                                                                                                                        <  1:178169/248‑1 (MQ=33)
               ccAACTGGCTGTCGCCATTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACg                                                                                                                                                                                                                    <  1:347938/250‑1 (MQ=18)
                            gCCAGTTGTCGAACCCCTGTCGGGGCGTCTCATCCCCCCGGTGTGTGAAATATACGAAAAAAAAGCCCGTACGTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGt                                                                                                                                                                                                      <  1:111620/251‑1 (MQ=16)
                             ccAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTGAA‑ACGCTTCACTCGCCCCGCGGGGACCCCCCCCCCTTCGGTCACTGTCTGTTCAACTTGCTGTCGCCAGGTTTCGAACCCCCGTCCGGGCAGCTCAATCCCCCCG‑TGTGTTGAATATACGaaaaaaaaaaccct                                                                                                                                                                                                      >  1:347096/1‑244 (MQ=2)
                                     cGNACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCGTCGCTCGCCCTTCGGGCAGCCCGCTCACTGCGGTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGt                                                                                                                                                                                              <  2:475734/250‑1 (MQ=12)
                                     cGAACCCCGGTCGGGGCTTTTCAGCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTCCTTTCGTACGAGCTCTTCTTTAAATATTGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGt                                                                                                                                                                                              <  2:377439/250‑1 (MQ=16)
                                         ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGTAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAg                                                                                                                                                                                         <  1:521503/251‑1 (MQ=14)
                                         ccccGGTCCGGGCGTCTTATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAg                                                                                                                                                                                         <  1:214203/251‑1 (MQ=9)
                                              gTCGGGGCTTCTCGTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCtctt                                                                                                                                                                                    <  1:439362/251‑1 (MQ=14)
                                                  gggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCttt                                                                                                                                                                                <  1:177935/251‑1 (MQ=14)
                                                    ggttcgtaTCCCCCCTGGGTGTTCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCTATCTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCGCA‑TCCCCCCGGGGTGTGCTATAGACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTaa                                                                                                                                                                              <  1:504074/244‑1 (MQ=11)
                                                       tctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTGCGTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCTCTGCGTTCACGGTCTGTCCAACTGGCGGTCGCCAGTTGTCGAACCCCGGTCGTGGCATCTCA‑TCCCCCCGGTGTGTGCAAGATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAat                                                                                                                                                                            <  1:465767/250‑1 (MQ=11)
                                                          cgtCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATgg                                                                                                                                                                        <  1:449158/249‑1 (MQ=14)
                                                          cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑GCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAatat                                                                                                                                                                          <  1:156692/249‑1 (MQ=14)
                                                             ccccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGTGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCGCGCTCACTGCTTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCgg                                                                                                                                                                     <  2:347096/251‑1 (MQ=14)
                                                               ccccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTACCTC‑CTTCGCCCGCCCCGCGGGGAGCCCGCTCACTGCGTTTCCGGTCTGTCCCACCGGCTGTGGCCCCGTGTACAACACCCGTCAGGGCCTCTCC‑TCCCCCCGGGTTGTTGAATATACGAAAAAAAAAACCGTACCTTTGTTAGAGCTCTTTTTTacatttttggggg                                                                                                                                                                   >  1:59720/1‑239 (MQ=2)
                                                                                               cGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCg                                                                                                                                 <  1:386412/251‑1 (MQ=39)
                                                                                                     ttCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGGC‑TCGCTTCGCTCGCCCTGCGGGCAGCCTGTTTACTGCGTTCACGGGCTACCCAACTTGCTGTCGCCCGTTGTCGAACCCAGGTCGGGGCTTCTCA‑CCCCCTCTGCGTGTGCAATATTCGAGAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGACCCCCC‑GATACGTTGCCGtatcc                                                                                                                            <  2:427881/250‑3 (MQ=14)
                                                                                                      tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATaca                                                                                                                           <  1:420914/250‑1 (MQ=255)
                                                                                                                               gCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCGTAACTTGGCCGT‑TACACA‑CTTTTCAGACGTGCTccctttgcc                                                                                                >  1:84163/1‑244 (MQ=11)
                                                                                                                                       ggggggATTGAC‑TCGCTTCGTTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCGCGGTCTGTCCAAGTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑GCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATACACT‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGac                                                                                        <  1:491506/251‑1 (MQ=255)
                                                                                                                                                                                 ccACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTCGTTCCTGTCTTGCTGGAACGGgcgc                                              <  2:333749/250‑1 (MQ=255)
                                                                                                                                                                                                                    cAGTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGGTTGGACCCCCC‑CA‑ACGTTTCCGTATACACAACTTTCCAG‑CGTGCTCCTTCAGCCCCCCGGGCACCCCCCCCAAATT‑TTGTTTCTTTCTTTCTGGGACGGGCGCCAAATTTTGGGAAAT‑ATAACCTGAgggtcacaaa           >  1:157871/1‑242 (MQ=2)
                                                                                                                                                                                                                      gTTGTCGAACCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACATCTCACC‑AAATTGTCGTTCCGGTCGGGCTGGAACGGGCGC‑TAATTTAGGGAAATCATGACCTGAGGTCAACAAACtt        <  1:312349/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaCCCCGGTCGGGGCATCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTGGAACACCC‑CATACGGTGCCGTATACACA‑CTTTCCCGGCGTGCTCCTTCAGCCCCCCCGCCCCCTCACC‑CAAATGTCGTTCCTGTCTTTCTGGAACGGGCGC‑TCATTTATGGAAATCATTACCTGAGGGCACAAAACATTTTGat  >  2:419571/1‑249 (MQ=18)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                
GCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAC‑TCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCA‑TCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCC‑GATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑AAATTGTCGTTCCTGTCTTGCTGGAACGGGCGC‑TAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAA  >  NC_000913/1097259‑1097728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑