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% 151.5 / ‑5.9 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 (22/21);  minor base C (0/1);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAG  >  NC_000913/1493086‑1493338
                                                                                                                                     |                                                                                                                       
gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    >  2:58598/1‑139 (MQ=255)
gcgcGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    <  1:58598/139‑1 (MQ=255)
    gTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                <  2:1340044/139‑1 (MQ=255)
        gtgCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                            <  2:1178963/139‑1 (MQ=255)
                   ggTCATGACGCTGATAACATTTATGGTGCTACGTCGGACCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGt                                                                                                 >  2:500476/1‑139 (MQ=255)
                         gACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTa                                                                                           >  1:1044108/1‑139 (MQ=255)
                                         aTGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGc                                                                           <  2:1490588/139‑1 (MQ=255)
                                           ggTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCaa                                                                         <  1:1216855/139‑1 (MQ=255)
                                               cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    <  1:335104/92‑1 (MQ=255)
                                               cTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    >  2:335104/1‑92 (MQ=255)
                                                       gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                <  1:997254/88‑1 (MQ=255)
                                                       gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCg                                                                                                   >  2:786161/1‑101 (MQ=255)
                                                       gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCg                                                                                                   <  1:786161/101‑1 (MQ=255)
                                                       gATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGACGATCTGACGATGAATGATGAGCCGGCGGGt                                                                                                                >  2:997254/1‑88 (MQ=255)
                                                         tCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGc                                                                           <  1:1317608/123‑1 (MQ=255)
                                                         tCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGc                                                                           >  2:1317608/1‑123 (MQ=255)
                                                             cATTCGTCCACTGAAACATGCCCAACAACGGTTTGAAAAAATCGCCAGTGGCGATCTGACGAGGAAGGAGGGCCCGGCGGGTGGTAATGAAATCGGTCCCTTAAGCGGTCATTACAGCAAATGCAGCATTCACTGGGGa                                                      <  1:691002/139‑1 (MQ=255)
                                                                       cTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTACAACAAATGGAGCATTCACTGGGGGTGGCAGTAgg                                            >  2:691002/1‑139 (MQ=255)
                                                                          cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                            <  1:747765/73‑1 (MQ=255)
                                                                          cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTa                                                                                                            >  2:747765/1‑73 (MQ=255)
                                                                          cAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGa                                          <  2:913757/139‑1 (MQ=255)
                                                                                 ccGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCg                                   >  1:559814/1‑139 (MQ=255)
                                                                                  cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                    <  1:31249/121‑1 (MQ=255)
                                                                                  cGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATg                                                    >  2:31249/1‑121 (MQ=255)
                                                                                    cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAgg                                            >  1:908399/1‑127 (MQ=255)
                                                                                    cacaACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAgg                                            <  2:908399/127‑1 (MQ=255)
                                                                                         cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    <  2:1033369/50‑1 (MQ=255)
                                                                                         cGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCcggc                                                                                                                    >  1:1033369/1‑50 (MQ=255)
                                                                                            aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                     >  2:712217/1‑110 (MQ=255)
                                                                                            aTTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGAt                                                     <  1:712217/110‑1 (MQ=255)
                                                                                              tGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                      <  2:559814/139‑1 (MQ=255)
                                                                                                     aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                      <  1:514845/132‑1 (MQ=255)
                                                                                                     aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGaa                      >  2:514845/1‑132 (MQ=255)
                                                                                                     aTCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGAttt               >  1:328875/1‑139 (MQ=255)
                                                                                                      tCGCCAGTGGCGATCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTa              >  2:1263145/1‑139 (MQ=255)
                                                                                                                  aTCTGACGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAg  >  1:433900/1‑139 (MQ=255)
                                                                                                                       aCGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCa      <  2:356467/130‑1 (MQ=255)
                                                                                                                       aCGATGAATGATGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCa      >  1:356467/1‑130 (MQ=255)
                                                                                                                          atgaatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGAt                 <  2:168351/116‑1 (MQ=255)
                                                                                                                          atgaatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGAt                 >  1:168351/1‑116 (MQ=255)
                                                                                                                             aatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGc                                                                           >  2:634255/1‑55 (MQ=255)
                                                                                                                             aatgaTGAGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGc                                                                           <  1:634255/55‑1 (MQ=255)
                                                                                                                                 atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAg  >  2:1389153/1‑124 (MQ=255)
                                                                                                                                 atgaGCCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAg  <  1:1389153/124‑1 (MQ=255)
                                                                                                                                     |                                                                                                                       
GCGCGTTTGTGCTTGCCCTGGTCATGACGCTGATAACATTTATGGTGCTACGTCGGATCGTCATTCGTCCACTGCAACATGCCGCACAACGGATTGAAAAAATCGCCAGTGGCGATCTGACGATGAATGATGAACCGGCGGGTCGTAATGAAATCGGTCGCTTAAGTCGTCATTTACAGCAAATGCAGCATTCACTGGGGATGACAGTAGGGACTGTTCGACAGGGCGCGGAAGAGATTTATCGTGGCACCAG  >  NC_000913/1493086‑1493338

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: