Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,707,883	G→A	intergenic (‑178/+733)	dppA ← / ← proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,707,883	0	G	A	100.0%	72.6 / NA	22	intergenic (‑178/+733)	dppA/proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/8);  total (14/8)

TCAGCATCCCTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATG  >  NC_000913/3707672‑3708116                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          tCAGCATCCCTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATAc                                                                                                                                                                                                                 <  2:82573/238‑1 (MQ=255) tCAGCATCCCTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATAc                                                                                                                                                                                                                 >  1:82573/1‑238 (MQ=255)  cAGCATCCCTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGAt                                                                                                                                                                                                    <  2:122537/250‑1 (MQ=255)          cTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCa                                                                                                                                                                                            >  1:109952/1‑250 (MQ=255)                              cgcatagttCTGCTACAAGTGTGCTGTTTGTAGTTTTAACCCTTTGAAGTAGGTTCTCGCTGCCTGACCTTTTTTGTTTTTTGCCAGGTCGGGTCAACGTTATGAGGTGGTGATGCCGTACTAATTAACATCAGTGTAACAGTTCGGCTCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGAGAGAATCcacaaaca                                                                                                                                                                        >  2:174742/8‑250 (MQ=255)                                         gCTCCAATTGTGATTTTTGTTTTTTTAACTCTTGGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTGTTGCCCGGTCGGGACAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGACAACCGGGATTTACACCAACGGTGAGAACCCACAAACAAAGATTaaaaaa                                                                                                                                                            <  2:1396/251‑1 (MQ=255)                                                    gATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACa                                                                                                                                                  <  1:47781/250‑1 (MQ=255)                                                           gttgttTTAACCCTTTGCAGTGAGTTGTCGCTTCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGAGGGGATGCAGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTCCAATCTATCAACAGAATGTGACAACGTCAATACAGCCAACCGGGCTTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTg                                                                                                                                          >  2:29708/1‑251 (MQ=255)                                                                 ttAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTa                                                                                                                                              >  1:214747/1‑241 (MQ=255)                                                                 ttAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTa                                                                                                                                              <  2:214747/241‑1 (MQ=255)                                                                  tAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTgcagc                                                                                                                                      >  1:38727/1‑248 (MQ=255)                                                                                                aCCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGCGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAACTGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCCACGGTGAGAATCAACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCAAACATGACCAGa                                                                                                      >  2:125993/1‑250 (MQ=255)                                                                                                         gttttttgCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACATCCAACCGGGATTTACATCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCt                                                                                            <  1:247555/251‑1 (MQ=255)                                                                                                         gttttttgCCAGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTATTTATTAACCTCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGACAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCt                                                                                            >  2:243218/1‑251 (MQ=255)                                                                                                                            cAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAAGCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACAttt                                                                          <  2:109952/250‑1 (MQ=255)                                                                                                                                  taagaggaggggcggCCGTACTAATTAACATCAGTGTCCCAGAGCGGATCCTACATAAAATGCCACTCGTGACCTACAATCTATCAACAGACTGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTAAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGAtt                                                                   >  2:60732/15‑251 (MQ=255)                                                                                                                                                           tAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGCTTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaaa                                           >  2:19248/1‑250 (MQ=255)                                                                                                                                                                   gtgtAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAgg                                  >  2:49430/1‑251 (MQ=255)                                                                                                                                                                     gtAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGa                                 >  2:79234/1‑250 (MQ=255)                                                                                                                                                                       aaCAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGaaaa                              <  2:214546/251‑1 (MQ=255)                                                                                                                                                                                   cctccATAAAATGCCCCTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAACAAATCTGAGGAAAAGGTGGATATCTg                  >  1:140875/1‑251 (MQ=255)                                                                                                                                                                                                   cTCGTGACCTACAATCTATCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATg  >  1:142679/1‑251 (MQ=255)                                                                                                                                                                                                                    |                                                                                                                                                                                                                                          TCAGCATCCCTGACTTTTTCAAGGAAATACGCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATG  >  NC_000913/3707672‑3708116

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