Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	3,788,632	C→G	M1M (ATG→ATC) †	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,788,632	0	C	G	100.0%	127.8 / NA	38	M1M (ATG→ATC) †	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (18/20);  total (18/20)

TCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGAATGATCAGCATCA  >  NC_002947/3788492‑3788760                                                                                                                                              |                                                                                                                                 tcgAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGaa                                                                                                                           >  8:368567/1‑148 (MQ=255)   gagaAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAAtt                                                                                                                         <  2:70021/148‑1 (MQ=255)   gagaAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAAt                                                                                                                          <  8:389630/147‑1 (MQ=255)      aaaGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTc                                                                                                                       >  3:297589/1‑147 (MQ=255)         gCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCt                                                                                                                    >  3:50349/1‑147 (MQ=255)                   tCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGtttt                                                                                                         >  2:93979/1‑148 (MQ=255)                   tCGATCCAGCGAGGGGCCCCTGGACTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGttt                                                                                                          >  3:422322/1‑147 (MQ=255)                         cAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTGCCTTCACTTGTTTTTATAAt                                                                                                   >  7:354051/1‑148 (MQ=255)                         cAGCGAGGGGCCCCTGGACTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTtata                                                                                                     <  4:422322/146‑1 (MQ=255)                          aGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAAtt                                                                                                  >  7:370623/1‑148 (MQ=255)                               ggggCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTc                                                                                              >  4:55285/1‑147 (MQ=255)                                gggCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTc                                                                                              >  5:416739/1‑146 (MQ=255)                                 ggCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTa                                                                                           >  8:227182/1‑148 (MQ=255)                                    cccTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACaa                                                                                        <  7:368567/148‑1 (MQ=255)                                        ggCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACt                                                                                     <  2:57824/147‑1 (MQ=255)                                         gCCTCCCCGGCCGGCTCTCTTCAAGCGGGCTGGGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTg                                                                                    >  1:57824/1‑148 (MQ=255)                                          ccTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTg                                                                                    <  4:50349/146‑1 (MQ=255)                                              cccGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATaa                                                                               >  6:255423/1‑147 (MQ=255)                                               ccggccggCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAAct                                                                             >  7:98449/1‑148 (MQ=255)                                                          tctTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACg                                                                    <  7:227182/146‑1 (MQ=255)                                                          tctTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGa                                                                   >  4:165989/1‑147 (MQ=255)                                                           ctTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACg                                                                    <  4:178069/145‑1 (MQ=255)                                                               aaGCGGGCT‑GGCCGCGTAGGCGGAATGGCCCAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTAGAATTGTTCTTACAAACTGGATAACTCTGTTAACGattatt                                                              >  5:196849/1‑147 (MQ=255)                                                               aaGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTCATGAACTGATTTTCCCTGAAGTTGTTTTTTTAATTGTTCTTACAAAATGTATAAAAGTGTTAAGGataatt                                                              >  4:240038/1‑147 (MQ=255)                                                                aGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTAtttt                                                            <  2:16284/148‑1 (MQ=255)                                                                   gggCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTg                                                          <  6:416739/147‑1 (MQ=255)                                                                       t‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGCTGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTaa                                                      <  4:297589/147‑1 (MQ=255)                                                                           ccGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCat                                                  <  8:370623/148‑1 (MQ=255)                                                                               gTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATgc                                              <  8:98449/148‑1 (MQ=255)                                                                                tAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATgc                                              <  4:378523/147‑1 (MQ=255)                                                                                    cGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCg                                          >  6:369931/1‑147 (MQ=255)                                                                                                    aaGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGcc                          >  3:335498/1‑147 (MQ=255)                                                                                                     aGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGcc                          <  5:255423/146‑1 (MQ=255)                                                                                                       cACCAAGATGGGCATGGATAACCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCgtg                       <  3:240038/147‑1 (MQ=255)                                                                                                               ggggCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGa               <  3:47832/147‑1 (MQ=255)                                                                                                                gggCAGGCATACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGaa              <  3:55285/147‑1 (MQ=255)                                                                                                                  gcaggcaTACCCCGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGAATGa           <  1:93979/148‑1 (MQ=255)                                                                                                                            cccGTTCAACGGCTTGATGATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGAATGATCAGCATCa  <  4:335498/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                 TCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCTCTTCAAGCGGGCT‑GGCCGCGTAGGCGGAATGGCACAGCGAAAGCACCATGAGGGGCAGGCATACCCCGTTCAACGGCTTGATCATCAGAATTTCCCTTCACTTGTTTTTATAATTGTTCTTACAAACTGGATAACTCTGTTAACGATTATTTTTGTTAACCATATGCGCCGTCAATAGCGTTTTGCCGTGTACAAGGAATGATCAGCATCA  >  NC_002947/3788492‑3788760

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