Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 381,859 T→A 11.8% D230V (GAC→GTC)  PP_0317 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947381,8590TA11.8% 143.5 / 10.9 51D230V (GAC→GTC) PP_0317methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (18/27);  new base A (2/4);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

TTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGGC  >  NC_002947/381714‑382004
                                                                                                                                                 |                                                                                                                                                 
ttctGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGc                                                                                                                                               <  1:174713/150‑1 (MQ=255)
 tctGCCGCCGAGGCCAGCTGTGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCg                                                                                                                                              <  1:782250/150‑1 (MQ=255)
   tGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGAc                                                                                                                                            <  2:1345849/150‑1 (MQ=255)
     ccgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGAcct                                                                                                                                          >  1:390910/1‑148 (MQ=255)
      cgccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACCTg                                                                                                                                         <  1:280732/150‑1 (MQ=255)
       gccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACCTga                                                                                                                                        <  2:750865/150‑2 (MQ=255)
       gccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGACGGCGACCTga                                                                                                                                        <  2:390910/150‑2 (MQ=255)
        ccgAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCt                                                                                                                                       >  1:1166670/1‑150 (MQ=255)
          gAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCg                                                                                                                                     <  1:174757/150‑1 (MQ=255)
                    gggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                           <  2:999394/150‑1 (MQ=255)
                     ggTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGACCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAg                                                                                                                          <  2:1099388/150‑1 (MQ=255)
                      gTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGt                                                                                                                         <  2:212506/150‑1 (MQ=255)
                               ccGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCg                                                                                                                >  2:957525/1‑150 (MQ=255)
                                       cggctgctAGGTATCGCGCAGCTTGTCCTGCATGCCCTCGATGTCGCGCAACAGCCGGGCGGCTTCGTCCGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCtt                                                                                                        <  1:1384535/148‑1 (MQ=255)
                                          ctgcAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                                     <  2:1357419/150‑1 (MQ=255)
                                           tgcAAGGTCGAGCGCAGCGTGTCCTTCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACGTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACg                                                                                                    >  2:1156283/1‑150 (MQ=255)
                                                ggTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACt                                                                                               >  2:795868/1‑150 (MQ=255)
                                                    tCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCgg                                                                                           <  1:926690/150‑1 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGc                                                                                             >  2:844744/1‑146 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGc                                                                                             <  1:844744/146‑1 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGt                                                                                         >  2:1069213/1‑150 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGTCGCCGACGGCGAcctga                                                                                                                                        >  2:633018/1‑98 (MQ=255)
                                                      gcgcAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGTCGCCGACGGCGACCTga                                                                                                                                        <  1:633018/102‑2 (MQ=255)
                                                           gCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGAGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTcagca                                                                                    <  1:1325277/150‑1 (MQ=255)
                                                             ttGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAgg                                                                                  >  1:1122008/1‑150 (MQ=255)
                                                                tCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATa                                                                                                 >  2:146015/1‑132 (MQ=255)
                                                                tCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATa                                                                                                 <  1:146015/132‑1 (MQ=255)
                                                                tCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAg                                                                               <  2:839601/150‑1 (MQ=255)
                                                                  cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGc                                                                             <  2:1166670/150‑1 (MQ=255)
                                                                  cTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGc                                                                             >  2:1319884/1‑150 (MQ=255)
                                                                   tGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGGCGGTGTCCTCTGCGCGGCTGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCGCCAGGACGGcc                                                                            >  2:487907/1‑149 (MQ=255)
                                                                           cgcCATGGCGCGCAACAGCCGGGCGGCTTCGTTGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGaa                                                                    <  1:1215599/150‑1 (MQ=255)
                                                                            gcCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAg                                                                   <  1:1319884/150‑1 (MQ=255)
                                                                                tGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCa                                                               <  2:1273966/150‑1 (MQ=255)
                                                                                     cgcAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACg                                                          <  2:1122008/150‑1 (MQ=255)
                                                                                        aaCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                           <  1:837081/82‑1 (MQ=255)
                                                                                        aaCAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGa                                                                                                                           >  2:837081/1‑82 (MQ=255)
                                                                                            gCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                                     >  2:547/1‑100 (MQ=255)
                                                                                            gCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGAc                                                                                                     <  1:547/100‑1 (MQ=255)
                                                                                                         gttggTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGaa                                                                                <  2:356837/105‑1 (MQ=255)
                                                                                                         gttggTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGaa                                                                                >  1:356837/4‑108 (MQ=255)
                                                                                                         gtcgGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGcc                                      <  1:1074078/150‑1 (MQ=255)
                                                                                                                cccTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggc                                                     >  2:1150591/1‑128 (MQ=255)
                                                                                                                cccTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTggc                                                     <  1:1150591/128‑1 (MQ=255)
                                                                                                                 ccTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCa                              >  1:320587/1‑150 (MQ=255)
                                                                                                                         cgGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGaaa                      <  1:795868/150‑1 (MQ=255)
                                                                                                                                 gTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGc              >  2:1395052/1‑150 (MQ=255)
                                                                                                                                 gTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGc              >  2:842171/1‑150 (MQ=255)
                                                                                                                                      tCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGc                          <  2:121018/133‑1 (MQ=255)
                                                                                                                                      tCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGc                          >  1:121018/1‑133 (MQ=255)
                                                                                                                                             gCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTggc  <  1:1069213/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCTGCCGCCGAGGCCAGCTGGGTGGCGGAGCCGGCGATCTGCTGCAAGGTATCGCGCAGCTTGTCCTGCATGCGCGCCATGGCGCGCAACAGCCGGGCGGCTTCGTCGGTGCCCTCTGCGCGGATGACGTGGGTCAGGTCGCCGTCGGCGACTTGCTCGGCGCATTTGAGTGCTTCTTCGATCGGCTTGACGATACTGCGGGTCAGCAGGAAGGCACAGGCGAAGGTCAGTACGGTGGCGGCAACCAGCAGGCCGATCACCAGGGCGAAAGCGCCGGCATACTGGTTGGC  >  NC_002947/381714‑382004

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