Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,936,919 A→T 14.2% intergenic (‑330/‑88) flhA ← / → PP_4345 flagellar export pore protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,936,9190AT14.2% 102.0 / 7.6 42intergenic (‑330/‑88)flhA/PP_4345flagellar export pore protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/19);  new base T (2/4);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

TCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGA  >  NC_002947/4936772‑4937064
                                                                                                                                                   |                                                                                                                                                 
tctACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAca                                                                                                                                                   >  3:149418/1‑148 (MQ=255)
    cGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatc                                                                                                                                                 <  4:343685/146‑1 (MQ=255)
     ggCAATGGCCAAAGGCTTAAGCTCTTTGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGc                                                                                                                                             >  8:136532/1‑149 (MQ=255)
         aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                         <  8:278293/149‑4 (MQ=255)
         aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                         <  1:69409/149‑4 (MQ=255)
         aTGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGATGTTCATGccttt                                                                                                                                         <  5:99754/149‑4 (MQ=255)
                         gCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCTAACTGCTCGATGGTAGctc                                                                                                                         <  6:89343/149‑1 (MQ=255)
                           tctTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatca                                                                                                                                                >  6:117281/1‑124 (MQ=255)
                           tctTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAAcatca                                                                                                                                                <  5:117281/124‑1 (MQ=255)
                             ttGGAAAGCAACAGAGAGGTTGGCGCAACGACAAAAAAAATCTGAGCTACTTGATCATGCAAAAGCTCACGGTTTCTATGCGTATGTAGAGAGCCAATTACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAg                                                                                                                      <  1:73858/148‑1 (MQ=255)
                               ggAAAGCCACAGAGATGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTTATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                     <  7:243833/147‑1 (MQ=255)
                               ggAAAGCCACAGAGAGGGTGGCGCCACGTCAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGt                                                                                                                     <  5:293307/147‑1 (MQ=255)
                                  aaGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTc                                                                                                                >  1:194757/1‑149 (MQ=255)
                                           agagGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATAcc                                                                                                       >  4:109533/1‑149 (MQ=255)
                                               ggTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCg                                                                                                   >  8:121465/1‑149 (MQ=255)
                                                       cACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCt                                                                                           >  2:322129/1‑149 (MQ=255)
                                                          gACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGATCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTTATACCTTCGCCAGCGCtttt                                                                                        <  2:74631/149‑1 (MQ=255)
                                                            cAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGa                                                                                      <  4:149418/149‑1 (MQ=255)
                                                             aaaaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGAc                                                                                     <  5:132430/149‑1 (MQ=255)
                                                             aaaaaaaaTCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGAc                                                                                     <  1:322129/149‑1 (MQ=255)
                                                                                 gATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGgag                                                                  <  7:136532/148‑1 (MQ=255)
                                                                                      tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAAc                                                               <  6:44818/146‑1 (MQ=255)
                                                                                      tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAAc                                                               >  5:44818/1‑146 (MQ=255)
                                                                                      tGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACAtg                                                            <  3:109533/149‑1 (MQ=255)
                                                                                           aaGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGc                                                       <  2:194757/149‑1 (MQ=255)
                                                                                            aGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCg                                                      <  5:111115/149‑1 (MQ=255)
                                                                                               tCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATCCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGATATTGCCTACATCGCTAGACCCGGAGAACATGGTGGCGAAc                                                   >  8:120621/1‑149 (MQ=255)
                                                                                                    gTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctct                                              <  7:121465/149‑1 (MQ=255)
                                                                                                    gTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAActctct                                              <  7:21395/149‑1 (MQ=255)
                                                                                                         tATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGa                                           >  3:101820/1‑147 (MQ=255)
                                                                                                                          ccAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGAc                        <  4:101820/149‑1 (MQ=255)
                                                                                                                           cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTa                                                                                   >  5:109350/1‑89 (MQ=255)
                                                                                                                           cAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTa                                                                                   <  6:109350/89‑1 (MQ=255)
                                                                                                                                    ctCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGc               >  5:18228/1‑148 (MQ=255)
                                                                                                                                     tCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCa             >  2:144876/1‑149 (MQ=255)
                                                                                                                                     tCAAAGGCATGAACATCATGCCAACTACTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAAc                                                               >  4:212370/1‑99 (MQ=255)
                                                                                                                                     tCAAAGGCATGAACATCATGCCAACTACTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAAc                                                               <  3:212370/99‑1 (MQ=255)
                                                                                                                                       aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGcc              <  5:168641/146‑1 (MQ=255)
                                                                                                                                       aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGcc              >  6:168641/1‑146 (MQ=255)
                                                                                                                                       aaaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCAt            >  5:88744/1‑148 (MQ=255)
                                                                                                                                        ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACtctctc                                             >  1:129292/3‑114 (MQ=255)
                                                                                                                                        ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACtctctc                                             <  2:129292/112‑1 (MQ=255)
                                                                                                                                        ttggCATGATGTTCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa           >  3:314252/3‑148 (MQ=255)
                                                                                                                                        aaGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATa           >  5:282655/1‑148 (MQ=255)
                                                                                                                                             aTGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACGTCGCCATAAAACcg     >  3:46361/1‑149 (MQ=255)
                                                                                                                                                aaCATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGa  >  5:100661/1‑149 (MQ=255)
                                                                                                                                                  catcatGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACaaa                     >  2:2633/1‑128 (MQ=255)
                                                                                                                                                  catcatGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACaaa                     <  1:2633/128‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTACGGCAATGGCCAAAGGCTTAAGCTCTTGGAAAGCCACAGAGAGGGTGGCGCCACGACAAAAAAAATCTGAGCGACTTGATCATGCAAAAGCTCACGGTTTCTATGCGGATGTAGAGAGCCAATGACTACTCAAAGGCATGAACATCATGCCAACTGCTCGATGGTAGCTCTAGTGAGTCTGTGATACCTTCGCCAGCGCTTTTGACTACATCGCTAGACCCGGAGAACATGGTGGCGAACTCTCTCGATAGCTTTATCGTTCCTGACAAACATCGCCATAAAACCGCGA  >  NC_002947/4936772‑4937064

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: