Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,294,6370TA53.8% 8.7 / 20.0 26T652T (ACT→ACAPP_2893hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/10);  new base A (14/0);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCTGTGT  >  NC_002947/3294503‑3294760
                                                                                                                                      |                                                                                                                           
aCGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGactggcccctc                                                                                                                            >  1:428270/1‑135 (MQ=255)
     gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCTCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACAGGCCACAGGCCgc                                                                                                                       >  1:979256/1‑135 (MQ=255)
     gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCgc                                                                                                                       >  2:687601/1‑135 (MQ=255)
     gCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCGa                                                                                                                       >  1:557644/1‑136 (MQ=255)
       aGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCgcaa                                                                                                                     >  2:705684/1‑133 (MQ=255)
       aGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGCCGaaa                                                                                                                     >  2:679304/1‑135 (MQ=255)
          ccAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCa                                                                                                                  <  2:606457/136‑1 (MQ=255)
            aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg                                                                                                                <  2:236439/136‑1 (MQ=255)
            aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg                                                                                                                <  2:162693/136‑1 (MQ=255)
            aaCGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAgg                                                                                                                <  2:686377/136‑1 (MQ=255)
                 cAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCTGGGCGCATACCGGcccgg                                                                                                           >  1:49406/1‑133 (MQ=255)
                 cAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTCGCCTGGATGGCCAGTTGATCAAAGGTGACTGGCGCGGCCGAGAGGCCACAGGCCGAATCCAGGCCggg                                                                                                           >  2:231665/1‑135 (MQ=255)
                      gCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACATGCCACACGCCGCAACCCCGCCGGTgggcc                                                                                                      >  2:400002/1‑133 (MQ=255)
                                       gCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACGGGCCCCAGGCCGAATACACGCCCGGGGCCAGGAGTGGCAGTTTCAgg                                                                                     >  1:320683/1‑136 (MQ=255)
                                          gtGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGAGGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGgcgc                                                                                  <  1:715816/135‑1 (MQ=255)
                                            gACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGAGAGGCCACAGGCCGAATACACGCCGGTGGGCCGGCGTGGCAGTTTCAGGcgccc                                                                                >  1:199021/1‑135 (MQ=255)
                                                   gcgcGGGACGGCCCGCAGCTCACGCGCGACCTTGGCCTGGAGGGCCAGTTGAACAATGGTGCCTGGCGCGGCCGACTGGCCACTGGCCGAATCCTGGCCGGTGGCCTGGAGTGGCTGTTGCAGGcgcctgcacgcc                                                                         <  2:591378/136‑1 (MQ=255)
                                                    cgcTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGAGGGGCCCCAGGGCGCATACCGGGCGGGGGGCCGGCGTGGCAGTTGGCGGcgccggcccgccc                                                                        >  1:728770/1‑135 (MQ=255)
                                                     gcTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCCCAGGCCGCATCCCGGCCGGGGGGCAGGCGTGGCAGGTGCAGGcgccggcccgccgg                                                                       >  1:282786/1‑134 (MQ=255)
                                                               ccGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTgg                                                             <  1:761031/136‑1 (MQ=255)
                                                                cGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCTACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTggc                                                            >  1:310568/1‑136 (MQ=255)
                                                                cGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTggc                                                            >  2:1005955/1‑136 (MQ=255)
                                                                             cGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACAGGGCGAATACACGCCGGGGGGCAGGAGTGGGAGTTTGCGGCGCCTGCCCGCCGGGCGCGGCGGGCCCGCGGGCCGCTc                                               >  1:890581/1‑136 (MQ=255)
                                                                                                   ttGAACAACGGTGACGGGCGCGGCCGACTGCCCACTGGCCGAATCCCGCCCGGTGGCCAGTAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGgc                         <  1:705684/136‑1 (MQ=255)
                                                                                                   ttGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAAACCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGgc                         <  1:679304/136‑1 (MQ=255)
                                                                                                              tGACTGGCGCGGCCGACGGGCCACAGGCCGAATACACGCCCGTGGCCAGGAGTGGCAGTTTCAGGCGCCTGCACGCCGGCAGCCGCAGGCCCGCCGCCACCGCGACTTCGGGGCCCCCATCAGGCGCTACCGcccc              >  2:547571/1‑134 (MQ=255)
                                                                                                                          ccGACGGGCCACTGGCCGAAGCCAGGCCGGGGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCtgtgt  <  1:440585/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                           
ACGTTGCAGGCCAACGGCAAGGGCGACATTCGCCAACAGGCCTTGACCCTGGCGCTGGACGGCCCGCAGCTCAAGCTCGACCTTGGCCTGGATGGCCAGTTGAACAAAGGTGACTGGCGCGGCCGACTGGCCACTGGCCGAATCCAGGCCGGTGGCCAGGAGTGGCAGTTGCAGGCGCCTGCACGCCTGCAGCGGCTGGCCAGCGGCCAGCTCGACTTCGGGGCCCACTGCTGGCGCTCCGGCCAGGCAAGCCTGTGT  >  NC_002947/3294503‑3294760

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: