Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,068,1940AT42.9% 1.0 / 12.3 14intergenic (‑96/+25)plsC/PP_18451‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base T (2/4);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTG  >  NC_002947/2068052‑2068340
                                                                                                                                              |                                                                                                                                                  
gtaaGAAAATGCAGCGCCATCAGGAACATACGCAATGAATAAAGCATGGTACCACTCACACCAGACGGTCGCGGGGCAGTGTACGTGTGTGCACTGGCTGGGGCAAACCTTGGTGCACGGCTGTTACGGCCCTTTCGCGGGCTTGCCcgc                                                                                                                                             <  2:112290/148‑1 (MQ=255)
    gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTGCGGGCTTGCCCGCGaaa                                                                                                                                         <  1:379004/150‑1 (MQ=255)
    gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGaaa                                                                                                                                         <  1:28398/150‑1 (MQ=255)
                       gAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGcc                                                                                                                      >  2:479887/1‑150 (MQ=255)
                                     aaTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGTGGGCTAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCACTATGCTGCTGAt                                                                                                        >  2:159147/1‑150 (MQ=255)
                                         aaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTCTACGGGTGTGCACTGGCTGTGGCAACCCTTGTTGCAGGCCTGTTCCGGCCCTTTCGCGGCCAAGCCCGCGAAAAGGGCCCGTGCCGCTAGCCAATATGCTGCTGATCCAg                                                                                                    <  1:253757/150‑1 (MQ=255)
                                          aGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGc                                                                                                   <  1:289252/150‑1 (MQ=255)
                                                         acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTcc                                                                                    >  1:274332/1‑150 (MQ=255)
                                                                         gggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGcc                                                                    <  2:29793/150‑1 (MQ=255)
                                                                                                                        cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTAGCTACGCTTGTCCAGCGCCCCCGAACGTCCTGGGCGGAGGTTCAGCGTGGTGTTCTTGTGCGCCAGCATGTTCAGc                     <  1:479887/150‑1 (MQ=255)
                                                                                                                                 cccTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCGGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTCTCAGCTTGGTGTTCTCCTGCGCCACCATGTTCCCCTTCTCCacc            >  2:172016/1‑148 (MQ=255)
                                                                                                                                  cttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                             >  1:558150/3‑67 (MQ=21)
                                                                                                                                  cttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                             <  2:558150/65‑1 (MQ=21)
                                                                                                                                           ggCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTg  >  2:497436/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTG  >  NC_002947/2068052‑2068340

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