Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,532,746	T→A	W311R (TGG→AGG)	ydcC →	H repeat‑associated putative transposase YdcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,532,746	0	T	A	100.0%	74.0 / NA	28	W311R (TGG→AGG)	ydcC	H repeat‑associated putative transposase YdcC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/12);  total (16/12)

CTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAA  >  NC_000913/1532610‑1532784                                                                                                                                         |                                       cTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                     <  2:1193503/140‑1 (MQ=255)  tttCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCg                                    >  2:1273005/1‑140 (MQ=255)       gTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTgg                               <  2:1243852/140‑1 (MQ=255)         ccATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             >  1:131109/1‑140 (MQ=255)          cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             >  1:1078223/1‑139 (MQ=255)          cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGAc                             <  2:1078223/139‑1 (MQ=255)          cATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                            <  1:1226186/140‑1 (MQ=255)           ataataGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                           >  2:414221/1‑140 (MQ=255)               taGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                       <  2:1339058/140‑1 (MQ=21)               taGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTa                       <  2:838171/140‑1 (MQ=16)                  cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                    >  2:973314/1‑140 (MQ=255)                  cAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                    >  1:908107/1‑140 (MQ=25)                      aCAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaatg                <  1:454189/140‑1 (MQ=16)                       cAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAatgaatga               >  2:761426/1‑140 (MQ=255)                       cAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAatgaatga               >  1:913450/1‑140 (MQ=255)                           agaaagaGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:1273005/140‑1 (MQ=16)                           agaaagaGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:1060593/140‑1 (MQ=16)                               agagCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACt       >  2:496773/1‑140 (MQ=25)                                 agCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACtac     >  1:947400/1‑138 (MQ=16)                                  gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAGGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:581659/1‑140 (MQ=255)                                  gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACtaca    <  1:414221/140‑4 (MQ=255)                                  gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:1067293/1‑140 (MQ=16)                                  gCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa    >  1:706981/1‑140 (MQ=16)                                    tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           <  1:993173/131‑1 (MQ=2)                                    tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac           >  2:993173/1‑131 (MQ=17)                                    tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACta      >  2:752146/1‑135 (MQ=12)                                    tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACta      <  1:752146/136‑2 (MQ=2)                                    tCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaa  >  1:245375/1‑140 (MQ=12)                                                                                                                                         |                                       CTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAA  >  NC_000913/1532610‑1532784

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