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,2190AG97.1% 116.6 / NA 35E250E (GAA→GAGtrgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (16/18);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGC  >  NC_000913/1493085‑1493350
                                                                                                                                      |                                                                                                                                   
ggCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcgg                                                                                                                                 <  2:1052368/139‑1 (MQ=255)
 gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                >  2:350854/1‑139 (MQ=255)
    cgTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCggg                                                                                                                             <  1:223725/139‑1 (MQ=255)
                 ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                >  1:239977/1‑139 (MQ=255)
                 ccTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATc                                                                                                                >  1:413576/1‑139 (MQ=255)
                    ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcgg                                                                                                                                 >  1:1014771/1‑119 (MQ=255)
                    ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcgg                                                                                                                                 <  2:1014771/119‑1 (MQ=255)
                          gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTa                                                                                                       >  2:162523/1‑139 (MQ=255)
                            cGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                     <  2:944522/139‑1 (MQ=255)
                            cGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAg                                                                                                     <  2:413576/139‑1 (MQ=255)
                               tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcg                                                                                                  >  2:1068914/1‑139 (MQ=255)
                               tGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcg                                                                                                  >  2:1040995/1‑139 (MQ=255)
                                     cATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCAttt                                                                                            <  1:350854/139‑1 (MQ=255)
                                               gCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                  >  2:1307903/1‑139 (MQ=255)
                                                cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtc                                                                                                   <  2:1020699/121‑1 (MQ=255)
                                                cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtc                                                                                                   >  1:1020699/1‑121 (MQ=255)
                                                           cGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTg                                                                      <  2:151017/139‑1 (MQ=255)
                                                                         tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                             <  1:1088557/118‑1 (MQ=255)
                                                                         tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                             >  2:1088557/1‑118 (MQ=255)
                                                                           cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                                      <  1:162523/139‑1 (MQ=255)
                                                                           cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCAGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                                      <  2:239977/139‑1 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCaa                                                                                     <  1:1199541/100‑1 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCaa                                                                                     >  2:1199541/1‑100 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGa                                                                  <  1:567381/119‑1 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGa                                                                  >  2:567381/1‑119 (MQ=255)
                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGgcg                                       >  1:1131039/1‑139 (MQ=255)
                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                   <  2:404186/92‑1 (MQ=255)
                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                   >  1:404186/1‑92 (MQ=255)
                                                                                              ttGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGa                                   <  1:1307903/139‑1 (MQ=255)
                                                                                               tGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                                  <  2:1131039/139‑1 (MQ=255)
                                                                                                      aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                                >  1:1194359/1‑102 (MQ=255)
                                                                                                      aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                                <  2:1194359/102‑1 (MQ=255)
                                                                                                                      tGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGa           <  1:1068914/139‑1 (MQ=255)
                                                                                                                              aatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTAACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAg   >  1:648904/1‑139 (MQ=255)
                                                                                                                               atgatgAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGc  >  2:679595/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGC  >  NC_000913/1493085‑1493350

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