Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,217,528 T→C 23.5% intergenic (+128/+26) yqjI → / ← aer PadR family putative transcriptional regulator/fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,217,5280TC23.5% 66.0 / 2.7 37intergenic (+128/+26)yqjI/aerPadR family putative transcriptional regulator/fused signal transducer for aerotaxis sensory component/methyl accepting chemotaxis component
Reads supporting (aligned to +/- strand):  ref base T (23/3);  new base C (2/7);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 5.26e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCAGTGATGCACAAATAAAAAAGATCATTGCGGTGATCGACCGCGCCGCTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTCCAGCGTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGT  >  NC_000913/3217324‑3217753
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
atatCAGTGATGCACAAATAAAAAAGATCATTGCGGTGATCGACCGCGCCGCTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCACCCCCCGGATCCGGCATAcc                                                                                                                                                                                                            >  2:187792/1‑228 (MQ=255)
         ttgcaaacaaaaaaaaGATCAGTGGGCACAGCGCCCGGCCCGATTTTGATACTCCGCAAGGAAAGTAGTCGCCGCACCCGCCAGTGGCGGGGTGCAATGCCCGGATCGGCCGCTGTACCACCCTTACCCGGCCTACCCCCGCTACACCCCCCGCAGGCCTGATAGTACGGGCCAGCATCGAACCGGGAATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGc                                                                                                                                                                                                   <  1:469846/218‑1 (MQ=255)
           gCACAAATAAAAAAGATCATTGCGGTGATCGACCGCGCCGCTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGCTAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCAGCATACCGATTAATGCa                                                                                                                                                                                                  >  2:426778/1‑227 (MQ=255)
                                      aTCGACCGCGCCGCTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGGACCGTCACCGCGGCTTCCCGTCGCCt                                                                                                                                                                      >  2:513965/1‑228 (MQ=255)
                                              cgccgcTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGCCGCGCCTTATCCGGCCTCCACCCGCTCCACACCCCGCCGGCCCGATAAGTTGCGCCACCCTCTCACCAGGCCATGTCCCACTACCGCCGGACTCGGCATCCCGATTGATGCGGTCCCGCCCCCCCGCCCTCCACCCGGCCTGcgcgc                                                                                                                                                               >  1:195675/1‑226 (MQ=19)
                                                              aCGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCCCGGTGTTTCACCATCGCCg                                                                                                                                               >  2:28875/1‑227 (MQ=255)
                                                              aCGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGGGCTCCAACCGCCGGATCCGGCAACCCGATTAATGCAGTACCTTCACCGCGCCTTCCGGTCGCCTGCCCCGGGGTTTCACCACTGCCg                                                                                                                                               >  2:473530/1‑227 (MQ=38)
                                                                 cAAGTGGATTAACCCCCAGCCCCCGCAGGTGGCCGGTGCAATTGGCCGGATCGACCTGTGGCGGCGCCTTGACCGGCTCCACCCGATTACCCCACCCGCGGGCCTTATAAGATGCGACAGCAGCTAAGCAGGCGTTGGGCACCAACCTCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGaca                                                                                                                                            <  1:168687/228‑1 (MQ=255)
                                                                       gATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGGTTCACCATCGCCGACCCCTgcg                                                                                                                                      >  2:14252/1‑227 (MQ=255)
                                                                              cgccgcATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGCTTAATGCAGGCCCGTCCCCCCGTCTTCCCGTCGGCTGCCGCGGGGTTTCACCACACCCGACCCCCgcgccccccct                                                                                                                              >  1:340499/1‑221 (MQ=38)
                                                                                 cgcATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCCCGTCTTCCAGTCGCCTGGCGCGGGGTTTCACCATACCCGACACCTTCGCACCCCCTTc                                                                                                                            >  2:280533/1‑227 (MQ=39)
                                                                                              tGGCGCGGTGCAATTCCCGGATGCGCCGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTCCTCCAACCGCCGGCTCCGGCATCCCGATTAATGCAGTACCGTCACCGCGTCGTCCAGTCGGCTGCCCCGCTGTTTCACCATCGCCCACACCTGCGCCCTCCCCTCCCCCACCCCGGCa                                                                                                               >  1:164486/1‑228 (MQ=255)
                                                                                                               ccGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCCCCAACCGCCGGATCCGGCATACCGa                                                                                                                                                                                                          >  2:205805/1‑121 (MQ=255)
                                                                                                                cGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCCCCAACCGCCGGATCCGGCATACCGAt                                                                                                                                                                                                         <  1:205805/121‑1 (MQ=255)
                                                                                                                 ggATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTCTCCTCCAACCGCCGGAACCGGCACACCGCTTAATGCAGCACCGTCACCGCTTCTTCCAGTCGGCTGGCGCGGCGGTTCACCCTCGCCGCCACCTGCGCACACTCTCCCACCAGCCCGGCGTTTTTCTGGGCGACCAg                                                                                              >  1:536138/1‑227 (MQ=255)
                                                                                                                                    gcCTTACCCGGCCTACCCCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCCTCGAATGAGGCATTGGGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGc                                                                            <  2:409489/227‑1 (MQ=255)
                                                                                                                                     ccTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCAc                                                                          >  2:487449/1‑228 (MQ=255)
                                                                                                                                               ccTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCCTCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCCCGTCTTCCAGTCGGCTGGGCCGGTGTTTCACCATCTCCGCCACCTGCCCACTCCCTCCCACCCACTCGCCATTGTTCCGGTGGATCAGGTTAAGCTCACCCCCCCGACGGGTccccct                                                                >  1:516467/1‑223 (MQ=255)
                                                                                                                                                  acccccGCTACACTCCCCGCAGGCCTGATAAGAGGCGCAAGACTCGCATCAGGCTTTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGa                                                             <  2:232546/225‑1 (MQ=255)
                                                                                                                                                           cacacccCCCACGTGCCTGAAAAGCTGCCCCAGCATCGCATCAGGCATTGTGCCCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCAt                                                     <  2:536138/227‑1 (MQ=255)
                                                                                                                                                             ccgcccccccaaGGCCTGTAAGAAGGGCCCAGCATCGCATCAGGAATTGTGCCCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCCCCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCCCCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCAt                                                     <  2:555866/216‑1 (MQ=255)
                                                                                                                                                                 ccccGCAAGGCCGTATAAAATGCCCCAGAACCGCATAAGGCATGTGCCACCATCCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGcc                                                <  2:340499/228‑1 (MQ=255)
                                                                                                                                                                 ccccGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGGGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGGGATCAGGTTAAGCTCCTCCACTGCCCGGGTCAGGCCGGAAAGCCCATCGGCCt                                               >  2:573284/1‑228 (MQ=255)
                                                                                                                                                                     gCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGCAACCTGCGCCACCTCTTCCACCCACTCGGCATTTTTCTGGGGGATCCGGTTAAGCCCATCCACTCCCCGGTTCAGCCTGTAAACCCCAGCGGCCGGTTc                                           >  1:488414/1‑228 (MQ=255)
                                                                                                                                                                     gCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCCCTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTTTTc                                           >  2:207271/1‑228 (MQ=255)
                                                                                                                                                                     gCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTCCACGGGTCAGGCTGGAAAGCCCATCGGCCCGTTc                                           >  2:390427/1‑228 (MQ=255)
                                                                                                                                                                     gCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTATCTGGGTGATCAGGTTCAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTCTTc                                           >  2:294720/1‑228 (MQ=255)
                                                                                                                                                                     gCAGGCCGGATGAGATGCGCCAGCCGCGCATCGGGCTTTGGTCCCCACCCGCCGGATCCGGCATACCGATTAATGCAGTTCCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTc                                           <  2:516467/228‑1 (MQ=255)
                                                                                                                                                                        aGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGCTTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGGTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGGGATCAGGTTAAGCTCATCCACTGCACGGGGCCGGCGGGGAACGCACTCGGCCTGTTccg                                         >  1:373752/1‑227 (MQ=255)
                                                                                                                                                                             tGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTTGGTGATCAGGTTAAGCTCATCCACTCCACGGGTCAGGCTGGCAAGCCCATCTGCCTGTTCCACCGt                                     >  2:664953/1‑227 (MQ=255)
                                                                                                                                                                                    aTGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGGTAACCTCCTCCACTGCACGGGTACGGATGGGAAGCCCATCGCCCTGTCCCCGCGTTGCATGGc                             >  1:587988/1‑228 (MQ=255)
                                                                                                                                                                                    aTGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGACTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGCCATTTTTCTGGGTGTTCAGGTTAAGCTCATCCACTGCACGGTCCAGCCGGGCAAGCCCATCGGCCTGTTCCACCGTTGAATgga                             >  2:91546/1‑227 (MQ=255)
                                                                                                                                                                                           cagctcccgctcAGGCTTTGTGCCCCAACCGCCGGATCCGGCATTCCGATTAATGCAGTACCGGCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTCCAGCGTTGAATGGCTAATCt                       <  2:373752/218‑1 (MQ=255)
                                                                                                                                                                                                             tgtgCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCCCCGCCACCTGCGCAACCTCTTCACCCCGCTCGGCATTTTTCTGGGTGCTCAGGTTAAGCCCATCCCCTGACCCGGTCAGGCTGGACAGCCACTCGGCCGTTTCCAACGTTGAATGGCTCATCTGGGCGGTCAACCggggg      >  1:438818/1‑225 (MQ=255)
                                                                                                                                                                                                             tgtgCCCCACCCGCCGGATCCGGCGTACGGTTTAATGCCGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTCCAGCGTTGAATGGCTAATCTGGGCGATCAACTGGGTGa     <  1:22323/228‑1 (MQ=255)
                                                                                                                                                                                                              gtgCCCCACCCGCCGGATCCGGCATGCCGAGTAAGGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCAATCTCTTCCTCCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTCCATCGTTGAATGGCTAATCTGGGCGATCAACTGGGTGa     <  2:438818/227‑1 (MQ=255)
                                                                                                                                                                                                                  cTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGCCTGGCAAGCCCATCGGCCTGTTCCAGCGTTGAATGGATAATCTGGGCGACCACCTGGGTGACGt  >  2:613835/1‑227 (MQ=255)
                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 
ATATCAGTGATGCACAAATAAAAAAGATCATTGCGGTGATCGACCGCGCCGCTTTTGATATTACGCAACTGGATTAATCGCCGCATCCGCCAGTGGCGCGGTGCAATTGCCGGATGCGACGCTTGACGCGCCTTATCCGGCCTACACCCGCTACACACCCCGCAGGCCTGATAAGATGCGCCAGCATCGCATCAGGCATTGTGCTCCAACCGCCGGATCCGGCATACCGATTAATGCAGTACCGTCACCGCGTCTTCCAGTCGGCTGGCGCGGTGTTTCACCATCGCCGACACCTGCGCACTCTCTTCCACCAGCTCGGCATTTTTCTGGGTGATCAGGTTAAGCTCATCCACTGCACGGGTCAGGCTGGAAAGCCCATCGGCCTGTTCCAGCGTTGAATGGCTAATCTGGGCGATCAACTGGGTGACGT  >  NC_000913/3217324‑3217753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: