Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,317,069	C→T	G50G (GGG→GGA)	phnM ←	RPnTP hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,317,069	0	C	T	100.0%	115.5 / NA	36	G50G (GGG→GGA)	phnM	RPnTP hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/15);  total (21/15)

TGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACG  >  NC_000913/4316951‑4317204                                                                                                                       |                                                                                                                                        tgATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCt                                                                                                                    >  2:502344/1‑140 (MQ=255)                      gcgcGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCGGGCCGCAACCAGCCGCGTTCGACGTCCATCGCCTCCGGc                                                                                               >  2:282891/1‑140 (MQ=255)                        gcgTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAg                                                                                             >  1:799103/1‑140 (MQ=255)                              tGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGct                                                                                       >  1:759958/1‑140 (MQ=255)                                        tCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCg                                                                             >  1:23186/1‑140 (MQ=255)                                              aGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAg                                                                       >  1:979764/1‑140 (MQ=255)                                              aGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAg                                                                       >  2:959105/1‑140 (MQ=255)                                                   gCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGcgcg                                                                  >  2:631403/1‑140 (MQ=255)                                                      ggCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGAt                                                               >  2:741905/1‑140 (MQ=255)                                                          aGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCg                                                           >  2:82430/1‑140 (MQ=255)                                                           gTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCgt                                                          <  1:830206/140‑2 (MQ=255)                                                            tCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGcc                                                         >  1:819523/1‑140 (MQ=255)                                                                 tttCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCt                                                    <  2:184953/140‑1 (MQ=255)                                                                      ggCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCAcc                                               <  1:819159/140‑1 (MQ=255)                                                                      ggCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGATTTAGGCAAAGGGGCGGCTTTCGCCGTTCTGCAcc                                               >  2:373027/1‑140 (MQ=255)                                                                            gggTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGc                                         >  1:618779/1‑140 (MQ=255)                                                                               tGAAGAATTTATCCATATTATCGGTCTGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGaa                                      <  1:401485/140‑1 (MQ=255)                                                                               tGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGaa                                      <  2:979764/140‑1 (MQ=255)                                                                                 aagaaTTTATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTTGCTTTCGGAAAAGGCGCGGATTTCGTCGTTATGCACCTCCAGCGAAcc                                    >  2:348425/1‑140 (MQ=255)                                                                                    aaTTTATCCAGATTAACGGTATCCACCTCCATCAGTCACGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCt                                 <  1:282891/140‑1 (MQ=255)                                                                                       ttATCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTac                              >  1:679435/1‑140 (MQ=255)                                                                                         aTCCAGATTATCGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTacca                            <  2:111598/140‑1 (MQ=255)                                                                                                    cGGTATGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTcc                 <  1:631403/140‑1 (MQ=255)                                                                                                        aTGCAGCTCAATCAGTCCCGGCGGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa             <  1:348425/140‑1 (MQ=255)                                                                                                        aTGCAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCa             <  2:618779/140‑1 (MQ=255)                                                                                                           cAGCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCCTTCGGCAAAGGCGCGGATTTCGGCGTTCTGCTCCTCCCGAGAACCGCTCACCAAATCGTTTTCCAGCACCa          >  2:439505/1‑140 (MQ=255)                                                                                                             gCTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGc        >  1:700162/1‑140 (MQ=255)                                                                                                             gCTCAATCAGTCCCGGCAGCAACCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGc        >  2:99713/1‑140 (MQ=255)                                                                                                              cTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCt       >  2:1008092/1‑140 (MQ=255)                                                                                                              cTCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCt       >  1:1008761/1‑140 (MQ=255)                                                                                                               tCAATCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCtt      <  2:759958/140‑1 (MQ=255)                                                                                                                   tCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg  <  1:82430/140‑1 (MQ=255)                                                                                                                   tCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg  <  2:840809/140‑1 (MQ=255)                                                                                                                   tCAGTCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCAGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg  >  1:742698/1‑140 (MQ=255)                                                                                                                   tCAGTCCCGGCAGCAACCAGCCGCCTTCGCCGTCCATTGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg  <  1:373027/140‑1 (MQ=255)                                                                                                                   tCAGTCCCGGCAGCAACCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACg  <  1:99713/140‑1 (MQ=255)                                                                                                                       |                                                                                                                                        TGATGCCGCTCGCCACCATCAGCGCGTCGTGGCTGCTCATCGCCGAGTGGGCAGGCCAGTCAACTTTCGGGCGCGGGGTGAAGAATTTATCCAGATTATCGGTATGCAGCTCAATCAGCCCCGGCAGCAGCCAGCCGCCTTCGCCGTCCATCGCCTCCGGCAGGCGGCTCTGGCTTTCGGCAAAGGCGCGGATTTCGCCGTTCTGCACCTCCAGCGAACCGCTTACCACCTCGTTTTCCAGCACCAGCTTAACG  >  NC_000913/4316951‑4317204

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: ‑