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.7% 152.8 / ‑7.2 44E250E (GAA→GAGtrgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/20);  minor base C (0/1);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTG  >  NC_000913/1493085‑1493352
                                                                                                                                      |                                                                                                                                       
ggCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcgg                                                                                                                                     >  2:703009/1‑139 (MQ=255)
 gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                    >  1:931885/1‑139 (MQ=255)
  cgcgTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggcg                                                                                                                                   <  1:210116/139‑1 (MQ=255)
    cgTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCggg                                                                                                                                 <  1:703009/139‑1 (MQ=255)
     gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                                <  2:952290/139‑1 (MQ=255)
     gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                                <  2:895978/139‑1 (MQ=255)
               gCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGaaa                                                                                                                      >  1:521530/1‑139 (MQ=255)
                  cTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCg                                                                                                                   <  1:733462/139‑1 (MQ=255)
                    ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                                 >  2:602315/1‑139 (MQ=255)
                          gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTa                                                                                                           >  1:951692/1‑139 (MQ=255)
                                gATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAgtcgt                                                                                                     <  2:951692/139‑1 (MQ=255)
                                          aTGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                    >  2:496897/1‑98 (MQ=255)
                                          aTGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                                    <  1:496897/98‑1 (MQ=255)
                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAAt                                                                                                                          >  2:201202/1‑107 (MQ=255)
                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAAt                                                                                                                          <  1:201202/107‑1 (MQ=255)
                                           tGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTAGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCa                                                                                          >  1:569120/1‑139 (MQ=255)
                                               gCTACGTCTGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATg                                                                                      <  2:986712/139‑1 (MQ=255)
                                                cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgc                                                                                     >  2:940068/1‑139 (MQ=255)
                                                        gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑Ca                                                                           >  1:622756/1‑139 (MQ=255)
                                                          tCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACt                                                                         <  1:929725/139‑1 (MQ=255)
                                                               aTTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGa                                                                    <  2:931885/139‑1 (MQ=255)
                                                                 tCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATg                                                                  <  1:53022/139‑1 (MQ=255)
                                                                        cTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAg                                                           >  2:1045597/1‑139 (MQ=255)
                                                                        cTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAg                                                           >  1:446930/1‑139 (MQ=255)
                                                                         tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGACCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAgg                                                          <  2:674244/139‑1 (MQ=255)
                                                                         tGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTAGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAgg                                                          <  2:569120/139‑1 (MQ=255)
                                                                           cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGa                                                        <  1:166574/139‑1 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                                 >  2:459838/1‑108 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATgcagca                                                                                 <  1:459838/108‑1 (MQ=255)
                                                                                   cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGa                                                >  2:902577/1‑139 (MQ=255)
                                                                                     cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACa                                              >  2:300385/1‑139 (MQ=255)
                                                                                     cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACa                                              <  2:521530/139‑1 (MQ=255)
                                                                                     cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACa                                              >  1:744877/1‑139 (MQ=255)
                                                                                        aaCGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAggg                                           >  2:795303/1‑139 (MQ=255)
                                                                                          cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAGGgcg                                         >  2:30554/1‑139 (MQ=255)
                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACt                                                                         <  2:876248/104‑1 (MQ=255)
                                                                                             aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACt                                                                         >  1:876248/1‑104 (MQ=255)
                                                                                                       tCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTa                            >  2:325763/1‑139 (MQ=255)
                                                                                                                 cGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCAcc                  <  1:300385/139‑1 (MQ=255)
                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                       >  2:1038129/1‑55 (MQ=255)
                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAAt                                                                                       <  1:1038129/55‑1 (MQ=255)
                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATg                                                                  >  2:426718/1‑74 (MQ=255)
                                                                                                                                  atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATg                                                                  <  1:426718/74‑1 (MQ=255)
                                                                                                                                   tgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTg  >  1:567681/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATT‑‑CACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAGCGAAATTTCAGCTG  >  NC_000913/1493085‑1493352

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: