Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	16,704	T→C	intergenic (+147/+47)	insL1 → / ← hokC	IS186/IS421 family transposase/protein HokC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	16,704	0	T	C	100.0%	55.1 / NA	23	intergenic (+147/+47)	insL1/hokC	IS186/IS421 family transposase/protein HokC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/4);  total (19/4)

TTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGA  >  NC_000913/16486‑16920                                                                                                                                                                                                                             |                                                                                                                                                                                                                         ttAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCACGAAGTGCCGGATCCGAAAAGCAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAACCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACt                                                                                                                                                                                          >  1:198420/1‑251 (MQ=9)          gacATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAACAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGt                                                                                                                                                                                 >  2:66739/1‑251 (MQ=21)                          cGCTGAATTCCCCCCCCAAAAGGGCGGGATCCGAAAAGAAGAACTAACTCGTGGGGGAGAAAAAAAAAAAGGGTCATCTGGACCTTACAGGTGGCCAT‑GCGTGGG‑CAAGATCCCCTGCCAGCCTAAAAAACGCAATTAAATAACGCGGGGCATCGTCTTAAGGAGCCACCGGGGCGCCGCATTTTACTGTGGGCTCAACCCTTAAGTAAGCGTTTATGGGACCACTCCCCGCCGTTGCTCTTACTCGGAtt                                                                                                                                                                 <  1:66739/251‑1 (MQ=255)                           gCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGCACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTc                                                                                                                                                                >  2:122380/1‑251 (MQ=255)                              ggATTTCCCCCC‑AGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCCTCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTa                                                                                                                                                             >  1:196955/1‑250 (MQ=255)                                                 gCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAACATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCAC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGCTGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGTCGTGACAACAGCAAcc                                                                                                                                          >  2:346976/1‑251 (MQ=255)                                                       tCCGAAAAGCAGAACTAAATCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTc                                                                                                                                    >  1:379762/1‑251 (MQ=255)                                                                aagaaCTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGtt                                                                                                                           >  2:131655/1‑251 (MQ=255)                                                                  agactacatttgtgtgGAAAATAACAAAAAGGGTCAGCGGGTCATTACGGGTGGC‑ATGTCGGGGCACAGTATCCT‑CGACGTCCTACAAAACGCAAGGAAAAAACGGGGGGCGTCTTCTTAACGAGGCACGGAGGCGTCGCTTTCTACAGATGGCTCACCCCTAAATTGAGCGTGTAGGGGATCACCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGTCCAGTTCg                                                                                                                         <  2:164445/239‑1 (MQ=255)                                                                                   agaATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCccgccg                                                                                                                                                                                  >  1:246942/1‑177 (MQ=37)                                                                                   aaaaaaaCAAAAATGTTCATCGGGAGCTTACAGGGGGCCAT‑TCTGGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATGAAAGAGCGCTAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTGCAGAGGGCTCAACCTGTAAGGTGGGGTTTAGGGGATCACTCCccgccg                                                                                                                                                                                  <  2:246942/172‑1 (MQ=2)                                                                                     aaTAACAAAAATGGTAATCTGGAGCTTACAGGTGGCCCT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAACCAGCAACCTCAGTCTGGCCAGTTCGGATGTGAACCTCACAGAgg                                                                                                      >  2:80562/1‑251 (MQ=255)                                                                                            aaaaTGGTCATCTGGAGATTACAGGTGGCCAT‑TCGTGCGACCGTATCAC‑CGACAGCCTACACAACGCAATTGAAGAACGCGCGGCATCGTATTACCGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTCTGGGATCACTCCACGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCACCCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTtc                                                                                               >  1:435815/1‑251 (MQ=21)                                                                                                   tCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTAc                                                                                        >  1:259737/1‑251 (MQ=255)                                                                                                            gcctacgggggggccttgggggggACATTACCCC‑GGACAGCCCACAAACGGCAATGAAAGACCGGGGGGCATCTTCTTACCGAGGCACGGAGGCTTCCCATTCTCCAGAGGGCTCACCCCTTAAGTTGGCGTTTAGGGGATCCCCCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGcgccg                                                                                <  2:263052/230‑1 (MQ=255)                                                                                                                  aGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCTAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTAc                                                                         >  1:318023/1‑251 (MQ=255)                                                                                                                       gCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCtt                                                                                                                                                                          >  1:436761/1‑149 (MQ=255)                                                                                                                           t‑tCGTGGGACAGTATCCC‑TGAAAGCCTCCAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCAAAGCGTCGCATTCTTCAGATGGCTCAACCATTAAGTTAGCGCTTATGGGATCACTCCCCGCAGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGACACTCACAGAGGACTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTgg                                                                 >  2:81510/1‑249 (MQ=25)                                                                                                                                 gggACAGTATCCC‑TGACAGACTACAAAACGCCATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGGCTTTTCTCGTTACCAGCGCCGCCACTAcggcggggcaacagg                                                           >  2:294476/1‑242 (MQ=255)                                                                                                                                                                      aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGCTGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATACAACTGACGATCAGGGCGGCAAGCATCACCTTATcctgct                       >  2:51506/1‑251 (MQ=37)                                                                                                                                                                                    tCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGCGATACAGATGAAGATCCGGGCGACAATAATTGCCTTAagctggctccttgttcccttt         >  2:417219/1‑235 (MQ=38)                                                                                                                                                                                      gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTtctc       >  1:102398/1‑251 (MQ=255)                                                                                                                                                                                           aaCGAGGCACCGAGGCGCCGCATTCTTCAGATGGCTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTAGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTAGTTACCAGCGCCGCCACTACGGCGGGGATAAAGATGAAGATAAGGGGGAAAATCATCGCCTTCTGCTGCTTCAtttgtcttttcttcttgg  >  1:207947/1‑235 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                         TTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCAT‑TCGTGGGACAGTATCCC‑TGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGA  >  NC_000913/16486‑16920

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