Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,049,923	C→A	L31I (CTC→ATC)	insA‑4 →	IS1 protein InsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,049,923	0	C	A	100.0%	23.2 / NA	14	L31I (CTC→ATC)	insA‑4	IS1 protein InsA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (7/7);  total (7/7)

CTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTTGGATGCCGGGCAACAGCCC  >  NC_000913/1049792‑1050043                                                                                                                                    |                                                                                                                         cTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGctct                                                                                                                  <  2:93463/140‑1 (MQ=17)   tACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCGCTGCCGGACATCAGCGCTATATCTGCTCTCa                                                                                                                <  2:243100/140‑1 (MQ=17)    aCTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCAc                                                                                                               <  1:514701/140‑1 (MQ=17)          ttAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGt                                                                                                         >  1:66001/1‑140 (MQ=17)                  tgatgGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATg                                                                                                 >  1:225999/1‑140 (MQ=17)                       gtgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAAc                                                                                            >  1:437633/1‑140 (MQ=17)                        tgtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACt                                                                                           >  2:620740/1‑140 (MQ=17)                         gtTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTg                                                                                          <  1:298152/140‑1 (MQ=17)                          tttttGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGc                                                                                         >  1:59071/1‑140 (MQ=18)                                  gTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACt                                                                                 >  1:229150/1‑140 (MQ=17)                                  gTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACt                                                                                 >  2:485536/1‑140 (MQ=17)                                                                                       ggcgtggtgCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGaa                            <  2:403322/137‑1 (MQ=1)                                                                                                                cTGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTTGGATGCCGGGCAACAGcc   <  1:253392/140‑1 (MQ=14)                                                                                                                 tGCCGGACATCAGCGCTATATCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTTGGATGCCGGGCAACAGccc  <  2:589336/140‑1 (MQ=11)                                                                                                                                    |                                                                                                                         CTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTTCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACTGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTTGGATGCCGGGCAACAGCCC  >  NC_000913/1049792‑1050043

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