Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,264,008 T→A 9.8% D180V (GAT→GTT) ‡ PP_2861 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,264,0080TA 9.8% 134.1 / 4.9 51D180V (GAT→GTT) ‡PP_2861methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (2/3);  total (24/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 = 1.00e+00

TGTTGCCCATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCGGGCGAGGGCGACC  >  NC_002947/3263862‑3264151
                                                                                                                                                  |                                                                                                                                               
tGTTGCCCATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATc                                                                                                                                                >  7:247849/1‑148 (MQ=255)
  ttGCCCATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCGTCGa                                                                                                                                              <  6:20771/148‑1 (MQ=255)
  ttGCCCATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGAc                                                                                                                                             <  5:21556/149‑1 (MQ=255)
      ccATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTa                                                                                                                                          >  4:137434/1‑148 (MQ=255)
           gCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGc                                                                                                                                     <  8:220585/148‑1 (MQ=255)
           gCCTGCACGGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCCCCAGCGCCATCAGCAGGGCAATGCCCCCGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCGATGACGTACAGGcc                                                                                                                                    <  3:13072/149‑1 (MQ=255)
              tGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGACCg                                                                                                                                  >  6:188775/1‑148 (MQ=255)
                   tGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCga                                                                                                                            >  1:72847/1‑149 (MQ=255)
                     cTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCgat                                                                                                                           >  1:64350/1‑148 (MQ=255)
                     cTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCga                                                                                                                            <  8:247849/147‑1 (MQ=255)
                      ttCCTGCAGCGGGCGGGCGATGCAGCGGGAAATCAGCATGACCACCAGCGCCCTCAGCACGGCACTGCCCACCCCAACCCGCGAAGCGTCTCGCATCTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTCCACGCCCGACCCgatga                                                                                                                         >  6:238415/1‑149 (MQ=255)
                                      gcgATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCtt                                                                                                         <  6:199373/149‑1 (MQ=255)
                                                        aGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACg                                                                                        <  2:64350/148‑1 (MQ=255)
                                                         gCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACTCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAg                                                                                      <  5:68178/149‑1 (MQ=255)
                                                         gCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAg                                                                                      >  3:117064/1‑149 (MQ=255)
                                                          cATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAgg                                                                                     >  2:15426/1‑149 (MQ=255)
                                                            tGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCCTCGACGTACACGCCCGCACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTc                                                                                   >  6:324391/1‑149 (MQ=255)
                                                               ccaccaGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTGCGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTg                                                                                <  6:32311/149‑1 (MQ=255)
                                                               ccaccaGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTg                                                                                >  8:62282/1‑149 (MQ=255)
                                                                caccaGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTAGGTCATCGACGTACACGCCCGAAcc                                                                                                                              >  1:99809/1‑102 (MQ=255)
                                                                caccaGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAAcc                                                                                                                              <  2:99809/102‑1 (MQ=255)
                                                                caccaGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTgg                                                                               <  3:137434/149‑1 (MQ=255)
                                                                   caGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCgatga                                                                                                                         <  2:143036/104‑1 (MQ=255)
                                                                   caGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCgatga                                                                                                                         >  1:143036/1‑104 (MQ=255)
                                                                             agcagGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCAc                                                                   <  4:117064/148‑1 (MQ=255)
                                                                                  ggCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGgcgc                                                             >  5:308155/1‑149 (MQ=255)
                                                                                  cgCAATGCCCACGCCAACCAGCGACGCGTCCCGCAACTGGCGAGCGAATTCGGCCTCTACGTCATCGACGAACACGCCCGAACCGATGAACCAGCCCCACGGCTAGAACAGCTGAATATACGAGGTCTTGGCTACCCGCTAACTGgcgc                                                             >  5:308152/2‑149 (MQ=255)
                                                                                          ccACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTAc                                                                           <  2:266898/127‑1 (MQ=255)
                                                                                          ccACGCCAACCAGCGAAGCGTCGCGCAACTCGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAAAAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTg                                                     >  4:86840/1‑149 (MQ=255)
                                                                                           cACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGt                                                                                                                                           <  1:8324/62‑1 (MQ=255)
                                                                                           cACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGt                                                                                                                                           >  2:8324/1‑62 (MQ=255)
                                                                                            aCGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTggg                                                   >  4:101549/1‑149 (MQ=255)
                                                                                                 aaCCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCga                                                                                                                            <  1:166689/71‑1 (MQ=255)
                                                                                                 aaCCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCga                                                                                                                            >  2:166689/1‑71 (MQ=255)
                                                                                                     aGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATa                                          <  1:15426/149‑1 (MQ=255)
                                                                                                             gTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCg                                  >  3:266961/1‑149 (MQ=255)
                                                                                                                gcgcAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCGCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCAGGc                               >  5:185675/1‑149 (MQ=255)
                                                                                                                    aaCTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGc                           <  4:266961/149‑1 (MQ=255)
                                                                                                                             gcgaATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCctgctg                  <  6:308155/149‑1 (MQ=255)
                                                                                                                             gcgaATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCctgctg                  <  7:62282/149‑1 (MQ=255)
                                                                                                                              tgaATTCGTCCTGTACGTCTTCGACGTAAACGCCCGAACCGATGATCCAGCCCCACGGGTTGAACAGCTGAATTTACGTGGTCTTGGATTCGGTCTCACTGGCGCTGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGGGTCctgctg                  <  6:308152/147‑1 (MQ=255)
                                                                                                                                    cGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTggcgggc            <  5:324391/148‑1 (MQ=255)
                                                                                                                                      gcgtGTACGTCGATGACGTACAGGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCgggcgag         >  8:123026/4‑149 (MQ=255)
                                                                                                                                       cgtGTACGTCGATGACGTACAGGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCgggcgagg        <  7:123026/147‑1 (MQ=255)
                                                                                                                                        gtGTACGTCGATGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCct                      >  3:143823/2‑134 (MQ=255)
                                                                                                                                        gtGTACGTCGATGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTCGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCct                      <  4:143823/133‑1 (MQ=255)
                                                                                                                                        cTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCgggcgaggg       >  2:6273/1‑149 (MQ=255)
                                                                                                                                          gTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCgggcgagg        <  1:147254/146‑1 (MQ=255)
                                                                                                                                          gTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACCGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCgggcgagg        >  2:147254/1‑146 (MQ=255)
                                                                                                                                             cGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCGGGCGAGGGCGAcc  <  3:101549/149‑1 (MQ=255)
                                                                                                                                             cGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCGGGCGAGGGCGAcc  <  1:6273/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGTTGCCCATCGCCTGCACTGCTTCCTGCAGCGGGCGGGCGATGCTGCGGGCAATCAGCATGACCACCAGCGCCATCAGCAGGGCAATGCCCACGCCAACCAGCGAAGCGTCGCGCAACTGGCGAGCGAATTCGGCCTGTACGTCATCGACGTACACGCCCGAACCGATGATCCAGCCCCACGGCTTGAACAGCTGAATATACGAGGTCTTGGCTACCGGCTCACTGGCGCCGGGCTTGGGCCAGCGATAGTTGACCGGGCCGGCGTCCTGCTGGCGGGCGAGGGCGACC  >  NC_002947/3263862‑3264151

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: