Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730602,4660GT43.8% 7.7 / 10.0 16intergenic (+56/‑101)hchA/araBpossible hydrolase/ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (3/6);  new base T (6/1);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 6.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AATTATTAACAGGTGACAGTCCTTTAGCTTCAAATGAGTTAGGAAAATTAGCAGTTAATGAAATGTTAAATGCAATACAAAATAAATAATTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGATTATGGAACTGCTTCAGGCCGTGT  >  CP000730/602322‑602616
                                                                                                                                                 |                                                                                                                                                      
aaTTATTAACAGGTGACAGTCCTTTAGCTTCAAATGAGTAAGGAAAATTAGCAGTTAATGAAATGTTAAATGCAATACAAAATAAATAATTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCaaaa                                                                                                                                                  <  2:419868/151‑1 (MQ=255)
                                aaTGAGTTAGGAAAATTAGCAGTTAATGAAATGTTAAATGCAATACAAAATAAATAATTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCtatta                                                                                                                  >  2:162457/1‑151 (MQ=255)
                                                               agTTAAATGAAATACAAAATAAAAAAATAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAAc                                                                                   <  1:335532/150‑1 (MQ=255)
                                                                                      taatTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATg                                                            <  1:33308/151‑1 (MQ=255)
                                                                                        atTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATgga                                                          <  2:554128/151‑1 (MQ=255)
                                                                                         ttAAATATTAATTAGAGGATCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGCATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATggag                                                         <  1:782585/151‑1 (MQ=255)
                                                                                           aaaTATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTTTCAAATTTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATggagga                                                       >  2:824107/1‑151 (MQ=255)
                                                                                               attaattaGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTTTTTTTTTTTAAGTGATAATTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAACACTAATAACAGCAAATAGAGATATGGAGGaaaaa                                                   >  2:609658/1‑149 (MQ=255)
                                                                                                   attaGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACatgat                                               <  1:162457/151‑1 (MQ=255)
                                                                                                      agagGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTTGTAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGGTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAATATGATGTc                                            >  2:920009/1‑151 (MQ=255)
                                                                                                            gCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTTTTTTTTTTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAg                                      >  2:973656/1‑151 (MQ=255)
                                                                                                            gCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTTTTATTATTTAATTGATACTTGTAAAATAGAACCTTTTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGTTGTCTTATAg                                      >  2:469820/1‑151 (MQ=255)
                                                                                                                     gTAAATGTATGAGGGCTCTTTTTTTTTTTTTTTTTTTTAGTGATACTTGTACAATAGAACCTATTATGTGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGTAACTTGTTGTCTTATAGCATTTGAAt                             >  2:900738/1‑151 (MQ=255)
                                                                                                                          tGTATGAGGGCTCTTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGAt                         >  1:826134/1‑151 (MQ=255)
                                                                                                                              tGAGGGATCTTTTTTTTTTGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGATTATgg                    <  1:184359/151‑1 (MQ=255)
                                                                                                                                        cttttttttGGCAAAATTTAAGTGATACTTGTAAAATAAAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGATTATGGAACAGCTTc           >  1:14055/2‑151 (MQ=255)
                                                                                                                                                 ggCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGATTATGGAACTGCTTCAGGCCgtgt  >  1:99164/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
AATTATTAACAGGTGACAGTCCTTTAGCTTCAAATGAGTTAGGAAAATTAGCAGTTAATGAAATGTTAAATGCAATACAAAATAAATAATTAAATATTAATTAGAGGAGCCTCATATGTAAATGTATGAGGGCTCTTTTTTTTTGGCAAAATTTAAGTGATACTTGTAAAATAGAACCTATTATGAGTATGATTTAAGAAAACGCTTGCAAAACTAATAACCGCAACTAGCGATATGGAGGAAACATGATGTCTTATAGCATTGGAATTGATTATGGAACTGCTTCAGGCCGTGT  >  CP000730/602322‑602616

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: