Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,762	A→G	C254R (TGC→CGC)	thiP ←	thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,762	0	A	G	100.0%	122.7 / NA	33	C254R (TGC→CGC)	thiP	thiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/17);  total (16/17)

GTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAA  >  NC_000913/73647‑73891                                                                                                                    |                                                                                                                                  gtgtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCa                                                                                                            >  2:139923/1‑139 (MQ=255) gtgtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCa                                                                                                            >  1:721167/1‑139 (MQ=255)  tgtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAg                                                                                                           >  2:766170/1‑139 (MQ=255)   gtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGc                                                                                                          >  2:246497/1‑139 (MQ=255)   gtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGc                                                                                                          >  2:146962/1‑139 (MQ=255)    tCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGcaca                                                                                                                              <  1:292101/118‑1 (MQ=255)    tCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGcaca                                                                                                                              >  2:292101/1‑118 (MQ=255)                cTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCggg                                                                                             <  1:136017/139‑1 (MQ=255)                       gACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATc                                                                                      <  1:146962/139‑1 (MQ=255)                             cGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTc                                                                                >  2:666602/1‑139 (MQ=255)                                      cgcgCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAg                                                                       <  1:766170/139‑1 (MQ=255)                                        cgcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAgcg                                                                     <  1:487401/139‑1 (MQ=255)                                        cgcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAgcg                                                                     >  2:350249/1‑139 (MQ=255)                                         gcCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAgcgc                                                                    <  1:182801/139‑1 (MQ=255)                                               cTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                              <  1:686504/139‑1 (MQ=255)                                                        gcGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAg                                                     <  1:139923/139‑1 (MQ=255)                                                             gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                              >  2:70675/1‑125 (MQ=255)                                                             gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATa                                                              <  1:70675/125‑1 (MQ=255)                                                             gtgCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGa                                                <  1:666602/139‑1 (MQ=255)                                                                       cAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGcaca                                                                                                                              <  2:616456/51‑1 (MQ=255)                                                                       cAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGcaca                                                                                                                              >  1:616456/1‑51 (MQ=255)                                                                            gCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCg                                 <  1:350249/139‑1 (MQ=255)                                                                                      aTCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGGCAGCATTGCCGCGCGGGGAGGATCGTAGTCGTAACTCAGCTCCTGATTGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccccc                       <  2:403171/139‑1 (MQ=255)                                                                                      aTCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGc                                                        <  2:711014/106‑1 (MQ=255)                                                                                      aTCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGc                                                        >  1:711014/1‑106 (MQ=255)                                                                                      aTCGCTGACTCAACAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccccc                       <  1:80753/139‑1 (MQ=255)                                                                                                 aaCAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccc                         <  1:762951/126‑1 (MQ=255)                                                                                                 aaCAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccc                         >  2:762951/1‑126 (MQ=255)                                                                                                   cAGCACCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACg          >  2:239217/1‑139 (MQ=255)                                                                                                       aCCAGCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTgg      >  1:711660/1‑139 (MQ=255)                                                                                                           gCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccccc                       >  1:409512/1‑118 (MQ=255)                                                                                                           gCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGccccc                       <  2:409512/118‑1 (MQ=255)                                                                                                           gCCCGAGGCGGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGaa  >  2:587497/1‑139 (MQ=255)                                                                                                                    |                                                                                                                                  GTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAA  >  NC_000913/73647‑73891

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