Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,403,591 T→C 17.0% intergenic (‑57/+24) PP_1227 ← / ← PP_1228 membrane protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,403,5910TC17.0% 91.4 / 14.2 49intergenic (‑57/+24)PP_1227/PP_1228membrane protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base C (5/3);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

ACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCC  >  NC_002947/1403453‑1403734
                                                                                                                                              |                                                                                                                                                
aCAGGGTGGCGGCGATCGAGGGGCGTAATAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGGTGGAATAGGCATAAGGTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCACGCCCGCGAag                                                                                                                                        <  3:188082/149‑1 (MQ=39)
 cAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGA‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCACACCCGCGaa                                                                                                                                         >  1:61674/1‑147 (MQ=255)
 cAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGA‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCACACCCGCGAaga                                                                                                                                       <  2:61674/149‑1 (MQ=255)
   gggTGGCTGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTTGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAgg                                                                                                                                     <  8:125031/149‑1 (MQ=255)
              aTCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGc                                                                                                                          >  7:153380/1‑149 (MQ=255)
               tCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGc                                                                                                                          >  2:24424/1‑148 (MQ=255)
                          taaaGCCCTTGTTCGCTGGTGATGTCGAAGAAAGGGGAGCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAAc                                                                                                              <  3:297677/148‑1 (MQ=255)
                            aaGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAAc                                                                                                              >  4:297677/1‑147 (MQ=255)
                                ccTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTc                                                                                                         >  8:77646/1‑148 (MQ=255)
                                 cTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCa                                                                                                       >  6:99743/1‑149 (MQ=255)
                                                  tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAGAAAAAACGAAGGTTCGGT‑‑GCCTGTAAGGGCCTCTTAGC‑‑GGGCGTGCCCGGGAATAGGCCCTGAAAGGCTAGAACCGAAAAATCAACATCAGCGCCTGCATCTgg                                                                                      >  1:228322/1‑149 (MQ=255)
                                                  tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGAGAAGATGCCCTGACAGGCAAGACGCGAAACTGATCCATCAGCGCCTGCAGCTgg                                                                                      >  2:286242/1‑149 (MQ=39)
                                                  tCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGA‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACAGCTCCATCAGCGCCAGCAGCTgg                                                                                      >  2:265278/1‑149 (MQ=255)
                                                    gaagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTGAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCCCctgcttcgtgtt                                                                                    >  6:4892/1‑142 (MQ=255)
                                                     aagaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGTGGGGCAGTGACGGCCACTTAGC‑‑GGGCGGGGCCGCGAAAAGGGCCGGTCAGGCGAGACGCGAAAATGCTCCCAAAGCGCCAGCGggtggt                                                                                     >  1:92591/1‑149 (MQ=255)
                                                       gaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGc                                                                                 <  3:312486/149‑1 (MQ=255)
                                                       gaaAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGc                                                                                 >  7:41254/1‑149 (MQ=255)
                                                                        gATTAGCCATAAGTTGTGCTCAATAATACCTATGGTGCGGT‑‑GCCTGTCAGTGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCtt                                                                <  1:286242/149‑1 (MQ=39)
                                                                          ttAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTgg                                                              <  8:41254/149‑1 (MQ=255)
                                                                             ggCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                  >  1:230559/1‑110 (MQ=255)
                                                                             ggCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAgc                                                                                                  <  2:230559/110‑1 (MQ=255)
                                                                              gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                         >  8:266191/1‑118 (MQ=255)
                                                                              gCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                         <  7:266191/118‑1 (MQ=255)
                                                                                 tAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑TCCTGTCAGGGCCTCTTCGCCACCGCGGGCCCCAGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGAt                                                         <  1:265278/149‑1 (MQ=255)
                                                                                   aGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGt                                                     <  5:99743/149‑1 (MQ=255)
                                                                                          cTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                         <  3:95881/106‑1 (MQ=255)
                                                                                          cTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCctgctgc                                                                                         >  4:95881/1‑106 (MQ=255)
                                                                                          cTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCtt                                                                <  3:80035/131‑1 (MQ=255)
                                                                                          cTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCtt                                                                >  4:80035/1‑131 (MQ=255)
                                                                                           tCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAg                                                                     >  5:199641/1‑125 (MQ=255)
                                                                                           tCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAg                                                                     <  6:199641/125‑1 (MQ=255)
                                                                                               gAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCg                                        <  8:153380/149‑1 (MQ=255)
                                                                                                aaTAACGTTGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCgg                                                                          <  1:145704/115‑1 (MQ=255)
                                                                                                aaTAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCGAGCAGCGAAACTGCACCCACAGCGCCAGCTGCAGGTTGGGCGGCCGGATCAGGGCGTGGCTGATGCGTGCCGAATGATCgg                                       >  2:171412/1‑149 (MQ=25)
                                                                                                  tAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGcc                                     >  6:61746/1‑149 (MQ=255)
                                                                                                    aCGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATGAATCGGcac                                    <  7:248670/148‑3 (MQ=255)
                                                                                                      gATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGcc                                 >  7:213419/1‑149 (MQ=255)
                                                                                                            gcggt‑‑gcCTGTCAGGGCCGCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCccagcca                           <  6:194935/149‑1 (MQ=255)
                                                                                                                   gcCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTCGCCGATTCATCGGCCTGCCCAGCCAGTGa                       <  1:24424/149‑1 (MQ=255)
                                                                                                                    ccTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGAt                      <  5:61746/148‑1 (MQ=255)
                                                                                                                       gTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTg                  >  3:258369/1‑149 (MQ=255)
                                                                                                                           gggCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCAc              >  8:169071/1‑149 (MQ=255)
                                                                                                                               ctctTCGC‑‑GGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGc                                      >  6:2078/1‑121 (MQ=255)
                                                                                                                               ctctTCGC‑‑GGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGc                                      <  5:2078/121‑1 (MQ=255)
                                                                                                                                 ctTCGC‑‑GGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGtctc         >  1:271144/1‑148 (MQ=255)
                                                                                                                                  ttCGC‑‑GGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg        >  2:118148/1‑148 (MQ=255)
                                                                                                                                  ttCGC‑‑GGGCACGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCg        >  8:75907/1‑148 (MQ=255)
                                                                                                                                     gc‑‑gGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATg    >  2:164749/1‑149 (MQ=255)
                                                                                                                                         gggCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGcc  <  2:92591/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
ACAGGGTGGCGGCGATCGAGGTGCGTAAAAGCCCTTGTTCGCTGGTGATGTCGAAGAAAGTGGATCGGGTTGGATTAGGCATAAGTTGTGCTCAAGAATAACGATGGTGCGGT‑‑GCCTGTCAGGGCCTCTTCGC‑‑GGGCGTGCCCGCGAAGAGGCCCTGACAGGCTAGACGCGAAACTGCTCCATCAGCGCCTGCTGCTGGTTGGCCAGCCGGTTCAGGGCTTGGCTGATGCGTGCCGATTCATCGGCCTGCCCAGCCAGTGATTCTGTCACGTCTCGAATGCC  >  NC_002947/1403453‑1403734

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