Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947284,5870TA42.9% ‑1.6 / 10.3 14intergenic (+33/+30)betT‑II/atsKcholine/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base A (5/1);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 1.38e-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.

AACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCC  >  NC_002947/284459‑284734
                                                                                                                                |                                                                                                                                                   
aaCTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                        >  1:173581/1‑142 (MQ=255)
aaCTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAagag                                                                                                                                        <  2:173581/142‑1 (MQ=255)
                                    atcaACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTc                                                                                            >  2:371317/1‑150 (MQ=255)
                                                      cAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCGGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgc                                                                          <  2:72300/150‑1 (MQ=255)
                                                                 gCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                            <  1:335481/72‑1 (MQ=25)
                                                                 gCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                            >  2:335481/1‑72 (MQ=25)
                                                                       gccgTACCTCCACCTGGTCCGCTCCCCAAGCTCTCCTGTCCTGGCCTCTTCGGGGGCTGGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatc                                                         <  1:6815/147‑1 (MQ=255)
                                                                           tACCTGCACCTGATCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGGGGGCTTGCCCGCGAAGAGGCCAGACGGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATcc                                                     <  1:119293/150‑1 (MQ=255)
                                                                                         cggCTGACCATGCTCTCCTGTCCTGGCCTCTTCGGGGGTTTGCCGGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                       <  1:236290/148‑1 (MQ=255)
                                                                                                                    ctctTCGCGGGCAAGCCCGCGAAGAGGCTAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc            >  2:119816/1‑150 (MQ=255)
                                                                                                                    ctctTCGCGGGCAAGCCCGCGAAGAGGCTAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCAc            >  2:119817/1‑150 (MQ=255)
                                                                                                                     tctTCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACg           >  1:10875/1‑150 (MQ=255)
                                                                                                                       ttCGCGGGCAAGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTc         >  2:295660/1‑150 (MQ=255)
                                                                                                                              gCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTcc  >  1:70030/1‑150 (MQ=255)
                                                                                                                                |                                                                                                                                                   
AACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGAGGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTCTCC  >  NC_002947/284459‑284734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: