Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,356,932	G→A	T60I (ACC→ATC)	sapA ←	putative periplasmic binding protein SapA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,932	0	G	A	100.0%	53.7 / NA	18	T60I (ACC→ATC)	sapA	putative periplasmic binding protein SapA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/10);  total (8/10)

TGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCAT  >  NC_000913/1356690‑1357169                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              tGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTTGAAGGTAAACACCTCATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCTAGTTCTTCCCAGCTTTCGGCAAGTTCCTGCATCAGGCGATATGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGGAATGATATCGAc                                                                                                                                                                                                                                        <  2:219490/250‑1 (MQ=255)                ggTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGTCAATTAACCCACTGCTc                                                                                                                                                                                                                        <  2:57133/251‑1 (MQ=255)                        cGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTgg                                                                                                                                                                                                                 >  2:238208/1‑250 (MQ=255)                                           aGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACt                                                                                                                                                                                             >  1:38222/1‑251 (MQ=255)                                                                                            gTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCg                                                                                                                                            >  2:25776/1‑251 (MQ=255)                                                                                                  tttttGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGga                                                                                                                                       <  2:55685/250‑1 (MQ=255)                                                                                                              aaCATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCCGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGACGCGGATATCAGCATGCGGGGGAGATTCCGGCgcgg                                                                                                                          >  2:166788/1‑251 (MQ=255)                                                                                                                   cgcgACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGAtt                                                                                                                     <  1:34499/251‑1 (MQ=255)                                                                                                                             gTGGAAGCGATAGGAAGCGCCGTTGTCGAGTCCTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGCTTGCCTGACCAc                                                                                                           >  1:241453/1‑251 (MQ=255)                                                                                                                                                        gagtgctgcccAGCTTTCGGCAAGTTCCGGCATCAGTCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAGCTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGGTAAAGGTGTTGACTTGGCCGCTTACGCAATAGGCAAAACCGCGGTCGCGGAGATCATCATGCTGGGGATATTCAGGGGCGGCGATTGCCTTTCCACGCTCAAGGCCAGCAATCACGCAAagag                                                                                <  2:28919/243‑1 (MQ=255)                                                                                                                                                           tACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGCCCACTCACAAGTCCAGCAATCACCAAAAGAGACg                                                                             >  2:240059/1‑251 (MQ=255)                                                                                                                                                                  cAGCTTTCGGAAAGTGCCGGCATCAGGCGATAGGTATAGGGAGCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAAGTAGCCTACGGCTCGCTGGGGATGGGTTAAAGGTGTTGACTTGCTCGTTGACGAAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACGCACAAGACCAGCAATCACCAAAAGAGACGagaagac                                                                      <  2:269308/251‑7 (MQ=255)                                                                                                                                                                                       aTCAGGCGATAGGTATAGGGATCGTCATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGa                                                 <  2:193893/251‑1 (MQ=255)                                                                                                                                                                                                             cGACATCCAGCAGTCGATCAGAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGATTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATAGCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCt                            <  2:267759/250‑1 (MQ=255)                                                                                                                                                                                                                cATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACt                        <  1:197388/251‑1 (MQ=255)                                                                                                                                                                                                                    cagcagTCGATCATAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAat                     >  1:161975/1‑250 (MQ=255)                                                                                                                                                                                                                                aTAAAACTGGGCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACAtt        >  2:92378/1‑251 (MQ=255)                                                                                                                                                                                                                                        gggCGGCAAGGATATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCAt  <  1:240059/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                              TGACGTTATGCCACGGGTTGTTGCGGTCAAAAATTCGCTGGAAGGTAAACACCACATCGTCGGCATTCATTTTACGAGTGGGAGTAAACCAGTCGGTTTTTTGAAACGGAACATCGCGACGCAGGTGGAAGCGATAGGTCGCGCCGTTGTCGAGTACTTCCCAGCTTTCGGCAAGTTCCGGCATCAGGCGATAGGTATAGGGATCGACATCCAGCAGTCGATCATAAAACTGGGCGGCAAGGGTATCGACAATTAACCCACTGCTCGCTTTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCAT  >  NC_000913/1356690‑1357169

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