Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,922,2040AC57.7% ‑5.3 / 11.3 26D271A (GAC→GCC) PP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base A (2/9);  new base C (2/13);  total (4/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCTTCGACC  >  NC_002947/5922085‑5922335
                                                                                                                           |                                                                                                                                    
ccGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGGCGGGCCCCGCACCGCCCGCGGCAAAgc                                                                                                                     >  1:2616547/1‑136 (MQ=255)
     tCGAGGGCCGGCTCCACACCGCCACCGCCTCCCCCGCCGGGGCCCCCCGGGCCCTCGCCCCCCGCCTCAAACCCCTGGGCCGCATGGACGTGGCAGAGACGCGCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCt                                                                                                                <  2:470338/136‑1 (MQ=255)
             ctgctgCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGGCCGGCGCCGGCCGGCCCGCCGCCAACCCCCCCAccccggg                                                                                                        >  2:3071880/1‑132 (MQ=255)
                ctgcACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGGCGGCCGCCGGCCCGCCCGCGGGCAAACCCCCCAACCCgcgggg                                                                                                      >  1:1239387/1‑131 (MQ=255)
                ctccaccCCGTCACCGCCTGCCCCGCCGGGGTCACCCGGGCCCTCGCCCGCCGCCTCAAACCCCTGGGCCGCTGGCCGGTCGCTGACACGCGCCGGCCGCAGGCGGGCCGGCTGCAAAGCCGCCTACCCGgtggtg                                                                                                     <  2:2389585/133‑1 (MQ=255)
                           ctccGCCTCCCCCGCGGGGGCCCCCGGCCCCCGCGCCCGCCCCCTCAACCACGGGGGCCGATGGGCGGTTCCAGCGACGCGCCGCCCGCAGGCGGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAg                                                                                          <  2:1026025/133‑1 (MQ=255)
                           cacccccTCCCCCGCCGGGGACACCCGGCCCCTCGCCCCCCCCCTCAAACCCGGGGCCCGCTTGGCGGTTGACCAGACCCGCCGGCCCCCGGCCGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAg                                                                                          <  2:485049/131‑1 (MQ=255)
                           cTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGCCCGCCGCAGGACGGCCCGCAGCAAAACCGCCAACCCCCGGGTGGGGGAGTgggg                                                                                          >  2:3102616/1‑134 (MQ=255)
                           cTACGCCTTCCCCGCCGGGGTCACCCTGACCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGCCGGCCACAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAg                                                                                          <  1:45924/136‑1 (MQ=255)
                             ccGCCTCCCCCGCGGGGGCCCCCCGGCCCCGCGCCCGCCGCCTCAACCCCGGGGGCCCCTGGCCCGTCGCACCTACGCCCCGCCCCCAGCACGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCt                                                                                        <  2:261224/135‑1 (MQ=255)
                                  ttCCCCGCGGGGGTCACCCTGGCCCTCGTCAGCCCCCTCAACCACCTGGGCCGAATGTACGTTCCAGATAAGCGCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGc                                                                                   <  2:2598488/136‑1 (MQ=255)
                                  ccccccGCCGGGTTCACCCTGCCCCGCGGCCCCCGCCTCAACCCCCTGGCCCGCATGGACGTTCCAGATACCCGCCGGCCGCCGGAGGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGc                                                                                   <  2:350710/134‑1 (MQ=255)
                                       ggccgGGGTCCCCCGGCCCCGCCGCCGCCGCATCACACACGTGGGCCGCTTCGAGTGTCCACAGACCGCGCGCCCGCAGCACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCg                                                                              <  2:254672/135‑1 (MQ=255)
                                       cgcggGGGCCCCCCGGCCCCTCGTCCCCCGCCTTAAACACCTGGCCCGCATGGACGTTGCAGACACGCGCCGGCCGCAGGACGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCg                                                                              <  2:1239387/136‑1 (MQ=255)
                                       ccggggggTCCCCCGGGCCCGCGTCCGCCGCCTCAACCCCGGGGCCCCCTGGGTGGTTCCAGACACGCCCGGCCGGCCGCCGGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCg                                                                              <  2:1750804/132‑1 (MQ=255)
                                        gccgGGGTCCCCCGGCCCCTCGCCCGCCGCCTCAACCCCCGGGGCCGCAGGGACGTTGCAGAGACGCGCCGGCCGCAGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGa                                                                             <  1:4214908/136‑1 (MQ=255)
                                                          cTCGTCAGCCGCCTGAAACACGTGGGCCGCATGGACGTTGCAGAGAAGCGCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTg                                                           <  2:1820213/136‑1 (MQ=255)
                                                                          aaCACCTGGGCCTCTGGGCGGTGGCAGACCCGCGCCGGCCGCTGGCGGGCCGGCTGCTAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGtt                                           <  2:2115882/136‑1 (MQ=255)
                                                                                ggggCCGCAGGGACGTTGCACCGACGCGCCGCCCGCAGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCt                                     <  1:3102616/135‑1 (MQ=255)
                                                                                          ggcccGTTGCAGACACGCGCCGGCCGCCGGCGGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGAccc                            <  2:2355628/132‑1 (MQ=255)
                                                                                                tttGCAGTCACGCGGCGGGCGCTGGCGGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGATGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCAcc                       <  1:4390771/135‑1 (MQ=255)
                                                                                                          cgcgcCGGCCGCTGGCCGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAaga             <  1:1478657/135‑1 (MQ=255)
                                                                                                                 ggccgccggccgccCGGCTGCAAAGCCGCCTACCCGGTGGGGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCtt       <  1:2885689/136‑1 (MQ=255)
                                                                                                                 ggccgCCGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCtt       <  2:3201947/136‑1 (MQ=255)
                                                                                                                 gcccgCAGGAGGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCtt       <  2:2897542/134‑1 (MQ=255)
                                                                                                                       cAGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCTTCGAcc  <  2:1598804/136‑1 (MQ=255)
                                                                                                                       cAGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCTTCGAcc  <  1:4333100/136‑1 (MQ=255)
                                                                                                                           |                                                                                                                                    
CCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGGCCGCATGGACGTTGCAGAGAAGCGCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGCAAGAGGGCTTCGACC  >  NC_002947/5922085‑5922335

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