Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294756,2960T.25.0% 37.7 / 10.8 16coding (1332/1467 nt)PP_0049hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base . (4/0);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 8.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAGGACTGGCTGCTGTCAGCAACTCATCGAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTCAGCAGACCACGGCTCGTAGCCATATTGCGGCCAGAGCTGACCGTAGT  >  NC_002947/56172‑56423
                                                                                                                            |                                                                                                                               
gaggaCTGGCTGCTGTCAGCAACTCATCGAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGTGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTg                                                                                                                      <  2:775359/136‑1 (MQ=255)
       ggCTGCTGGAAGAAACTCATCTAGTGACGGAACGTTCGGGTGTGCACTGAGGACGAGTGTTCTGGATCGTAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACt                                                                                                               <  1:389040/136‑1 (MQ=255)
                         aTCGAGTGACGTAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTa                                                                                             <  1:451478/136‑1 (MQ=255)
                         aTCGAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTa                                                                                             <  1:102217/136‑1 (MQ=255)
                         aTCGAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTa                                                                                             >  1:178516/1‑136 (MQ=255)
                            gAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGc                                                                                          <  1:373983/136‑1 (MQ=255)
                                gACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCg                                                                                      >  2:600119/1‑136 (MQ=255)
                                   ggAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGTATACATCAAATAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGt                                                                                   <  2:797893/136‑1 (MQ=255)
                                         ttAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTg                                                                             <  1:236697/136‑1 (MQ=255)
                                                          gACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAAt                                                            <  1:207015/136‑1 (MQ=255)
                                                          gACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAAt                                                            <  1:440501/136‑1 (MQ=255)
                                                                     tGGATCGGAGTGTGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTA‑‑‑CTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTcagc                                              >  2:661470/1‑136 (MQ=255)
                                                                     tGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTA‑‑‑CTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTcagc                                              >  2:210608/1‑136 (MQ=255)
                                                                     tGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTA‑‑‑CTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTcagc                                              >  2:806177/1‑136 (MQ=255)
                                                                     tGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTA‑‑‑CTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGCTATAGTGCCATAGGGTTATGGAATACCTCCACTTcagc                                              >  2:665258/1‑136 (MQ=255)
                                                                                                                    ccGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTCAGCAGACCACGGCTCGTAGCCATATTGCGGCCAGAGCTGACCGTAGt  <  1:110343/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                               
GAGGACTGGCTGCTGTCAGCAACTCATCGAGTGACGGAACCTTAGCGTTTGCACTGAGGACGAGTGTTCTGGATCGGAGTATGGATACATCAAAGAAAGATTCGCACTCGACTCTTCCGTTTACTTCTCTCCAGTGGGTGACTTCTGGTATCAACTCATTAGGCACCGGGTATAGTGCCATAGGGTTATGGAATACCTCCACTTCAGCAGACCACGGCTCGTAGCCATATTGCGGCCAGAGCTGACCGTAGT  >  NC_002947/56172‑56423

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: