Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,195,356 A→G 7.1% S103S (TCT→TCC)  PP_1940 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,195,3560AG 7.1% 157.3 / 2.7 56S103S (TCT→TCC) PP_1940methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (27/25);  new base G (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

TGCCTAGCGCGGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGATGTGG  >  NC_002947/2195217‑2195497
                                                                                                                                           |                                                                                                                                                
tGCCTAGCGCGGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTg                                                                                                                                          >  5:124384/1‑148 (MQ=255)
   cTAGCGCGGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGAATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCGGATGGCTGAAAt                                                                                                                                      >  1:73693/1‑149 (MQ=255)
      gcgcGGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATgg                                                                                                                                    <  6:58707/148‑1 (MQ=255)
        gcgGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCGTTCTCTAAACGCAGATGGCTGAAATGGAc                                                                                                                                  >  2:82908/1‑148 (MQ=255)
            tATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGAAGATGGCTGAAATGGACTAttt                                                                                                                             >  6:7530/1‑149 (MQ=255)
                  caacaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTACCAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCAcc                                                                                                                       <  4:303135/149‑1 (MQ=255)
                  caacaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCAcc                                                                                                                       >  6:102662/1‑149 (MQ=255)
                   aacaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGAATTGCTTGTAAAACTCTCTTGGCTTGGAAGACTATTCTTTCACTAAACGCAGATGGCTGAAATGGACGATTTTCCAGCt                                                                                                                      >  4:280192/1‑149 (MQ=255)
                    acaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGa                                                                                                                                   >  6:352742/1‑135 (MQ=255)
                    acaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGa                                                                                                                                   <  5:352742/135‑1 (MQ=255)
                    acaaATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTg                                                                                                                     >  8:39422/1‑149 (MQ=255)
                      aaaTATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGt                                                                                                                   >  7:274913/1‑149 (MQ=255)
                       aaTATTTCATCTTTCCCATCATGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGtt                                                                                                                  <  1:179844/149‑1 (MQ=255)
                              cATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAACTCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCATGGCTTGGAAGCCTATGCTTGCTCTAACCGCAGAGGGCTGACATGGAGTAGTTTCCACCTGGTTTAGCAGa                                                                                                           >  4:256961/1‑149 (MQ=255)
                                        aTCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCGGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCa                                                                                                 <  2:73693/149‑1 (MQ=255)
                                        aTCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGAATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCACATGGATGAAATGGAATATTTTCAACCTGGTTTAGCATAAgctgatcaa                                                                                                 >  1:263562/1‑144 (MQ=255)
                                               gCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑cg                                                                                       >  7:67076/1‑149 (MQ=255)
                                                cATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑cgc                                                                                      <  3:40921/149‑1 (MQ=255)
                                                cATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑cg                                                                                       <  7:267681/148‑1 (MQ=255)
                                                      gCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGa                                                                                                                                   <  7:204621/101‑1 (MQ=255)
                                                      gCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGa                                                                                                                                   >  8:204621/1‑101 (MQ=255)
                                                      gCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGaaaa                                                                                 <  5:102662/148‑1 (MQ=255)
                                                          gAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCAt                                                                            <  6:124384/149‑1 (MQ=255)
                                                           aaTACCAAATAAATCCCTAACAGCGCCTTATTTGCTTGTAAAACTGTCTTGGCTTTGAATACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCAGCTGGTTTAGCAGAAGGTGCTGCATAT‑‑GT‑CGCGAAAAGCCAt                                                                            <  2:263562/148‑1 (MQ=255)
                                                                cAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTAtttt                                                                                                                            >  2:303451/1‑98 (MQ=255)
                                                                cAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTAtttt                                                                                                                            <  1:303451/98‑1 (MQ=255)
                                                                cAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGACATAGcag                                                                       >  4:9044/1‑146 (MQ=255)
                                                                 aaaGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTa                                                                     <  8:274913/149‑1 (MQ=255)
                                                                     aaaTCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTgg                                                                 <  1:73720/149‑1 (MQ=255)
                                                                     aaaTCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTgg                                                                 >  2:73720/1‑149 (MQ=255)
                                                                                aGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGAGTTAGCTTTCTCTAAACGCAGATGGCTGAAATGTACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGt                                                      <  3:280192/149‑1 (MQ=255)
                                                                                    ccGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTg                                                                          >  7:307590/1‑125 (MQ=255)
                                                                                    ccGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTg                                                                          <  8:307590/125‑1 (MQ=255)
                                                                                              ttGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGga                                        >  7:143427/1‑149 (MQ=255)
                                                                                               tGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCGGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTAc                                                                    <  1:125340/120‑1 (MQ=255)
                                                                                               tGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCGGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTAc                                                                    >  2:125340/1‑120 (MQ=255)
                                                                                                   aaaCTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTAGTTTAGCAGAAGCTGCTGCATATGAGTGCGGGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGgat                                       >  1:115770/1‑148 (MQ=255)
                                                                                                     aCTGTCTTGGCTTGGAAGAATATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACTTGGTTTAGCAGAAGCTGCTTCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCaa                                 <  7:230874/149‑1 (MQ=255)
                                                                                                           ttggcttggAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTg                           >  1:308612/1‑149 (MQ=255)
                                                                                                            tggcttggAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTg                           >  2:101847/1‑148 (MQ=255)
                                                                                                              gcttggAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTa                         <  8:143427/148‑1 (MQ=255)
                                                                                                                  ggAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTTACGGCGGATGATCAAGGTTTGTATTTCa                    >  4:29651/1‑149 (MQ=255)
                                                                                                                  ggAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCa                    >  4:185138/1‑149 (MQ=255)
                                                                                                                   gAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTGCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCa                    <  5:7530/148‑1 (MQ=255)
                                                                                                                         tATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCa                    >  6:128721/1‑142 (MQ=255)
                                                                                                                         tATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCa                    <  5:128721/142‑1 (MQ=255)
                                                                                                                          aTGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCg                                                        >  1:215530/1‑105 (MQ=255)
                                                                                                                          aTGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCg                                                        <  2:215530/105‑1 (MQ=255)
                                                                                                                             cTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGAcggc                                           >  4:339144/1‑115 (MQ=255)
                                                                                                                             cTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGAcggc                                           <  3:339144/115‑1 (MQ=255)
                                                                                                                              tttCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCAt                                                          <  1:346038/99‑1 (MQ=255)
                                                                                                                              tttCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCAt                                                          >  2:346038/1‑99 (MQ=255)
                                                                                                                                 ctctAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGAt      <  2:308612/148‑1 (MQ=255)
                                                                                                                                 ctctAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGAt      <  1:101847/148‑1 (MQ=255)
                                                                                                                                    tAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGATGTgg  >  4:111558/1‑149 (MQ=255)
                                                                                                                                      aaCGCAGATGGCTGAAATGGCCTATTTTCCACCTGGGTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTCCTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGATGTgg  <  3:256961/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                
TGCCTAGCGCGGTATTACCAACAAATATTTCATCTTTCCCATCAAGTGCATTTGGCTCGAATACCAAAGAAATCCCTAACAGCCCCGGATTTGCTTGTAAAACTGTCTTGGCTTGGAAGACTATGCTTTCTCTAAACGCAGATGGCTGAAATGGACTATTTTCCACCTGGTTTAGCAGAAGCTGCTGCATAT‑‑GT‑CGCGAAAAGCCATTGCCGTACTGGTAGGCATCGGTGAACTGACGGCGGATGATCAAGGTTTGTATTTCACCTTGCGATTCGATGTGG  >  NC_002947/2195217‑2195497

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