Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,771 G→A 8.0% intergenic (‑127/+87) PP_3414 ← / ← PP_3415 methyl‑accepting chemotaxis transducer/sensory box protein/LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,7710GA 8.0% 143.7 / 2.4 50intergenic (‑127/+87)PP_3414/PP_3415methyl‑accepting chemotaxis transducer/sensory box protein/LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

GATGGAAATCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCT  >  NC_002947/3865629‑3865919
                                                                                                                                              |                                                                                                                                                    
gATGGAAATCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCt                                                                                                                                                <  4:64559/149‑1 (MQ=255)
 aTGGAAATCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCtt                                                                                                                                               >  5:210899/1‑149 (MQ=255)
     aaaTCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCa                                                                                                                                           <  2:203158/149‑1 (MQ=255)
     aaaTCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCa                                                                                                                                           <  6:115409/149‑1 (MQ=255)
      aaTCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAg                                                                                                                                          >  1:303474/1‑149 (MQ=255)
                            ccTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGtt                                                                                                                    >  8:173660/1‑149 (MQ=255)
                              tATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCAc                                                                                                                           >  1:25105/1‑140 (MQ=255)
                              tATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCAc                                                                                                                           <  2:25105/140‑1 (MQ=255)
                              tATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTc                                                                                                                  >  2:230955/1‑149 (MQ=255)
                                tGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCGCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCg                                                                                                                 <  7:173660/148‑1 (MQ=255)
                                                cATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGc                                                                                                                                            >  1:166719/1‑105 (MQ=255)
                                                cATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGc                                                                                                                                            <  2:166719/105‑1 (MQ=255)
                                                 aTAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTgg                                                                                               >  2:100053/1‑149 (MQ=255)
                                                     aTCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGt                                                                                           >  1:190083/1‑149 (MQ=255)
                                                     aTCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCACCTGCTTCGCAGGAAAACCCGCTGCACCTTc                                                                                                                       >  6:21514/1‑121 (MQ=255)
                                                     aTCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCACCTGCTTCGCAGGAAAACCCGCTGCACCTTc                                                                                                                       <  5:26198/121‑1 (MQ=255)
                                                     aTCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCACCTGCTTCGCAGGAAAACCCGCTGCACCTTc                                                                                                                       <  5:21514/121‑1 (MQ=255)
                                                     aTCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCACCTGCTTCGCAGGAAAACCCGCTGCACCTTc                                                                                                                       >  6:26198/1‑121 (MQ=255)
                                                          cGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTTGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTcgcgc                                                                                      >  8:244110/1‑149 (MQ=255)
                                                              ggCAAGTCGTGGCATCTACAGTGCCTGTAGGAGCGGGTTCACCCGCGAACACGGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCgag                                                                                   <  8:271270/148‑1 (MQ=255)
                                                                         gCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTTGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCg                                                                       <  7:244110/149‑1 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCtgtg                                                                                                <  6:65655/122‑1 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCtgtg                                                                                                >  5:65655/1‑122 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCgag                                                                                   >  3:39901/1‑135 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCgag                                                                                   <  4:39901/135‑1 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCggg                                                                     >  6:243872/1‑149 (MQ=255)
                                                                           aTCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCggg                                                                     <  2:303474/149‑1 (MQ=255)
                                                                              tACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCAc                                                                  >  8:84040/1‑149 (MQ=255)
                                                                                cAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCActc                                                                <  5:243872/149‑1 (MQ=255)
                                                                                cAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCAct                                                                 <  3:26374/148‑1 (MQ=255)
                                                                                  gcgcCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCActct                                                               >  3:146059/1‑148 (MQ=255)
                                                                                    gcCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAg                                                            <  2:190083/149‑1 (MQ=255)
                                                                                     ccTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAgg                                                           <  8:186403/149‑1 (MQ=255)
                                                                                      cTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTgg                                                                                               >  2:18724/1‑112 (MQ=255)
                                                                                        gTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTgg                                                                         >  1:130157/1‑132 (MQ=255)
                                                                                        gTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTgg                                                                         <  2:130157/132‑1 (MQ=255)
                                                                                         tAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGct                                                       >  4:40623/1‑149 (MQ=255)
                                                                                                 ggTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCa                                               <  6:210899/149‑1 (MQ=255)
                                                                                                  gTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCacac                                              >  5:231424/1‑149 (MQ=255)
                                                                                                           cgAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTc                                                                                                                  <  6:242343/72‑1 (MQ=255)
                                                                                                           cgAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTc                                                                                                                  >  5:242343/1‑72 (MQ=255)
                                                                                                                acCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGc                                >  1:100754/1‑149 (MQ=255)
                                                                                                                   ggCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAAACATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGTGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCa                             <  1:100053/149‑1 (MQ=255)
                                                                                                                                ccATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTgcg                 >  4:33229/1‑148 (MQ=255)
                                                                                                                                  aTGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCAt              <  7:84040/149‑1 (MQ=255)
                                                                                                                                     ccccGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCtctc                                                    <  5:131187/108‑1 (MQ=255)
                                                                                                                                     ccccGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCtctc                                                    >  6:131187/1‑108 (MQ=255)
                                                                                                                                     ccccGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTc            <  3:40623/148‑1 (MQ=255)
                                                                                                                                     ccccGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCACACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTc            <  6:231424/148‑1 (MQ=255)
                                                                                                                                             cGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCg    >  6:297888/1‑148 (MQ=255)
                                                                                                                                              gCCTGCTTAGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCAATCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATACTGCTGCGCATTCGTCACCCGCt  >  8:212561/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GATGGAAATCGAACATGTCGCGGGGCTCCCTATGCGCGGTTTTCTGAGCATAGATCGCCGTTGGCAAGTCGTGGCATCTACAGCGCCTGTAGGAGCGGGTTCACCCGCGAACACCGGCGCAGCCAGTGCCATGCCCCGCGTCGCCTGCTTCGCAGGAAAACCCGCTGCACCTTCGTTTCGGGCCAATCATTCCTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCT  >  NC_002947/3865629‑3865919

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