Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,922,549	T→G	T146P (ACG→CCG)	yqcA ←	putative flavodoxin YqcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,922,549	0	T	G	100.0%	85.3 / NA	25	T146P (ACG→CCG)	yqcA	putative flavodoxin YqcA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/14);  total (11/14)

GCAAACGCCAGTAAAAATCCGCGCTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAAC  >  NC_000913/2922429‑2922667                                                                                                                         |                                                                                                                       gCAAACGCCAGTAAAAATCCGCGCTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTc                                                                                                            >  2:23901/1‑133 (MQ=255)   aaaCGCCAGTAAAAACCCGCGCTAAGAATCAGATTTTACTGCCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGa                                                                                                          <  2:173299/133‑1 (MQ=255)     aCGCCAGTAAAAATCCGCGCTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGGGGGCCCCACGGTTCGAcc                                                                                                        >  1:190300/1‑133 (MQ=255)                     gcgcTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCg                                                                                        >  1:131380/1‑133 (MQ=255)                      cgcTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGt                                                                                       >  1:307642/1‑133 (MQ=255)                      cgcTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGt                                                                                       >  1:512081/1‑133 (MQ=255)                        cTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTCACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGttt                                                                                     <  2:126880/133‑1 (MQ=255)                              aGCAGTTTCTACTGGAATCTGCCTGAAAAGCTTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCC‑CGGATTTGATTCCGTttccggttc                                                                              <  1:178122/133‑1 (MQ=255)                               gCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGTACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTttccggttc                                                                              <  1:17633/133‑1 (MQ=255)                                cAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTttccggttcc                                                                             <  2:426300/133‑1 (MQ=255)                                  gTTTTTACCGGCATTCGCCTGACAAGATTCTATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTttccggttccgg                                                                           <  1:160632/133‑1 (MQ=255)                                     tttACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGtt                                                                        <  2:411941/133‑1 (MQ=255)                                               tttGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCg                                                              <  2:470050/133‑1 (MQ=255)                                                ttGCCTGAGAAGATTTGATTCAACACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTggtgt                                                             <  1:208217/133‑4 (MQ=255)                                                  gCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGtca                                                           <  2:307642/133‑1 (MQ=255)                                                                   tcAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGcc                                          >  1:142963/1‑133 (MQ=255)                                                                            cTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCGCGGATTTGATTCCGTTTCCAGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTg                                 >  1:81151/1‑133 (MQ=255)                                                                                aaCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAgcg                             <  1:23901/133‑1 (MQ=255)                                                                                 aCGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAgcgc                            <  1:64703/133‑1 (MQ=255)                                                                                             gctcgctGAACATATCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGca                >  2:74720/1‑133 (MQ=255)                                                                                                cgctGAACATTTCAGCACAACAGCTGGCACCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGcaaca             <  2:183061/133‑1 (MQ=255)                                                                                                          atCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGaa   >  2:130519/1‑133 (MQ=255)                                                                                                          atCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGaa   >  1:130269/1‑133 (MQ=255)                                                                                                           tCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAAc  <  2:142963/133‑1 (MQ=255)                                                                                                           tCAGGACAACAGCGGGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAAc  >  2:346614/1‑133 (MQ=255)                                                                                                                         |                                                                                                                       GCAAACGCCAGTAAAAATCCGCGCTATGAAGCAGTTTTTACTGGCATTTGCCTGAAAAGATTCGATTCAGCACCGCTAAAACGACATTTACCGCTCGCTGAACATATCAGGACAACAGCGTGCCCCACTGTTCGACCCACGGATTTGATTCCGTTTCCGGTTCCGGGTTCTCACTGGCGTCAATCAACAGCATTTCGCCAACCCGCTGAGCGCTCTGTTCCTGCAACAAGGCATCGAAC  >  NC_000913/2922429‑2922667

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