Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,677,049	T→C	A285A (GCA→GCG)	pntA ←	pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,677,049	0	T	C	100.0%	83.3 / NA	24	A285A (GCA→GCG)	pntA	pyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/12);  total (12/12)

TAAAGCTGTGAGGATTGCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCTGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGACAATGATATCGACCTCTTT  >  NC_000913/1676924‑1677165                                                                                                                              |                                                                                                                     tAAAGCTGTGAGGATTGCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGgccgcc                                                                                                                   <  2:534537/129‑1 (MQ=255) tAAAGCTGTGAGGATTGCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGgccgcc                                                                                                                   >  1:534537/1‑129 (MQ=255)              aTTGCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACt                                                                                                  >  1:583023/1‑133 (MQ=255)                tGCGTCGGCACACGGTCCGGAAGACCGGCATACCCAATCACGTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGc                                                                                                <  2:162666/133‑1 (MQ=255)                 gCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGcc                                                                                               <  1:101281/133‑1 (MQ=255)                 gCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGcc                                                                                               <  2:583023/133‑1 (MQ=255)                           aCGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATg                                                                                     >  1:447516/1‑133 (MQ=255)                                         cGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTc                                                                       >  1:87811/1‑133 (MQ=255)                                                     tCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACg                                                                    >  2:7946/1‑124 (MQ=255)                                                     tCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACg                                                                    <  1:7946/124‑1 (MQ=255)                                                                 cATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTTATGGAGTCGACCATCTGGCGGGTAATTAGCTTCGGCGCtcg                                               >  2:300641/1‑131 (MQ=255)                                                                  aTTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGt                                              <  2:437316/133‑1 (MQ=255)                                                                       ccGTAGTGAATATTTGACCCGGCAAGGTATATTCACAGTTGACGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGcc                                         <  2:269647/133‑1 (MQ=255)                                                                             tGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGc                                          >  1:155573/1‑126 (MQ=255)                                                                             tGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGc                                          <  2:155573/126‑1 (MQ=255)                                                                             tGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAAt                                   >  1:588317/1‑133 (MQ=255)                                                                                 gATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAgc                               >  2:449997/1‑133 (MQ=255)                                                                                    ttCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGcgcg                            <  2:447516/133‑1 (MQ=255)                                                                                    ttCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGcgcg                            <  1:449997/133‑1 (MQ=255)                                                                                       aCCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCggtg                         >  1:97004/1‑133 (MQ=255)                                                                                       aCCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCggtg                         <  1:528719/133‑1 (MQ=255)                                                                                            gCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGAc                    >  1:577932/1‑133 (MQ=255)                                                                                            gCACGGTGAATTCACTGTTGCCGCCGTTTTAGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCggtggctgt                    <  1:300641/133‑5 (MQ=255)                                                                                                       tCACAGTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGACAATGATATCGa         <  2:588317/133‑1 (MQ=255)                                                                                                            gTTGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGCCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGACAATGATATCGACctct    >  1:46823/1‑133 (MQ=255)                                                                                                              tGCCGCCGTTTTGGGCCGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGACAATGATATCGACCTCttt  >  2:4486/1‑133 (MQ=255)                                                                                                                              |                                                                                                                     TAAAGCTGTGAGGATTGCGTCGGCAGACGGCCCGGAAGATCGGTATAACCAATCACTTTGACACCATTTTCCGTAGTGAAGATTTCACCCGGCACGGTGTATTCACAGTTGCCGCCGTTTTGGGCTGCCAGGTCGACAATCACACTGCCCGCCTTCATGGAGTCAACCATTTCACGGGTAATTAGCTTCGGCGCTGGTTTGCCTGGAATAAGCGCGGTGGTGACAATGATATCGACCTCTTT  >  NC_000913/1676924‑1677165

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