Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,102,890	T→C	R8G (AGA→GGA)	insH1 ←	IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,102,890	0	T	C	100.0%	109.8 / NA	31	R8G (AGA→GGA)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/13);  total (18/13)

TCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTT  >  NC_000913/2102811‑2103025                                                                                |                                                                                                                                        tctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATc                                                                              >  1:994730/1‑139 (MQ=11)  ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATcc                                                                             >  2:44984/1‑139 (MQ=14)  ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATcc                                                                             <  2:1114941/139‑1 (MQ=14)    tttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCt                                                                           >  1:574453/1‑139 (MQ=18)               aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGa                                                                <  2:97968/139‑1 (MQ=255)                     cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAGACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAAt                                                          >  2:803769/1‑139 (MQ=255)                     cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATg                                                                 >  1:365551/1‑132 (MQ=255)                     cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATg                                                                 <  2:365551/132‑1 (MQ=255)                                aTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGaaa                                               >  2:134930/1‑139 (MQ=255)                                       gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGac                                                               <  1:334127/116‑1 (MQ=255)                                       gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGac                                                               >  2:334127/1‑116 (MQ=255)                                        tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGa                                       >  2:440299/1‑139 (MQ=255)                                         cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGaa                                      >  1:536423/1‑139 (MQ=255)                                         cGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGaa                                      >  2:446849/1‑139 (MQ=255)                                          ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGa                                                                >  1:401978/1‑112 (MQ=255)                                          ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGa                                                                <  2:401978/112‑1 (MQ=255)                                                  tAAGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTa                             <  2:774468/139‑1 (MQ=255)                                                  tAAGTTGATGACTCATGATGAACCCTGTTCCATGGCGCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTa                             <  1:705685/139‑1 (MQ=255)                                                    aGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTa                             <  1:206068/137‑1 (MQ=255)                                                    aGTTGATGACTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTa                             >  2:206068/1‑137 (MQ=255)                                                       tgatgaCTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaa                        <  2:574453/139‑1 (MQ=255)                                                           gaCTCATGATGAACCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGaataat                    >  1:489967/1‑139 (MQ=255)                                                                      aaCCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAt         <  1:134930/139‑1 (MQ=255)                                                                       aCCCTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAtt        >  1:274508/1‑139 (MQ=255)                                                                        cccTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCt                               >  2:1105458/1‑115 (MQ=255)                                                                        cccTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCt                               <  1:1105458/115‑1 (MQ=255)                                                                        cccTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAttt       >  1:672122/1‑139 (MQ=255)                                                                          cTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAAc                                                                        <  2:740941/72‑1 (MQ=255)                                                                          cTGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAAc                                                                        >  1:740941/1‑72 (MQ=255)                                                                           tGTTCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGc    >  1:894490/1‑139 (MQ=255)                                                                             ttCCATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCtt  <  1:1030201/139‑1 (MQ=255)                                                                                |                                                                                                                                        TCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTT  >  NC_000913/2102811‑2103025

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