Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,507,463	G→A	intergenic (+12/+3)	insO → / ← insI1	pseudogene, IS911 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related/IS30 transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,507,463	0	G	A	100.0%	50.7 / NA	18	intergenic (+12/+3)	insO/insI1	pseudogene, IS911 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related/IS30 transposase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/11);  total (7/11)

CTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACAT  >  NC_000913/4507332‑4507595                                                                                                                                    |                                                                                                                                     cTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                               <  1:405320/139‑1 (MQ=255) cTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAAt                                                                                                                               <  1:648564/139‑1 (MQ=255)                     cAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCtt                                                                                                           <  2:577540/139‑1 (MQ=255)                               aaaGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACACCTTTCAATTAtct                                                                                                 >  1:818773/1‑139 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                           <  2:833353/86‑1 (MQ=255)                                                          aTCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGt                                                                                                                           >  1:833353/1‑86 (MQ=255)                                                                              gAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGc                                                  <  1:753929/139‑1 (MQ=255)                                                                                             gATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGa                                   >  2:819337/1‑139 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                            >  1:202532/1‑129 (MQ=255)                                                                                              aTTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGagca                                            <  2:202532/129‑1 (MQ=255)                                                                                                       gagaATGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCAt                         <  1:627081/139‑1 (MQ=255)                                                                                                            tGAAATATTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGa                    <  1:1220/139‑1 (MQ=255)                                                                                                                  atTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTc                                                                             <  2:575621/76‑1 (MQ=255)                                                                                                                  atTAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTc                                                                             >  1:575621/1‑76 (MQ=255)                                                                                                                  atTAAAAAGGCTACTGTAAAGTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTatat              <  2:96833/139‑1 (MQ=255)                                                                                                                    tAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                       >  1:768335/1‑112 (MQ=255)                                                                                                                    tAAAAAGGCTACTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAAc                                       <  2:768335/112‑1 (MQ=255)                                                                                                                              aCTGTAAATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTATCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATCTTTTGCACGACAt  >  1:224797/1‑139 (MQ=11)                                                                                                                                    |                                                                                                                                     CTGCGTGATGAGCGTCAGGGCAAAACACCAAAAGCCTCTCCGATAACACCAGAACAAATCGAAATACGTAAGCTGAGGAAAAAGCTACAACGGATTGAAATGGAGAATGAAATATTAAAAAGGCTACTGTAGATTCAATCTGTCAATGCAACACCCCTTTCAATTATCTCTTTCGGTGTTTTGAACTTCAGTGTCTTTCTCGGTCTGTTGTTTAGCTGAGCAGCAACCAGATCTAGTTCATGTTGAGTATATTGGGCAAGACAT  >  NC_000913/4507332‑4507595

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