Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,019,896	G→T	A72A (GCC→GCA)	PP_4423 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,019,896	0	G	T	100.0%	82.0 / NA	25	A72A (GCC→GCA)	PP_4423	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/11);  total (14/11)

GTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAG  >  NC_002947/5019753‑5020034                                                                                                                                                |                                                                                                                                           gTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCa                                                                                                                                        >  1:132084/1‑148 (MQ=255)     cGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCtt                                                                                                                                    >  3:21365/1‑148 (MQ=255)          gCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGa                                                                                                                               <  6:72353/148‑1 (MQ=255)           ccTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAg                                                                                                                              >  6:92135/1‑148 (MQ=255)                 cTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAAt                                                                                                                        <  8:9528/148‑1 (MQ=255)                           aCGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTc                                                                                                               >  3:133545/1‑147 (MQ=255)                              aTTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATAc                                                                                                           <  1:164444/148‑1 (MQ=255)                              aTTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATAc                                                                                                           <  4:54484/148‑1 (MQ=255)                              aTTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGAGAGTTCATGGAGGGCAAAATAATGCTTCGCCCTCTGACTGTTTCGAGGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTGCAGAc                                                                                                           >  4:130730/1‑148 (MQ=255)                                tCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGCTTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTc                                                                                                         >  7:469/1‑148 (MQ=255)                                tCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTc                                                                                                         >  7:2748/1‑148 (MQ=255)                                     ttCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGtctc                                                                                                    <  7:144627/148‑1 (MQ=255)                                        ccGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTg                                                                                                                                >  7:4246/1‑117 (MQ=255)                                         cgcgATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGa                                                                                                                               <  8:4246/117‑1 (MQ=255)                                                aaTCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTgccg                                                                                         >  4:17896/1‑148 (MQ=255)                                                    ggAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAg                                                                                     >  4:312/1‑148 (MQ=255)                                                       gagGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTga                                                                                  <  1:60147/148‑1 (MQ=255)                                                          gTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCTAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTgagga                                                                               >  2:128595/1‑148 (MQ=255)                                                                                   gATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGtttt                                                       >  3:91159/1‑147 (MQ=255)                                                                                        ttCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGa                                                 <  4:149162/148‑1 (MQ=255)                                                                                           aTGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACg                                               <  4:21365/147‑1 (MQ=255)                                                                                               aGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGAt                                          <  5:92135/148‑1 (MQ=255)                                                                                                                   ccTCTGACTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCAcccc                      >  5:128851/1‑148 (MQ=255)                                                                                                                          cTGTTTCGATGTCCAGTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCaga               <  3:312/148‑1 (MQ=255)                                                                                                                                         gTTGATTTGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAg  >  2:119096/1‑146 (MQ=255)                                                                                                                                                |                                                                                                                                           GTCCCGTCAGCCTTACCTGGGAAGGTACGATTCAGATTCCCGCGATCAATCGGAGAGGTTCTGGAGGCGTTTTTAAAAGCAGGATAGTTCATGAAGGGCAAAATAATGATTCGCCCTCTGACTGTTTCGATGTCCAGTTGATTGGCCAGCTTGGTGAGTGCAATCGGGCCTTCATACTCGTCTCCATGATTGCCGCCAGTGAGGAGCACAGTCGGGCCGATCCCGTTTTTCATGACGGTGATAGGAATCATGATCGCACCCCACGCAGAGTCATCACGTGAG  >  NC_002947/5019753‑5020034

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