Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,493,774 A→G 20.0% L435L (TTA→TTG 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,7740AG20.0% 120.5 / 29.1 65L435L (TTA→TTGtrgmethyl‑accepting chemotaxis protein III, ribose and galactose sensor receptor
Reads supporting (aligned to +/- strand):  ref base A (23/29);  new base G (5/8);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 7.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

ATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAGGAAATCG  >  NC_000913/1493678‑1493871
                                                                                                |                                                                                                 
aTTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGAtt                                                                                                >  1:2018483/1‑100 (MQ=255)
aTTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAAtt                                                                                                <  2:1142923/100‑1 (MQ=255)
 tttGCCGTTGTCGCGAGCGAAGTACGGACACGCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTg                                                                                               <  2:1484984/100‑1 (MQ=255)
   tGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGACGGCTTGATCAGTGCATCAGTCAGGTTAATTGAc                                                                                             >  2:1261818/1‑100 (MQ=255)
       gTTGTCGCCAGCGAAGTACGGACACTCGCAAGTGGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTgg                                                                                         >  2:1584311/1‑100 (MQ=255)
          gTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGt                                                                                      >  1:659429/1‑100 (MQ=255)
          gTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGt                                                                                      >  1:1374083/1‑100 (MQ=255)
           tCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTc                                                                                     <  2:298082/100‑1 (MQ=255)
             gCCAGCGAAGAACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCgg                                                                                   >  2:1689492/1‑100 (MQ=255)
              ccAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGa                                                                                  <  2:2083437/100‑1 (MQ=255)
              ccAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGa                                                                                  <  2:1898419/100‑1 (MQ=255)
                aGCGAAGTCAGGACACTCGGAAGTCGCAGTGCTCATGCCGCGAAAGAGATTGAAGGCTTGATCAGTGTATCAGTCAGGTTAATTGACCTGGGGTCGGa                                                                                  <  2:824357/98‑1 (MQ=255)
                aGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATg                                                                                <  2:1979588/100‑1 (MQ=255)
                aGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATg                                                                                <  2:1983048/100‑1 (MQ=255)
                aGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGAt                                                                                 >  2:1883332/1‑99 (MQ=255)
                 gCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATg                                                                                <  2:2401234/99‑1 (MQ=255)
                  cGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAg                                                                              <  2:1877826/100‑1 (MQ=255)
                      gTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTgg                                                                          <  1:421312/100‑1 (MQ=255)
                       tACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGc                                                                         <  1:2360038/100‑1 (MQ=255)
                           gACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACg                                                                     >  2:1588529/1‑100 (MQ=255)
                           gACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACg                                                                     >  2:565876/1‑100 (MQ=255)
                           gACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACg                                                                     >  2:1512369/1‑100 (MQ=255)
                             cacTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGc                                                                   <  1:1968706/100‑1 (MQ=255)
                               ctctCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTTAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGcc                                                                  <  2:334868/99‑1 (MQ=255)
                                 cGCAAGTCGCAGCGCTCAGGCGGCGATAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCggg                                                               >  1:1462318/1‑100 (MQ=255)
                                 cGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCggg                                                               <  1:1291099/100‑1 (MQ=255)
                                 cGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCggg                                                               <  2:2071062/100‑1 (MQ=255)
                                       tCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAAcc                                                         <  2:2039658/100‑1 (MQ=255)
                                        cGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCa                                                        >  1:1350940/1‑100 (MQ=255)
                                                   ggcggcGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCGGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTAtt                                             >  2:1785479/1‑100 (MQ=255)
                                                   ggcggcGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTAtt                                             >  2:244426/1‑100 (MQ=255)
                                                      ggcgAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTCAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTAttgt                                           <  2:1678199/99‑1 (MQ=255)
                                                        cgAAAGGGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGa                                        <  2:1856347/100‑1 (MQ=255)
                                                         gAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGAt                                       >  1:1191449/1‑100 (MQ=255)
                                                          aaaGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATg                                      <  2:264019/100‑1 (MQ=255)
                                                          aaaGAGATTGAAGGCTTGATCAGTGAATCAGGCAGGTTAATTGGCCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGa                                        <  1:143392/98‑1 (MQ=255)
                                                           nagagATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGc                                     >  2:2213338/2‑100 (MQ=255)
                                                           aaGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGc                                     >  1:373634/1‑100 (MQ=255)
                                                           aaGAGATTGAAGGCTTGATAAATGAATAAATCAGGTTAATTGACATGAGGTCGGATAAGATAGCAACGACCGAGAAAACCATGAGCACTATTGTTAAt                                       >  2:713479/1‑98 (MQ=255)
                                                             gagaTTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCg                                   <  1:1588549/100‑1 (MQ=255)
                                                               gaTTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATAcc                                    >  1:1902772/1‑97 (MQ=255)
                                                                    aaGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAg                            >  1:1314834/1‑100 (MQ=255)
                                                                      ggCTTGATCAGTTAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAgtg                          <  2:2313618/100‑1 (MQ=255)
                                                                        cTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTc                        <  2:373631/100‑1 (MQ=255)
                                                                        cTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTc                        >  2:136211/1‑100 (MQ=255)
                                                                          tGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcac                      >  1:1759575/1‑100 (MQ=255)
                                                                           gATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcaca                     <  1:220900/100‑1 (MQ=255)
                                                                            aTCAGTGAATCAGTCAGGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcacac                    <  2:1660834/100‑1 (MQ=255)
                                                                            aTCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGTAACGGCCTGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcacac                    <  2:47796/100‑1 (MQ=255)
                                                                             tCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcacaca                   >  2:50060/1‑100 (MQ=255)
                                                                             tCAGTGAAGCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTcacaca                   <  1:453166/100‑1 (MQ=255)
                                                                              cAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGCGTGTCACACat                  >  2:1193784/1‑100 (MQ=255)
                                                                              cAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACat                  >  2:2264199/1‑100 (MQ=255)
                                                                               aGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACata                 <  1:1286450/100‑1 (MQ=255)
                                                                               aGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACata                 >  2:1237324/1‑100 (MQ=255)
                                                                               aGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACata                 >  1:502290/1‑100 (MQ=255)
                                                                                 tGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATc               <  1:861265/100‑1 (MQ=255)
                                                                                     tcagtcagGTTGATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGc           <  2:659423/100‑1 (MQ=255)
                                                                                     tcagtcagGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGc           <  2:1994262/100‑1 (MQ=255)
                                                                                     tcagtcagGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGc           <  1:785378/100‑1 (MQ=255)
                                                                                      cagtcagGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCa          >  2:2393348/1‑100 (MQ=255)
                                                                                        gtcagGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAgg        <  1:513099/100‑1 (MQ=255)
                                                                                            ggTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAGGAAAt    <  1:1316046/100‑1 (MQ=255)
                                                                                              ttAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAGGAAATcg  <  1:1644484/100‑1 (MQ=255)
                                                                                              ttAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAGGAAATcg  <  1:1883354/100‑1 (MQ=255)
                                                                                                |                                                                                                 
ATTTGCCGTTGTCGCCAGCGAAGTACGGACACTCGCAAGTCGCAGCGCTCAGGCGGCGAAAGAGATTGAAGGCTTGATCAGTGAATCAGTCAGGTTAATTGACCTGGGGTCGGATGAGGTGGCAACGGCCGGGAAAACCATGAGCACTATTGTTGATGCCGTCGCGAGTGTCACACATATCATGCAGGAAATCG  >  NC_000913/1493678‑1493871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: