Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,624,8450AT20.4% 78.1 / 23.5 44P156P (CCA→CCT) yhhIputative transposase
Reads supporting (aligned to +/- strand):  ref base A (21/14);  new base T (9/0);  total (30/14)
Fisher's exact test for biased strand distribution p-value = 3.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGC  >  NC_000913/3624707‑3624981
                                                                                                                                          |                                                                                                                                        
aaaaaCGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAg                                                                                                                                         <  1:1273467/140‑1 (MQ=34)
       cTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAActtct                                                                                                                                  >  2:632414/1‑140 (MQ=34)
         ccGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTa                                                                                                                                >  1:986263/1‑140 (MQ=34)
         ccGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTa                                                                                                                                <  1:1435445/140‑1 (MQ=34)
         ccGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGAGATTCATGTCATTAGTGAGTTCTCAACAATGAACAGTCTGGTCATCGGACAGATCAAGAAGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTa                                                                                                                                <  2:1251405/140‑1 (MQ=255)
              aCTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATg                                                                                                                           <  2:50421/140‑1 (MQ=33)
                tgttATGACAAGAGTAGCAGCAGGGGAGCGATTCATGTCATTAGTGAGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAATACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTATCATGCt                                                                                                                         <  2:1654406/138‑1 (MQ=11)
                tctTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCt                                                                                                                         <  1:354897/140‑1 (MQ=32)
                actTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCTTTCGTGCGTTCTCAACACTGCACAGACTGGTCATCGGACAGTTCAAGACGGATGAGCACTCTAATGAGCTTACAGCTCTCCCAGAACTTCTTAACATGCt                                                                                                                         <  2:1654402/139‑1 (MQ=11)
                       aCAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATACCAGAACTTCTTAACATGCTGGATAtt                                                                                                                  >  1:1129638/1‑140 (MQ=14)
                                                                                                                            atcaCAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAc             >  1:1393908/4‑140 (MQ=11)
                                                                                                                            atcaCAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAc             >  2:622438/4‑140 (MQ=11)
                                                                                                                            atcaCAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAc             >  2:1556731/4‑140 (MQ=11)
                                                                                                                                  gCTATCCCAGAACTTCTTAACATGCTTGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAggg       <  2:310231/140‑1 (MQ=17)
                                                                                                                                   cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg      <  2:572675/140‑1 (MQ=32)
                                                                                                                                   cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg      <  2:260511/140‑1 (MQ=32)
                                                                                                                                   cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg      <  2:883332/140‑1 (MQ=32)
                                                                                                                                   cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg      <  1:1105001/140‑1 (MQ=32)
                                                                                                                                   cTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAgggg      <  1:14249/140‑1 (MQ=32)
                                                                                                                                     aTCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  1:182045/1‑140 (MQ=17)
                                                                                                                                     aTCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  2:877277/1‑140 (MQ=17)
                                                                                                                                     aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTTCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  1:80992/1‑140 (MQ=17)
                                                                                                                                     aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  1:78016/1‑140 (MQ=32)
                                                                                                                                     aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  1:1745145/1‑140 (MQ=32)
                                                                                                                                     aTCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcg    >  1:1573009/1‑140 (MQ=32)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAATATCATCACAACTGTTGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAATTGGAACCCAGGggcgg   >  1:1121033/1‑140 (MQ=11)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAATATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1671827/1‑140 (MQ=17)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:290240/1‑140 (MQ=17)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1023595/1‑140 (MQ=17)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:560120/1‑140 (MQ=17)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1612367/1‑140 (MQ=17)
                                                                                                                                      tCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGAGGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTACAAGGAACCCCGGggcgg   >  1:290178/1‑140 (MQ=255)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGTAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1340717/1‑140 (MQ=17)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  2:1621875/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1678373/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1680486/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  2:1123472/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1511937/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  2:417467/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:1511693/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:957159/1‑140 (MQ=32)
                                                                                                                                      tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcgg   >  1:915583/1‑140 (MQ=32)
                                                                                                                                       cccAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGggcggc  <  1:58840/140‑1 (MQ=32)
                                                                                                                                       cccAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTAATTCGCGGTAAAAGGAACCCAGGggcggc  >  1:1774515/1‑140 (MQ=17)
                                                                                                                                          |                                                                                                                                        
AAAAACGCTCCGGCACTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGACAGATCAAGACGGATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGC  >  NC_000913/3624707‑3624981

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