Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,624,8450AT22.6% 46.0 / 13.8 31P156P (CCA→CCT) yhhIputative transposase
Reads supporting (aligned to +/- strand):  ref base A (17/7);  new base T (7/0);  total (24/7)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGG  >  NC_000913/3624706‑3624980
                                                                                                                                            |                                                                                                                                       
gAAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCa                                                                                                                                         <  1:531221/140‑1 (MQ=34)
gAAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCa                                                                                                                                         <  1:531285/140‑1 (MQ=34)
   aaaCGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGaa                                                                                                                                      <  2:690146/140‑1 (MQ=34)
   aaaCGCTCCGGCACTCTTATGACAAGAGTCGCCGAAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGaa                                                                                                                                      <  2:753352/140‑1 (MQ=21)
        ctacgtcACTCTTAAGACAATAGACGCCGAAGGGGAGCGATTCATGTCATTATTGAGTTCTCAA‑AAAACAAACTAATGGTCATATGACAGATAAAGACGGATCAGAAATCTAATGAGATTACAGCTATCCAAGAActtca                                                                                                                                 >  2:18626/7‑139 (MQ=2)
          ccGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTa                                                                                                                               <  1:613450/140‑1 (MQ=34)
           cGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGTGATTACAGCTATCCCAGAAATTCTTaa                                                                                                                              >  2:1463066/1‑140 (MQ=21)
           cGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTaa                                                                                                                              >  2:665998/1‑140 (MQ=34)
                  ctTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTg                                                                                                                       <  1:1143845/140‑1 (MQ=32)
                     aTGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGat                                                                                                                    >  2:1521869/1‑140 (MQ=32)
                     aTGACAAGAGTCGCCGCAGGGGAGCAATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGat                                                                                                                    >  2:919133/1‑140 (MQ=17)
                                                                                                                              aTTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAc            >  2:1035031/1‑140 (MQ=11)
                                                                                                                                 aCAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTACAAGGAACCCa         >  2:637698/1‑140 (MQ=17)
                                                                                                                                 aCAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCa         >  2:637673/1‑140 (MQ=32)
                                                                                                                                   aGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgg       <  2:1263630/140‑1 (MQ=32)
                                                                                                                                   aGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgg       <  2:1069787/140‑1 (MQ=32)
                                                                                                                                     cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg     <  2:881270/140‑1 (MQ=32)
                                                                                                                                       aTCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg   >  2:888769/1‑140 (MQ=17)
                                                                                                                                       aTCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATATATTCGCGGTAAAAGGAACCCAGGggcg   >  1:40642/1‑140 (MQ=17)
                                                                                                                                       aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg   >  1:802468/1‑140 (MQ=32)
                                                                                                                                        tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:954535/1‑140 (MQ=17)
                                                                                                                                        tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:195572/1‑140 (MQ=17)
                                                                                                                                        tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:1171128/1‑140 (MQ=17)
                                                                                                                                        tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTAAAAATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:795209/1‑140 (MQ=255)
                                                                                                                                        tCCCTGAACATCTTAACATGCTGGATATTCAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGACCCCAGGggcgg  >  1:392310/1‑140 (MQ=255)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAATATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:223174/1‑140 (MQ=17)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:1102099/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:1273058/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:1416622/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  2:621892/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:825848/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg  >  1:635533/1‑140 (MQ=32)
                                                                                                                                        tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGAGGTAAAAGGAACACAGGggcgg  >  2:39621/1‑140 (MQ=17)
                                                                                                                                            |                                                                                                                                       
GAAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGT‑CTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGG  >  NC_000913/3624706‑3624980

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: ‑