Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,504,878	T→A	intergenic (+13/+27)	insQ → / ← ydcO	putative insertion element transposase InsQ/putative transport protein YdcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,504,878	0	T	A	100.0%	59.1 / NA	18	intergenic (+13/+27)	insQ/ydcO	putative insertion element transposase InsQ/putative transport protein YdcO
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (9/9);  total (9/9)

CGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTT  >  NC_000913/1504754‑1505010                                                                                                                             |                                                                                                                                     cgcTCGCAACATTTTAGCGGCGGGGCACCCCGGTCTTTCCTGTGGGGGGAGGGGGGAGTCAGGCCGCCCCCTGAAGCAGGGACCCACCCAAAAGATTTAGGCGGCAGCCCGGACGGAGCAGGGGACCACGGCCCTCAggg                                                                                                                       >  2:156470/1‑140 (MQ=255)   acagCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCg                                                                                                                     <  1:213590/137‑1 (MQ=255)        aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTggag                                                                                                                <  1:263588/140‑1 (MQ=255)        aaCATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTggag                                                                                                                <  1:59077/140‑1 (MQ=255)           aTTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGAt                                                                                                             <  2:102670/140‑1 (MQ=255)             tttAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGt                                                                                                           <  1:50394/140‑1 (MQ=255)                                 ttCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGc                                                                                       >  1:185801/1‑140 (MQ=255)                                           tgGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTc                                                                             >  2:474902/1‑140 (MQ=255)                                               gagaTGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAaca                                                                         <  2:185801/140‑1 (MQ=255)                                                          gTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACaa                                                              >  1:189147/1‑140 (MQ=255)                                                            cAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAAc                                                            >  1:163939/1‑140 (MQ=255)                                                             aGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACg                                                           >  1:389420/1‑140 (MQ=255)                                                               gccgccCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGcc                                                         >  2:84537/1‑140 (MQ=255)                                                                                       ccGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGc                                 <  2:389420/140‑1 (MQ=255)                                                                                        cGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCa                                >  1:142928/1‑140 (MQ=255)                                                                                                     gCGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGa                   >  1:202614/1‑140 (MQ=255)                                                                                                      cGACAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGAc                  <  2:142928/140‑1 (MQ=255)                                                                                                        aCAGCCTGAACGTAGCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCa                >  2:251251/1‑140 (MQ=255)                                                                                                                      gCAGGGAACCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACtt  <  1:474902/140‑1 (MQ=255)                                                                                                                             |                                                                                                                                     CGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGGATGTCAATATCGGTTTCTGTCAGCGATTAAATTCAACACCACGTAACAAACGCCTCCGGCAATTAATCCCCAAAACGCAGAACCAATCCCGACCAGCGTCAATCCACTT  >  NC_000913/1504754‑1505010

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