Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,294,403	Δ1 bp	intergenic (+259/+23)	tdk → / ← insZ	thymidine/deoxyuridine kinase/putative transposase InsZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,294,403	0	C	.	100.0%	99.7 / NA	24	intergenic (+259/+23)	tdk/insZ	thymidine/deoxyuridine kinase/putative transposase InsZ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (8/16);  total (8/16)

CACCGACGTTGACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATACGGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGC  >  NC_000913/1294276‑1294536                                                                                                                                |                                                                                                                                      cACCGACGTTGACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAggg                                                                                                                          <  2:304934/140‑1 (MQ=255)           gACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGtt                                                                                                                <  2:365802/140‑1 (MQ=255)           gACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGtt                                                                                                                <  1:371227/140‑1 (MQ=255)           gACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGtt                                                                                                                <  2:249692/140‑1 (MQ=255)                   aCCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAag                                                                                                        >  2:339191/1‑140 (MQ=255)                        tAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAgaggag                                                                                                   <  1:453490/140‑1 (MQ=255)                             tACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCa                                                                                              <  2:66503/140‑1 (MQ=255)                             tACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCa                                                                                              <  2:365241/140‑1 (MQ=255)                                aaTATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTa                                                                                           <  1:339191/140‑1 (MQ=255)                                                                        gatgaAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCTGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGttttt                                                                   <  2:401910/124‑1 (MQ=255)                                                                        gatgaAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCTGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGttttt                                                                   >  1:401910/1‑124 (MQ=255)                                                                            gaGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCAAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCgg                                               <  2:67110/139‑1 (MQ=255)                                                                               gCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCt                                            <  2:28704/140‑1 (MQ=255)                                                                               gCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCt                                            >  1:178523/1‑140 (MQ=255)                                                                                 tAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTg                                          <  1:461051/140‑1 (MQ=255)                                                                                         tCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTg                                          >  2:314672/1‑132 (MQ=255)                                                                                         tCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTg                                          <  1:314672/132‑1 (MQ=255)                                                                                         tCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCgg                                  >  2:389886/1‑140 (MQ=255)                                                                                              ttATGGAAAAGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCg                             >  1:175630/1‑140 (MQ=255)                                                                                                     aaaGTTGCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGtt                      >  2:124349/1‑140 (MQ=255)                                                                                                           gCTTTGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAac                >  2:335586/1‑140 (MQ=255)                                                                                                              ttGGGTAAACAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAAcacc             <  1:389886/140‑1 (MQ=255)                                                                                                                     aaCAAAAAATA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCa      <  1:124349/140‑1 (MQ=255)                                                                                                                         aaaaaTA‑GGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGc  <  2:175630/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                      CACCGACGTTGACTTGTGACCTGTAAAGTACAATATCCCTGTGTTTAGGCGTTATACATCGTCGCAAATATGATGAAGGCTAATGCTGTCGGTTTATGGAAAAGTTGCTTTGGGTAAACAAAAAATACGGCCCCAGAAGGGCAATGCCGTTCACTTAAGAGGAGCGGCACTATGTTTCACAGGATAACGGGTTTTTGATATCTTAACCGACCTCGGCCTTGATGGTCGGGGGCGTTTTGTTATGAACACCACTTCCAGAGC  >  NC_000913/1294276‑1294536

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