Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	316,453	G→A	*289* (TAG→TAA)	insF1 →	IS3 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	316,453	0	G	A	100.0%	126.4 / NA	41	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/23);  total (18/23)

CGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGA  >  NC_000913/316346‑316584                                                                                                            |                                                                                                                                    cGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCAt                                                                                                      >  2:1528852/1‑139 (MQ=11)      aaCGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGt                                                                                                 >  2:94899/1‑139 (MQ=21)             tttAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCggg                                                                                          >  1:54153/1‑139 (MQ=255)                aaTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                       <  2:1650739/139‑1 (MQ=255)                 aTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGt                                                                                      <  1:1012556/139‑1 (MQ=255)                       tCGAATGTTATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                <  1:1399821/139‑1 (MQ=255)                       tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                       <  2:329475/132‑1 (MQ=255)                       tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCg                                                                                       >  1:329475/1‑132 (MQ=255)                       tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                >  2:715524/1‑139 (MQ=255)                       tCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                <  1:1528852/139‑1 (MQ=255)                                 ttACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGc                                                                      >  1:288481/1‑139 (MQ=255)                                      aTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTc                                                                 <  1:94899/139‑1 (MQ=255)                                       tCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                >  2:1101427/1‑123 (MQ=255)                                       tCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                <  1:1101427/123‑1 (MQ=255)                                        cGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg                                                               <  2:556072/139‑1 (MQ=255)                                         ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg                                                               <  2:778303/138‑1 (MQ=255)                                         ggTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAg                                                               >  1:778303/1‑138 (MQ=255)                                          gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             <  2:54153/139‑1 (MQ=255)                                          gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             <  1:715524/139‑1 (MQ=255)                                          gTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             <  2:1326905/139‑1 (MQ=255)                                           tGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                <  1:1485432/119‑1 (MQ=255)                                           tGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                >  2:1485432/1‑119 (MQ=255)                                                              gtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTcc                                                                                               >  2:1308887/1‑85 (MQ=255)                                                              gtgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTcc                                                                                               <  1:1308887/85‑1 (MQ=255)                                                               tgGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCa                                        <  2:1447685/139‑1 (MQ=255)                                                                      cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             <  1:1612227/111‑1 (MQ=255)                                                                      cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             >  2:1612227/1‑111 (MQ=255)                                                                      cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGttt                                 <  1:771037/139‑1 (MQ=255)                                                                      cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGttt                                 >  2:771037/1‑139 (MQ=255)                                                                       tCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTc                                >  2:91202/1‑139 (MQ=255)                                                                                     ttagAAAACAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAAt                  >  1:785673/4‑139 (MQ=255)                                                                                          aaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                                                                                 >  2:1117264/1‑71 (MQ=255)                                                                                          aaaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTcc                                                                                 <  1:1117264/71‑1 (MQ=255)                                                                                           aaCAAGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCg            <  1:452515/139‑1 (MQ=255)                                                                                               aGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                >  1:1466385/1‑67 (MQ=255)                                                                                               aGAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCa                                                                                <  2:1466385/67‑1 (MQ=255)                                                                                                gAACCTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGAtttt       >  2:588414/1‑139 (MQ=255)                                                                                                   ccTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             >  2:749774/1‑82 (MQ=255)                                                                                                   ccTCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCa                                                             <  1:749774/82‑1 (MQ=255)                                                                                                     tCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGa  <  1:91202/139‑1 (MQ=255)                                                                                                     tCGCTTAAGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGa  <  2:288481/139‑1 (MQ=255)                                                                                                            |                                                                                                                                    CGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATAAAGTCCTCGGGTCGTTGTCCACTCTTTTGGCATGTCAGCATATTCACATACCGTACCAGCAAAAGTTTCAACCGTGAGCAAATCGACCGATTTTGTTGA  >  NC_000913/316346‑316584

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