Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,503,797	A→G	intergenic (‑366/‑261)	insG ← / → yjhB	IS4 transposase/putative MFS transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,503,797	0	A	G	100.0%	50.7 / NA	15	intergenic (‑366/‑261)	insG/yjhB	IS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (3/12);  total (3/12)

AAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTA  >  NC_000913/4503657‑4503930                                                                                                                                             |                                                                                                                                      aaGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACt                                                                                                                             <  1:430780/151‑1 (MQ=255)  aGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACtg                                                                                                                            >  2:92583/1‑151 (MQ=255)   gAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACtgt                                                                                                                           <  2:672814/151‑1 (MQ=255)           gCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAgtg                                                                                                                    <  2:339982/150‑1 (MQ=255)                             aaCGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACa                                                                                                 <  2:361551/151‑1 (MQ=255)                                    gCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGtttt                                                                                          <  1:37853/151‑1 (MQ=255)                                       aaaGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGttttt                                                                                         <  1:753420/149‑1 (MQ=255)                                                          cTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGt                                                                     >  1:574204/1‑150 (MQ=255)                                                                            gggAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTAt                                                  <  1:8997/151‑1 (MQ=255)                                                                               aGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCtt                                               <  1:92584/151‑1 (MQ=255)                                                                                   ttATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCa                                            <  2:72652/150‑1 (MQ=255)                                                                                       tCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATaga                                        <  1:503232/150‑1 (MQ=255)                                                                                                  ggACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGtt                            >  2:342139/1‑151 (MQ=255)                                                                                                                  cacaATGAATAAGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATggagga            <  2:420454/151‑1 (MQ=255)                                                                                                                            aaGTCAATTTTTGAAGCGGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTa  <  2:547318/151‑1 (MQ=255)                                                                                                                                             |                                                                                                                                      AAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTCAATTTTTGAAGCAGGGTTGAACTGTGCCAGTGGTGGAATATCTGGCGAACATTGTTTTTTGTTGGTGACCCAAAACTGTAGGTGGGTATCAGTTTTATCTTTCATAGAGTGAAATATGTTAAGAAGAAATGGAGGAAAAGATTGTA  >  NC_000913/4503657‑4503930

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