Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8580TA69.7% ‑6.0 / 28.8 33intergenic (+89/+19)rlmH/tpnhypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base A (9/14);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATA  >  CP000730/35731‑35956
                                                                                                                               |                                                                                                  
gTAGAAGCTATTCAGTATAAAGTTATTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGAAAAAAAAAAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCtt                                                                                <  2:351717/148‑1 (MQ=25)
      gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTg                                                                          <  2:90029/148‑1 (MQ=255)
                                          cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGGTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                      <  1:279458/148‑1 (MQ=255)
                                          cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                      <  1:158655/148‑1 (MQ=255)
                                          cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                      <  1:333145/148‑1 (MQ=255)
                                          cTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                      <  2:42741/148‑1 (MQ=255)
                                                            caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                  <  1:60801/135‑1 (MQ=11)
                                                            caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                  <  1:331039/135‑1 (MQ=11)
                                                            caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                  <  2:33664/135‑1 (MQ=11)
                                                             aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  1:190793/136‑1 (MQ=11)
                                                             aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  1:71905/136‑1 (MQ=11)
                                                             aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  1:206277/136‑1 (MQ=11)
                                                             aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc                 <  1:190792/136‑1 (MQ=11)
                                                               aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAGAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATATGAAGAGACCTGCgg               >  1:320595/11‑148 (MQ=11)
                                                               aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg               >  2:9587/11‑148 (MQ=11)
                                                               aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg               >  1:224004/11‑148 (MQ=11)
                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  1:132189/139‑1 (MQ=11)
                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  2:7246/139‑1 (MQ=11)
                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  1:85028/139‑1 (MQ=11)
                                                                agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACAAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt              <  1:24603/139‑1 (MQ=11)
                                                                  agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc            >  2:184058/8‑148 (MQ=11)
                                                                     attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt         >  1:359239/5‑148 (MQ=11)
                                                                     attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttt         >  1:369677/5‑148 (MQ=11)
                                                                      atttaacaaAAATAGAAGGTTCCGTTGCAAAGTAAAAAAAAATAGCTAACCACTAATTTATCATGTCAGTGGTCGCCTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAACTCCGTAGACGTGAAGAGACCTGCGGCTCtt          >  2:229311/4‑148 (MQ=11)
                                                                      aTCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt          >  2:205607/1‑148 (MQ=255)
                                                                       tttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt       <  2:104065/146‑1 (MQ=11)
                                                                         taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt       <  2:115761/148‑1 (MQ=40)
                                                                         taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt       <  2:213786/148‑1 (MQ=40)
                                                                           ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatat   >  1:18264/1‑148 (MQ=11)
                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata  >  2:27321/1‑148 (MQ=11)
                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata  <  1:2206/148‑1 (MQ=11)
                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata  >  2:164503/1‑148 (MQ=11)
                                                                            taaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtatata  <  2:54256/148‑1 (MQ=11)
                                                                                                                               |                                                                                                  
GTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATA  >  CP000730/35731‑35956

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