Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,974,681 C→T 10.8% intergenic (‑209/‑132) PP_3503 ← / → PP_3504 sigma‑54 dependent transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,974,6810CT10.8% 93.9 / 3.3 37intergenic (‑209/‑132)PP_3503/PP_3504sigma‑54 dependent transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (2/2);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

CAGTAGCAACGAAGTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCA  >  NC_002947/3974538‑3974813
                                                                                                                                               |                                                                                                                                     
cAGTAGCAACGAAGTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAAtg                                                                                                                                  >  8:42537/1‑149 (MQ=255)
     gCAACGAAGTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTAAAAACCGCCAAGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAAtgtgtg                                                                                                                              >  1:204099/1‑148 (MQ=255)
      cAACGAAGTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGAc                                                                                                                            >  7:85760/1‑149 (MQ=255)
             gTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATGACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCa                                                                                                                     <  8:85760/149‑1 (MQ=255)
                            tACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATc                                                                                                      <  6:134966/149‑1 (MQ=255)
                                    gtctgCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACa                                                                                             >  7:25166/1‑149 (MQ=255)
                                        gCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAAc                                                                                         >  8:52691/1‑149 (MQ=255)
                                           gAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAAtgtgt                                                                                                                               <  4:241993/109‑1 (MQ=255)
                                           gAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAAtgtgt                                                                                                                               >  3:241993/1‑109 (MQ=255)
                                                        aCAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGAc                                                                                              <  5:146149/128‑1 (MQ=255)
                                                        aCAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGAc                                                                                              >  6:146149/1‑128 (MQ=255)
                                                         cAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTcg                                                                         <  7:52691/148‑1 (MQ=255)
                                                          aaaCCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTcgc                                                                        >  4:258662/1‑148 (MQ=255)
                                                                 cACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATg                                                                <  7:42537/149‑1 (MQ=255)
                                                                   cGGGAATTTTCGCACTTTATCGAAATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACACACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAg                                                              <  1:73642/149‑1 (MQ=255)
                                                                     ggAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGc                                                            >  6:24573/1‑149 (MQ=255)
                                                                      gAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGc                                                            >  8:25090/1‑148 (MQ=255)
                                                                            tCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTTAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCACGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACtt                                                       <  1:191584/148‑1 (MQ=255)
                                                                              gCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGaa                                                   <  8:25166/149‑1 (MQ=255)
                                                                                    tATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAAc                                              <  5:24573/148‑1 (MQ=255)
                                                                                      tCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAAcgcg                                           >  1:183768/1‑149 (MQ=255)
                                                                                      tCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAAcgc                                            <  2:204099/148‑1 (MQ=255)
                                                                                         tcATCACATGGCGGCGCATGAAAGAACAGCGTAGACCCTGAAAAATATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGaa                                                   <  7:223337/137‑1 (MQ=255)
                                                                                         aCATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGac                                        >  6:263099/1‑149 (MQ=255)
                                                                                         aCATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCAGAAAAAAATTTACAGCCAATCAGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCCCTGCACTTCGCCATCAATGAGGCGAAATTGAAgatatagg                                           >  8:223337/1‑139 (MQ=255)
                                                                                              acaTGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGacttca                                    <  7:25090/148‑1 (MQ=255)
                                                                                                  ggcggcGCGGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGacttcacttc                                >  8:232547/1‑148 (MQ=255)
                                                                                                          aGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTgg                       <  2:183768/149‑1 (MQ=255)
                                                                                                            gCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGa                     >  4:171157/1‑149 (MQ=255)
                                                                                                              aaGAACAGCGTCGACCCTGAAAAAAATTTACAGTCAATGTGTGACTTACCCACGaacaac                                                                                                             >  8:87636/1‑60 (MQ=255)
                                                                                                              aaGAACAGCGTCGACCCTGAAAAAAATTTACAGTCAATGTGTGACTTACCCACGaacaac                                                                                                             >  6:171770/1‑60 (MQ=255)
                                                                                                              aaGAACAGCGTCGACCCTGAAAAAAATTTACAGTCAATGTGTGACTTACCCACGaacaac                                                                                                             <  7:87636/60‑1 (MQ=255)
                                                                                                              aaGAACAGCGTCGACCCTGAAAAAAATTTACAGTCAATGTGTGACTTACCCACGaacaac                                                                                                             <  5:171770/60‑1 (MQ=255)
                                                                                                                  aCAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATa                <  5:263099/148‑1 (MQ=255)
                                                                                                                             ccTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAAc    <  3:258662/149‑1 (MQ=255)
                                                                                                                               tGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGATGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACCAGGACTTCACTTCGATGTTTGGGACCATATAGAGACCGAACCa  >  7:43796/1‑149 (MQ=255)
                                                                                                                                       aTTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCCATGTTTgg                       <  8:190226/120‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                     
CAGTAGCAACGAAGTTGACGGTTCCCAGTACCGTCTGTCTGCCGAATCATTTTTTTACAAACCGCCACGGGAATTTTCGCACTTTATCGACATCACATGGCGGCGCAGGCAAGAACAGCGTCGACCCTGAAAAAAATTTACAGCCAATGTGTGACTTACCCACGAACAACGAGCATCAGTAGACATAACCCCGCGCTGCACTTCGCGATCAATGAGGCGACTTTGAAGCAACGCGGACTTCACTTCGATGTTTGGGACCATAGAGAGACCGAACCA  >  NC_002947/3974538‑3974813

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: