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_0009131,493,2190AG100.0% 148.8 / NA 42E250E (GAA→GAG) trgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/21);  total (21/21)

GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGC  >  NC_000913/1493086‑1493350
                                                                                                                                     |                                                                                                                                   
gcgcGTTTGTGCTTGCCCTGGTCATGA‑‑CTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggcgg                                                                                                                              >  1:128088/1‑139 (MQ=255)
    gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                            <  1:904839/139‑1 (MQ=255)
                ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                >  1:255650/1‑139 (MQ=255)
                   ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                            <  2:193577/124‑1 (MQ=255)
                   ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                            >  1:193577/1‑124 (MQ=255)
                          aCGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                            <  1:1080984/117‑1 (MQ=255)
                          aCGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                            >  2:1080984/1‑117 (MQ=255)
                           cGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                     <  2:379672/139‑1 (MQ=255)
                              tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAATAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCg                                                                                                                          >  2:765724/1‑115 (MQ=255)
                              tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCg                                                                                                                          <  1:765724/115‑1 (MQ=255)
                               gATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcgt                                                                                                 <  1:141972/139‑1 (MQ=255)
                                     aTTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTa                                                                                           <  1:986104/139‑1 (MQ=255)
                                       ttATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                        >  2:1080126/1‑108 (MQ=255)
                                       ttATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                        <  1:1080126/108‑1 (MQ=255)
                                          tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                >  1:290960/1‑97 (MQ=255)
                                          tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                <  2:290960/97‑1 (MQ=255)
                                             tGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                   <  1:672856/139‑1 (MQ=255)
                                               cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgc                                                                                 >  2:1247738/1‑139 (MQ=255)
                                                  cGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagc                                                                              >  1:102791/1‑139 (MQ=255)
                                                  cGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagc                                                                              >  1:153181/1‑139 (MQ=255)
                                                     cGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCCGCAtt                                                                           <  2:128088/139‑1 (MQ=255)
                                                              aTTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGa                                                                  >  2:565895/1‑139 (MQ=255)
                                                                tCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                                <  2:255650/139‑1 (MQ=255)
                                                                       cTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACt                                                                       <  1:625608/125‑1 (MQ=255)
                                                                       cTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACt                                                                       >  2:625608/1‑125 (MQ=255)
                                                                        tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAgg                                                        <  1:1247738/139‑1 (MQ=255)
                                                                               tGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGtt                                                 >  1:1034555/1‑139 (MQ=255)
                                                                                gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                        >  1:891720/1‑67 (MQ=255)
                                                                                gCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                                        <  2:891720/67‑1 (MQ=255)
                                                                                  cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGa                                              >  2:696870/1‑139 (MQ=255)
                                                                                    cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                             <  2:462447/106‑1 (MQ=255)
                                                                                    cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                             >  1:462447/1‑106 (MQ=255)
                                                                                            aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCgg                                    <  1:565895/139‑1 (MQ=255)
                                                                                              tGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                                  <  1:696870/139‑1 (MQ=255)
                                                                                                 aaaaaTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAagag                               >  1:586545/1‑139 (MQ=255)
                                                                                                   aaaTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGAt                             <  2:8909/139‑1 (MQ=255)
                                                                                                           aGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCgg                                    <  1:947180/124‑1 (MQ=255)
                                                                                                           aGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCgg                                    >  2:947180/1‑124 (MQ=255)
                                                                                                                            gaatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCa    >  2:994958/1‑139 (MQ=255)
                                                                                                                             aatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCa                                                                         >  1:642636/1‑69 (MQ=255)
                                                                                                                             aatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCa                                                                         <  2:642636/69‑1 (MQ=255)
                                                                                                                              atgatgAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGc  >  1:1295180/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                   
GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGC  >  NC_000913/1493086‑1493350

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