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%	183.1 / NA	50	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 (24/26);  total (24/26)

GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGG  >  NC_000913/1493086‑1493353                                                                                                                                      |                                                                                                                                       gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                   >  2:548495/1‑139 (MQ=255)         gtgCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                           <  2:1438962/139‑1 (MQ=255)                 ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                   <  1:1616925/139‑1 (MQ=255)                  cTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCg                                                                                                                  <  2:1293428/139‑1 (MQ=255)                    ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                <  2:1558128/139‑1 (MQ=255)                                       tttATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTAc                                                                                             <  1:573717/139‑1 (MQ=255)                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggcg                                                                                                                                  >  1:863307/1‑98 (MQ=255)                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggcg                                                                                                                                  <  2:863307/98‑1 (MQ=255)                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCa                                                                                         >  2:1022959/1‑139 (MQ=255)                                            ggTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCaa                                                                                        <  1:306070/139‑1 (MQ=255)                                              tGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                      <  1:1022959/139‑1 (MQ=255)                                                cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgc                                                                                    >  1:661838/1‑139 (MQ=255)                                                   cGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                <  2:1114266/108‑1 (MQ=255)                                                   cGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                >  1:1114266/1‑108 (MQ=255)                                                       ggATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCt                                                                                                            <  2:807648/108‑1 (MQ=255)                                                       ggATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCt                                                                                                            >  1:807648/1‑108 (MQ=255)                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                               <  2:605243/88‑1 (MQ=255)                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                               >  1:605243/1‑88 (MQ=255)                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcg                                                                                                     >  2:1065293/1‑114 (MQ=255)                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcg                                                                                                     <  1:1065293/114‑1 (MQ=255)                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCa                                                                            >  2:672840/1‑139 (MQ=255)                                                                     ccACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAg                                                               >  2:99512/1‑139 (MQ=255)                                                                           cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                                         >  1:987274/1‑139 (MQ=255)                                                                                tGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGtt                                                    <  1:99512/139‑1 (MQ=255)                                                                                 gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGt                                                              >  1:1332642/1‑128 (MQ=255)                                                                                 gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGt                                                              <  2:1332642/128‑1 (MQ=255)                                                                                 gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCACTTAAGTCGTCATTTACAGCAAAt                                                                                      >  1:7101/1‑104 (MQ=255)                                                                                 gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCACTTAAGTCGTCATTTACAGCAAAt                                                                                      <  2:7101/104‑1 (MQ=255)                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                     >  2:466304/1‑103 (MQ=255)                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                     <  1:466304/103‑1 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                               <  2:1050064/54‑1 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                        <  2:632442/77‑1 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                        >  1:632442/1‑77 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGa                                                                     >  1:1087008/1‑112 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGa                                                                     <  2:1087008/112‑1 (MQ=255)                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCAGACGATGAATGATGAGCCGGCGGGt                                                                                                                               >  1:1050064/1‑54 (MQ=255)                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                      >  2:1178595/1‑92 (MQ=255)                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                      <  1:1178595/92‑1 (MQ=255)                                                                                                  aaaaaTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                <  1:1530955/61‑1 (MQ=255)                                                                                                  aaaaaTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                >  2:1530955/1‑61 (MQ=255)                                                                                                      aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCAttt                                                                                               <  1:787891/74‑1 (MQ=255)                                                                                                      aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCAttt                                                                                               >  2:787891/1‑74 (MQ=255)                                                                                                              tGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                   >  1:1398764/1‑46 (MQ=255)                                                                                                              tGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                   <  2:1398764/46‑1 (MQ=255)                                                                                                              tGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                     >  1:857086/1‑76 (MQ=255)                                                                                                              tGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                     <  2:857086/76‑1 (MQ=255)                                                                                                               ggCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCa                     <  2:744607/139‑1 (MQ=255)                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCgg                                       >  2:1150802/1‑102 (MQ=255)                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTgg  >  2:1221283/1‑139 (MQ=255)                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGCGACTGTTCGACAGGGCACgg                                       <  1:1150802/102‑1 (MQ=255)                                                                                                                                      |                                                                                                                                       GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTGG  >  NC_000913/1493086‑1493353

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