Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,186,190	A→C	I2L (ATT→CTT)	insC1 →	IS2 repressor TnpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,186,190	0	A	C	100.0%	69.0 / NA	22	I2L (ATT→CTT)	insC1	IS2 repressor TnpA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/11);  total (11/11)

ATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCG  >  NC_000913/3186042‑3186244                                                                                                                                                     |                                                       attaaCCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGAtgctt                                                      <  1:893313/151‑1 (MQ=255)        acaATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCt                                               >  2:667671/1‑151 (MQ=255)          aaTAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCtt                                              >  2:850124/1‑150 (MQ=255)                  aaTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCgg                                      >  2:581008/1‑150 (MQ=255)                   aTTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCgg                                      >  1:744751/1‑149 (MQ=255)                     tGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGgaga                                   <  2:477272/150‑1 (MQ=255)                     tGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCAAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGgaga                                   >  1:797833/1‑150 (MQ=255)                      gATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGaa                                  <  2:470378/150‑1 (MQ=255)                         tAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACg                               <  1:443230/150‑1 (MQ=255)                            aaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGa                           <  1:136052/151‑1 (MQ=255)                              atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACg                         <  2:39939/151‑1 (MQ=255)                               tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACgg                        >  2:778209/1‑151 (MQ=255)                                atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                       >  1:298520/1‑151 (MQ=255)                                   tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGt                     >  1:745427/1‑150 (MQ=255)                                    attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTAc                   >  2:167280/1‑151 (MQ=255)                                       gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                 <  2:48800/150‑1 (MQ=255)                                         tCAGTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCGCTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACcaca               <  2:224753/150‑1 (MQ=255)                                                ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaa        <  2:91345/150‑1 (MQ=18)                                                 cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAg      <  1:584270/151‑1 (MQ=18)                                                  gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaa        >  2:420841/1‑148 (MQ=17)                                                  gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGa     <  2:292314/151‑1 (MQ=17)                                                     attaGACTGGCCCCCTGAATCTCCAGACAAGCAGTATCACTTAAATAGGTGATAGGCTTAATGCTAGTTTTTAGACTAGGCATTGGAGAACAGAGGCTTGATGTCGTAGTGCCGGAGAAACGCAGACGGCGTACCACACAGGaaa        >  1:445589/1‑145 (MQ=11)                                                     attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCg  <  2:146899/151‑1 (MQ=11)                                                                                                                                                     |                                                       ATTAACCACAATAAGCGAATTGATTAAAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCG  >  NC_000913/3186042‑3186244

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑