Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,956	G→T	A194S (GCG→TCG)	yhhI →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,956	0	G	T	97.7%	93.9 / NA	45	A194S (GCG→TCG)	yhhI	putative transposase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (20/23);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATT  >  NC_000913/3624840‑3625093                                                                                                                     |                                                                                                                                          tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACAGGGGCGGCta                                                                                                                 >  1:311056/1‑143 (MQ=17)       aaattattgaaaTGCAGGATATTAAAGGAATAAACATAACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCAGCtaaataaa                                                                                                           <  1:207448/132‑1 (MQ=11)              tAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAACAGGAACCCAGGGGCGGCTAAATAAAGCCTTTg                                                                                                    >  2:253762/1‑143 (MQ=17)                   tgccgggtaggatgggAAAAATCATGACAAAGCAGGCGATCCGGAGCCGGATAGATATTTTAGAGAAGAGACAAAAACAGGGAGGTGATTATTTATTCTCGGGAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTGGAGGaa                                                                                               <  1:380535/129‑1 (MQ=2)                     cTGGATATTAAAGGAAAAATCATCACAACTGATGCGCTGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGaaaa                                                                                             >  1:355706/1‑143 (MQ=17)                        gATAATAAAGGAAAAGTCATCAAAACAGATGAGAAGGATTGCGAAAAAGTTATTGCAGAGAAAAAACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAAtt                                                                                          <  1:338317/143‑1 (MQ=11)                             tAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAACAGGAACCCAGGGGCGGCCAAATCAAGCCTTTGAGGAAAAATTTCCGc                                                                                     >  2:930478/1‑143 (MQ=255)                             tAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGc                                                                                     >  1:427775/1‑143 (MQ=17)                                 ggAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaa                                                                                 >  2:499912/1‑143 (MQ=17)                                  gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                <  2:243903/143‑1 (MQ=17)                                  gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                >  2:15401/1‑143 (MQ=17)                                  gAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACCGGGCCACCCAAAGAAAGCCTTTGAGGAAAAATTTCCGCTGaaa                                                                                >  1:208692/1‑143 (MQ=11)                                    aaaaTCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGa                                                                              <  2:668545/143‑1 (MQ=17)                                     acgccagcccgaCTGATGCGACGGGTGGCCAGAAGGATATTGTAGAGAAGATACAAAAACAGGGAGGGTGTTAATACTTCTGGGTCACGCGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGaa                                                                             <  2:425939/132‑1 (MQ=2)                                       atcctaacaacaGATGCGATAGGATGAGAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAtt                                                                           <  1:90573/137‑1 (MQ=11)                                       aacagaacaggagATGCGATGGGGTGCCAGAAAGGTATTGCAGGGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAAtt                                                                           <  1:954273/131‑1 (MQ=11)                                         catcaCACCTGATGCGATGGGATGTCAGAATAATATTGCAGAAAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGTAACCCAGTGGCGGCTAAATAAAGTCTTTGAGGAAAAATTTCCGCTGAAAGAATTaa                                                                         <  1:141974/143‑1 (MQ=11)                                           tcaCAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaat                                                                       >  2:241262/1‑143 (MQ=17)                                              caACTGAAGCGAGGGGGTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAaataat                                                                    <  2:335290/143‑1 (MQ=255)                                                aCTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATcc                                                                  >  2:707668/1‑143 (MQ=17)                                                   gATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCaga                                                               >  1:421944/1‑143 (MQ=18)                                                    aagCGAGGGGATGCCAGAAAGATATGGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCagag                                                              <  1:730156/141‑1 (MQ=255)                                                       cGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGTGGAAAAATTTCCGCTTAAAGAATTAAATAATCCAGAGCAt                                                           <  2:249868/143‑1 (MQ=255)                                                       cGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCAt                                                           >  2:88632/1‑143 (MQ=25)                                                           gggTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACa                                                       <  2:100798/143‑1 (MQ=25)                                                                  cAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAGATAATCCAGAGCATGACAGTTACGc                                                >  2:371754/1‑143 (MQ=25)                                                                   aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCa                                               >  1:470830/1‑143 (MQ=25)                                                                       aGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCa                                           <  2:355706/143‑1 (MQ=38)                                                                          tatTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAACTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTg                                        <  2:811876/143‑1 (MQ=255)                                                                            ttGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACACAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGaa                                      <  2:848917/143‑1 (MQ=255)                                                                             tGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCTCGGTAATATGATCCTAGGTGCGTCTATATTAAGGCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGaaa                                     <  2:831094/143‑1 (MQ=11)                                                                               ttgagaAGATACAACAAGAGGGAGGTGATGCTTTATGCGAGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCTCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAag                                   <  2:673606/141‑1 (MQ=11)                                                                                    aagaTACAATAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCCGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              >  1:643948/1‑143 (MQ=255)                                                                                    aagaTACAAAAACAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              <  2:606434/143‑1 (MQ=38)                                                                                    aagaTACAAAAACAGGGAGGCGACTATTTATTCTCGGTAAAAGGAACGCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAACTTTCCGCTGACAGAATGAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCa                              <  2:482204/143‑1 (MQ=11)                                                                                           aaaatatgtgAGTTGATTAGTTATTGTCGGTAAAAGTAACCCAGGGGCGGCTAAATAAAGCCTTTTAGGAAAAATTTCCGCTGAAAGAATTAAATAATCTAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCaga                       <  1:598618/134‑1 (MQ=11)                                                                                           aaaaaCAGGGAGGTGATTATTTATTATCGGCAAAAGGGACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCaga                       >  2:697488/1‑143 (MQ=255)                                                                                            aaaaCAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCagag                      <  1:460151/143‑1 (MQ=38)                                                                                             aaaCAGGGAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAgaga                     >  1:628699/1‑143 (MQ=255)                                                                                                    gAGGTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTGACGGCAGAGAAGaaacc              >  1:573811/1‑141 (MQ=255)                                                                                                       gTGATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGGTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGt           <  2:531879/143‑1 (MQ=255)                                                                                                        tGATAATTTATTTACGACAAAAGAAACCCAAGGGCAGCTAAATAAAGTCTTTGAGGAAAAATTTCCGATGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTc          <  1:462598/143‑1 (MQ=11)                                                                                                         gATTATTTATTCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGGAGTATATAATCCAGAGCATGACAGGTACGCAATCAGTGAAAAGAGTCACGGCAGAGGAGAAATCGGtgt         >  2:38168/1‑141 (MQ=11)                                                                                                               ttattCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCatat   >  1:286921/1‑143 (MQ=255)                                                                                                                tattCTCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATAtt  <  2:39561/143‑1 (MQ=38)                                                                                                                     |                                                                                                                                          TCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAACCCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGAGCATGACAGTTACGCAATCAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATT  >  NC_000913/3624840‑3625093

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