Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,638,0410TG53.8% 1.3 / 17.2 26V174G (GTG→GGG) PP_3201BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (8/4);  new base G (0/14);  total (8/18)
Fisher's exact test for biased strand distribution p-value = 3.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGA  >  NC_002947/3637907‑3638164
                                                                                                                                      |                                                                                                                           
ggcagCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTg                                                                                                                            >  2:457435/1‑136 (MQ=255)
 gcagCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTgg                                                                                                                           >  2:295316/1‑136 (MQ=255)
 gcagCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTgg                                                                                                                           >  1:240466/1‑136 (MQ=255)
 gcagCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTAGTCGGTGCCGACGGGGTgg                                                                                                                           >  1:499465/1‑136 (MQ=255)
                gCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACcgacg                                                                                                            >  1:15385/1‑136 (MQ=255)
                  aGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTAGGTGCCTACGGGGTGGCCAATCGCACCgacgac                                                                                                          >  1:58970/1‑136 (MQ=255)
                           tCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACAGATGCCCGCAACGGCCTGGTCCTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGgcgc                                                                                                 >  2:287319/1‑136 (MQ=255)
                              gcggggggcgccgCAACGGCTGCGGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTTAATTTCACCGATGCCCGCACCGGCCTGGTCGTCGGGGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGc                                                                                              <  2:267826/130‑1 (MQ=255)
                              ccggggggcggcgCAACGGCTGCGGGCCGACGGCGCGTACAAGCCGTTGTTGGCCCTGAATTTCACCGATCCCCGCAACGGCCTGGGCGTCGGTGCCTACGGGGTGCCCATGCGCACCGACGACGGCGGGCGCAGc                                                                                              <  2:43759/130‑1 (MQ=255)
                                         cgcAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTAGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGGGGGCGCCGCTGGCAGGCggg                                                                                   >  1:191095/1‑135 (MQ=255)
                                                        cctaCGGCGCTTACAAGCGGTTGTTGGCCCTCAATTACACCGAGCCCCGCAACGGCCGGGTCGGCGGTGCCTACGGGGTGGCCATGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                    <  2:15385/133‑1 (MQ=255)
                                                        ccGAGGGCGGGTACACGCGGTTGTTGGCCTTGATTTTCACCGATGCCCGCACCGGCCGGGGCGTCGGTCCCTACGGGGGGGCCAGGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                    <  2:209955/136‑1 (MQ=255)
                                                        ccGACGGCGCGGACACCCCGTGGTTGGCCCTGAATTTCACCGATCCCCGCAACGGCCTGGGCGGCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCTGACTGcc                                                                    <  2:105111/136‑1 (MQ=255)
                                                        ccGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCACCGGCCTGGTCGTCGGGGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                    <  1:473863/136‑1 (MQ=255)
                                                        ccGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGGCGGCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                    <  1:36143/136‑1 (MQ=255)
                                                        ccGACGGCGCGCACAAGCCGTTGTTGGCCCTCAATTTCACCGATCCCCGCAACGGCCGGGTCGGGGGGGCCTACGGGGGGGCCAGGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                                    <  2:334561/136‑1 (MQ=255)
                                                             ggCGGGTACAAGCGGTTGTTGGCCCTTATTTTCACCGATGCCCGAAACGGCGTGGTCGGCGGGGCCTAGGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACg                                                               <  1:25713/136‑1 (MQ=255)
                                                              gcgcGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGGCGGCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACgg                                                              <  1:105252/136‑1 (MQ=255)
                                                                       aGCCGTTGTTGGCCCTTAATTTCCCCGATGCCCGCAACGGCGTGGGCGTCGGTGCCTCCGGGGGGGCCATCCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACt                                                     <  2:140295/136‑1 (MQ=255)
                                                                                            acACCGAGCCCCCCACCGGCCTGGGCGGCGGTGCCTCCGGGGGGGCCATGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGgcgcg                                <  1:121970/135‑1 (MQ=255)
                                                                                              cccGATGCCCGCACCGGCGTGGGCGTGGGGGCCTAGGGGGTGGCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGcc                              <  2:205054/135‑1 (MQ=255)
                                                                                                 gATGCCCGCACCGGCCTGGTCGGCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAg                           <  2:415339/136‑1 (MQ=255)
                                                                                                   tGCCCGCAAGGGCGTGGTCGTGGGTGCCTCGGGGGGGGCCTTGCGCACGGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                         <  2:58970/136‑1 (MQ=255)
                                                                                                   ggCCCGCAACGGCCTGGGCGTCGGTGCCTACGGGGGGCCCATCCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                         <  1:217011/135‑1 (MQ=255)
                                                                                                                    gcgcgggTGCCTAGGGGGGGGCCTTGCGCACCTCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGc        <  2:240466/131‑1 (MQ=255)
                                                                                                                          gggCCTCCGGGGGGGCCCTGCGCCCCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGa  <  2:362303/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                           
GGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGA  >  NC_002947/3637907‑3638164

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