Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,825,254	G→T	Y14* (TAC→TAA)	mutS ←	DNA mismatch repair protein MutS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,825,254	0	G	T	100.0%	145.3 / NA	41	Y14* (TAC→TAA)	mutS	DNA mismatch repair protein MutS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (16/25);  total (16/25)

CCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAGGCGCGTG  >  NC_002947/1825121‑1825376                                                                                                                                      |                                                                                                                           ccgACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTa                                                                                                                          <  2:707509/136‑1 (MQ=255)    tcTGACCGAGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActg                                                                                                                       <  1:60332/135‑1 (MQ=255)    aCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActg                                                                                                                       <  1:994614/136‑1 (MQ=255)    aCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActg                                                                                                                       <  1:17035/136‑1 (MQ=255)    aCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActg                                                                                                                       <  1:439887/136‑1 (MQ=255)      tGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGTAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActgct                                                                                                                     <  1:501262/136‑1 (MQ=255)      tGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTActgct                                                                                                                     <  1:290484/136‑1 (MQ=255)        aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACtgctgc                                                                                                                   <  1:238433/136‑1 (MQ=255)                 ggTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCggg                                                                                                          >  1:1267372/1‑136 (MQ=255)                     aGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGgtgt                                                                                                      >  2:303537/1‑136 (MQ=255)                     aGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCAACGGGgtgt                                                                                                      >  1:836956/1‑136 (MQ=255)                     aGGGTGATATCCAGCAGTATCGCGGCCTTTTTCGCATCTTCGTTGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGgtgt                                                                                                      <  2:797403/136‑1 (MQ=255)                      gggTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGgtgtg                                                                                                     >  2:873224/1‑136 (MQ=255)                       cataatatACAGCAGATTCGCGGCCTTTTTATAATCATAGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGgtgtgt                                                                                                    >  1:819138/5‑136 (MQ=255)                               ccAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGAGTGTGCTGaaat                                                                                            >  2:476218/1‑135 (MQ=255)                                cagcagTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGa                                                                                           >  1:1002957/1‑136 (MQ=255)                                    agTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTg                                                                                       >  1:109873/1‑136 (MQ=255)                                    agTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTg                                                                                       >  1:1053032/1‑136 (MQ=255)                                        tCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACAt                                                                                   <  2:693105/136‑1 (MQ=255)                                        tCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACAt                                                                                   <  2:109873/136‑1 (MQ=255)                                        tCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACAt                                                                                   <  1:865604/136‑1 (MQ=255)                                           cgGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCa                                                                                <  1:461415/136‑1 (MQ=255)                                           cgGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCa                                                                                >  1:1026751/1‑136 (MQ=255)                                            ggCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGAt                                                                                          >  2:1010059/1‑125 (MQ=255)                                            ggCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGAt                                                                                          <  1:1010059/125‑1 (MQ=255)                                              ccTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAggg                                                                             >  1:596505/1‑136 (MQ=255)                                                              gTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATc                                                             <  2:943633/136‑1 (MQ=255)                                                                    caTCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTg                                                       <  1:303862/135‑1 (MQ=255)                                                                                  ggCCCTTGCGTTATAACTTCTGCTGGTCCGGGTGCGGGTTTTTCTGCTTCCATTACTGCTGCTGCACGGGGGGGTGTGCGTAATTATCTGACATTCAGGGCCTAACAGGGGGTATCTGGTTGGCGTTGGTAAAccc                                         <  2:1230410/135‑1 (MQ=255)                                                                                  cGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAAccc                                         <  2:1053032/136‑1 (MQ=255)                                                                                        tGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCAACGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGTATCTGGTTGGCGATGGTAAACCCGTAATg                                   <  2:836956/136‑1 (MQ=255)                                                                                                         tGGTCCGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGtt                  <  2:1367577/136‑1 (MQ=255)                                                                                                             ccGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGa              >  2:484733/1‑136 (MQ=255)                                                                                                              ggggTGCTGGTTTTTGAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGAGATTCAGGGCCTTACAGGGGGGTTCTGGTTGGCGATGGTAAACCCGTAATAGTACAGGCTTTTTCGTTGCGAt             <  1:919566/135‑1 (MQ=255)                                                                                                              cGGGTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGAt             >  1:647942/1‑136 (MQ=255)                                                                                                                ggTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGc           >  2:242415/1‑136 (MQ=255)                                                                                                                 gTGCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCa          <  1:453736/136‑1 (MQ=255)                                                                                                                   gCTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAATGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAgg        <  2:450440/136‑1 (MQ=255)                                                                                                                    cTGGTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAGgc       <  1:242415/136‑1 (MQ=255)                                                                                                                      ggTTTTTCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAGgcgc     <  1:24024/136‑1 (MQ=255)                                                                                                                        tttttCAGCTTCCATTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACA‑GGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAGGCGCGTg  >  1:1278452/1‑136 (MQ=255)                                                                                                                                      |                                                                                                                           CCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAGAAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGTGCGGAAAGATCTGACATTCAGGGCCTTACAGGGGGGATCTGGTTGGCGATGGTAAACCCGTAATGGTACAGGCTTTTTCGTTGCGATGCAGGCGCGTG  >  NC_002947/1825121‑1825376

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