Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,493,219	A→G	E250E (GAA→GAG)	trg →	methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,493,219	0	A	G	100.0%	115.7 / NA	33	E250E (GAA→GAG)	trg	methyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/19);  total (14/19)

GTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGG  >  NC_000913/1493090‑1493353                                                                                                                                  |                                                                                                                                       gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                               <  2:280036/139‑1 (MQ=255) gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                               <  2:566450/139‑1 (MQ=255)  tttGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTc                                                                                                                              <  2:244836/139‑1 (MQ=255)     gtgCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                           <  1:136061/139‑1 (MQ=255)                      gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTa                                                                                                          >  2:601802/1‑139 (MQ=255)                      gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTa                                                                                                          >  2:611731/1‑139 (MQ=255)                        cGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                        <  1:673452/139‑1 (MQ=255)                           tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcg                                                                                                     <  1:267837/139‑1 (MQ=255)                                    ttATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACa                                                                                            <  2:669853/139‑1 (MQ=255)                                       tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTGAGTCGTCATTTACAGCa                                                                                         >  1:475189/1‑139 (MQ=255)                                       tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                        <  1:230085/124‑1 (MQ=255)                                       tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                        >  2:230085/1‑124 (MQ=255)                                       tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCa                                                                                         >  2:61742/1‑139 (MQ=255)                                       tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCa                                                                                         >  2:421835/1‑139 (MQ=255)                                               cGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagc                                                                                 <  1:601802/139‑1 (MQ=255)                                                gTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                                <  2:66644/139‑1 (MQ=255)                                                    gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCa                                                                            <  2:648336/139‑1 (MQ=255)                                                    gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGCCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCa                                                                            <  2:475189/139‑1 (MQ=255)                                                       cGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTg                                                                         <  1:640736/139‑1 (MQ=255)                                                          cATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTgggg                                                                      <  1:446486/139‑1 (MQ=255)                                                                       cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                                   >  2:157757/1‑129 (MQ=255)                                                                       cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                                   <  1:157757/129‑1 (MQ=255)                                                                         aCATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTg                                                                         >  1:638151/1‑121 (MQ=255)                                                                         aCATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTg                                                                         <  2:638151/121‑1 (MQ=255)                                                                             gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcg                                                                                                                                     >  2:210200/1‑57 (MQ=255)                                                                             gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcg                                                                                                                                     <  1:210200/57‑1 (MQ=255)                                                                              ccGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCg                                                  >  1:167876/1‑139 (MQ=255)                                                                               cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGa                                                 >  1:262321/1‑139 (MQ=255)                                                                                 cacaACGGATTGAAATAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                                    <  1:484607/118‑1 (MQ=255)                                                                                 cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                                    >  2:484607/1‑118 (MQ=255)                                                                                                                      gATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTgggg                                                                      >  2:140632/1‑79 (MQ=255)                                                                                                                      gATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTgggg                                                                      <  1:140632/79‑1 (MQ=255)                                                                                                                              atgaGCCGGCAGGTCATAATGAAATCGGTCGCTTAAGACGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTgg  >  2:515936/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                       GTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGG  >  NC_000913/1493090‑1493353

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