Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,876	C→T	34.1%	intergenic (+107/+1)	rlmH → / ← tpn	hypothetical protein/IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,876	0	C	T	34.1%	69.7 / 27.3	43	intergenic (+107/+1)	rlmH/tpn	hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base T (8/7);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAC  >  CP000730/35737‑36015                                                                                                                                            |                                                                                                                                            gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCt                                                                                                                                            >  2:102897/1‑141 (MQ=255)     ttCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAAc                                                                                                                                        <  1:135827/141‑1 (MQ=255)        aGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTg                                                                                                                                     <  1:247013/141‑1 (MQ=255)                 tttttGGAGAAAAATGATTACTTGAACAAGAGCTAAAATAAAACTGTTATCTAATATACATAGTAGGTTTAATAAAAAAAAAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCatga                                                                                                                            <  2:68041/141‑1 (MQ=18)                    ttGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGc                                                                                                                         <  1:165555/141‑1 (MQ=255)                     tgacataaaaTGATTACGTGAAGAAGAGCTAAAATAAAATTGTGAACAAATAAAAAGAGATTGTGCTGTTGAAAAGTAAAGAAAGATAAATAACCACTTACGTACCATGTCAGTGTGTGCTTAGATTGCAAGTAAGAAGCt                                                                                                                        <  2:136756/135‑1 (MQ=255)                      ggCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTa                                                                                                                       >  2:259498/1‑141 (MQ=255)                       gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAAGTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                                                                                      <  2:84671/141‑1 (MQ=255)                               tGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTgg                                                                                                              <  2:106618/141‑1 (MQ=255)                                     cTTGAACAAGAGCTAGAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATTCTAATTTCGTGTCATGGc                                                                                                        >  2:26323/1‑141 (MQ=255)                                      ttGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCg                                                                                                       <  2:165007/141‑1 (MQ=255)                                       tGAACAAGAGATAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGa                                                                                                      <  2:242453/141‑1 (MQ=255)                                             agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                                                >  2:98718/1‑141 (MQ=255)                                              gagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCg                                                                                               >  1:32546/1‑141 (MQ=255)                                                    aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                         <  1:284414/141‑1 (MQ=255)                                                      acataagcgtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATATCTAATCACTAATATATCTTGTCAGTGTTTGTTTAACTTGCTAGCATTATCCTAATTTCGTGGCATGGCGACATTCCGTAGTCCTGaa                                                                                     >  2:34872/15‑141 (MQ=2)                                                      acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                     >  2:99850/15‑141 (MQ=11)                                                           atTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag                                                                                  >  1:208316/1‑141 (MQ=255)                                                              gtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                             >  2:230409/7‑141 (MQ=11)                                                                gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg                                                                             <  1:55424/141‑1 (MQ=255)                                                                attgtaataaAAAT‑GATGATTCTTTTGCAAAGT‑AAAAAATATAGCTAGTCACTAATATATCTAGTCAGTGTTCGTTTACCTTGCTAGCATGATGATAATTTCGTGCCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg                                                                           <  2:185032/137‑1 (MQ=11)                                                                 ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                          >  1:289171/4‑141 (MQ=11)                                                                 ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg                                                                          <  1:15779/138‑1 (MQ=11)                                                                   ttaataaAAAT‑GATGGTACTGTTTCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGATGAGACCTGCGGtt                                                                        <  2:232123/140‑1 (MQ=11)                                                                    taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                         <  1:324997/141‑1 (MQ=40)                                                                     aataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                      <  2:57244/141‑1 (MQ=11)                                                                      ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                     >  1:132431/1‑141 (MQ=11)                                                                      ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCtt                                                                     >  1:77376/1‑141 (MQ=11)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  >  1:228216/1‑141 (MQ=34)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  >  2:27155/1‑141 (MQ=34)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  <  2:327748/141‑1 (MQ=34)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  >  1:21386/1‑141 (MQ=34)                                                                           aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                  <  2:306748/141‑1 (MQ=34)                                                                              aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTtata                                                               >  1:260123/1‑141 (MQ=32)                                                                                                                aatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTa                             <  1:79717/138‑1 (MQ=0)                                                                                                                                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt         >  1:304316/1‑141 (MQ=9)                                                                                                                                    aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt         <  1:210161/141‑1 (MQ=14)                                                                                                                                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGATAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt        <  2:227820/141‑1 (MQ=11)                                                                                                                                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt        >  2:57682/1‑141 (MQ=17)                                                                                                                                       tttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTg      <  2:71292/140‑1 (MQ=255)                                                                                                                                       gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg      <  2:208316/141‑1 (MQ=32)                                                                                                                                         tCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGAt    >  2:207943/1‑141 (MQ=11)                                                                                                                                           gCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc  >  1:174953/1‑141 (MQ=17)                                                                                                                                            |                                                                                                                                            GCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAC  >  CP000730/35737‑36015

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