Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,014,552 A→T 8.1% intergenic (+29/‑93) PP_3542 → / → PP_3543 hypothetical protein/iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,014,5520AT 8.1% 141.1 / 2.8 49intergenic (+29/‑93)PP_3542/PP_3543hypothetical protein/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base T (2/2);  total (24/25)
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.96e-01

CTGGCCGGCCAACGCGTGCCCATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCA  >  NC_002947/4014410‑4014678
                                                                                                                                              |                                                                                                                              
cTGGCCGGCCAACGCGTGCCCATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGt                                                                                                                          <  3:12301/149‑1 (MQ=255)
       gccAACGCGTGCCCATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTacc                                                                                                                   <  2:259529/149‑3 (MQ=255)
               gTGCCCATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATc                                                                                                           >  6:152153/1‑149 (MQ=255)
                   ccATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGaa                                                                                                       >  3:54629/1‑149 (MQ=255)
                    cATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGaa                                                                                                       <  6:95081/148‑1 (MQ=255)
                    cATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGaa                                                                                                       >  5:95081/1‑148 (MQ=255)
                    cATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGaa                                                                                                       >  3:53047/1‑148 (MQ=255)
                     aTTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAAcc                                                                                                     >  1:301044/1‑149 (MQ=255)
                     aTTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAAcc                                                                                                     >  8:11083/1‑149 (MQ=255)
                          cGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGc                                                                                                <  4:89661/149‑1 (MQ=255)
                           gAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAAc                                                                                                      >  6:28829/1‑142 (MQ=255)
                           gAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAAc                                                                                                      <  5:28829/142‑1 (MQ=255)
                             aaCCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGa                                                                                             <  2:301044/149‑1 (MQ=255)
                                    gCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGc                                                                                      <  5:152153/149‑1 (MQ=255)
                                    gCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCGGTGCTGATAACAGc                                                                                      >  7:23759/1‑149 (MQ=255)
                                      cGATGCACTGGCAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCtt                                                                                    >  2:266381/1‑149 (MQ=255)
                                                aCAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATAAGGTTTTTATAGTAAAATCTGAACATGTGCTGATAACAGCTACCGCGCCa                                                                            >  3:177345/1‑147 (MQ=255)
                                                  aGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGt                                                                        >  5:50212/1‑149 (MQ=255)
                                                          gTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGc                                                                <  1:266381/149‑1 (MQ=255)
                                                             cGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTTCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCt                                                              <  7:171572/148‑1 (MQ=255)
                                                              ggTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATATCAGCTTCAGCGCCATAGTGCTCATGCCTGa                                                            >  6:108504/1‑149 (MQ=255)
                                                                    tAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTAGGTGATATGAGGCCTGCCGGGGCTATGAGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAgg                                                      >  1:205700/1‑149 (MQ=255)
                                                                      gCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGc                                                     <  4:54629/148‑1 (MQ=255)
                                                                        cgaaGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATTCGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCAATAGTGCTCATGCCTGATCGAGGCCtg                                                  <  4:177345/146‑1 (MQ=255)
                                                                        cGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGttt                                                                                                                        >  8:31826/1‑79 (MQ=255)
                                                                        cGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGttt                                                                                                                        <  7:31826/79‑1 (MQ=255)
                                                                             gCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCt                                                   <  2:64321/143‑1 (MQ=255)
                                                                             gCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCt                                                   >  1:64321/1‑143 (MQ=255)
                                                                             gCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCt                                                   >  5:63920/1‑143 (MQ=255)
                                                                             gCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCt                                                   <  6:63920/143‑1 (MQ=255)
                                                                                    aaGGCCATCGCGCCAGGCTCGGTTATATGAGGCCTTCCGTGGCTATGCTGCAAATCTGATATTGTTTTTATCGAAAAATCTGAACCTGTGCTTATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCAc                                          <  2:206271/145‑1 (MQ=255)
                                                                                    aaGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCAc                                          <  2:304363/145‑1 (MQ=255)
                                                                                    aaGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCAc                                          <  2:299590/145‑1 (MQ=255)
                                                                                           gcgcgcCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGaa                               <  4:53047/149‑1 (MQ=255)
                                                                                                ccAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCtgat                                                >  4:293416/1‑127 (MQ=255)
                                                                                                ccAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCtgat                                                <  3:293416/127‑1 (MQ=255)
                                                                                                ccAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATa                          <  7:11083/149‑1 (MQ=255)
                                                                                                        ggTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCGGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATcc                    <  8:23759/147‑1 (MQ=255)
                                                                                                             tatGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTAc             <  6:50212/149‑1 (MQ=255)
                                                                                                                gAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACa            >  2:42260/1‑147 (MQ=255)
                                                                                                                     cTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCg     >  7:116971/1‑149 (MQ=255)
                                                                                                                      tGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCATGATCGAGGCCTGATGAGCACCGACCCATGAACGATAGAATACCTTTTACAGAGATCga    >  8:116017/1‑148 (MQ=255)
                                                                                                                        ccGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCa  >  8:119974/1‑149 (MQ=255)
                                                                                                                            ggCTATGCGGCAAATCTGTTATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCtt                 >  2:44037/1‑130 (MQ=255)
                                                                                                                            ggCTATGCGGCAAATCTGTTATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCtt                 >  2:48396/1‑130 (MQ=255)
                                                                                                                            ggCTATGCGGCAAATCTGTTATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCtt                 <  1:48396/130‑1 (MQ=255)
                                                                                                                            ggCTATGCGGCAAATCTGTTATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCtt                 <  1:44037/130‑1 (MQ=255)
                                                                                                                                          tCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCtgatg                                               >  7:133006/1‑86 (MQ=255)
                                                                                                                                          tCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCtgatg                                               <  8:133006/86‑1 (MQ=255)
                                                                                                                                              |                                                                                                                              
CTGGCCGGCCAACGCGTGCCCATTACCGAAACCGAGGCCGATGCACTGACAGTGGCGGGTGCGGTGGATAGCCGCAGGCATCTGAAGGCCAGCGCGCCAGGCTCGGTGATATGAGGCCTGCCGGGGCTATGCGGCAAATCTGATATGGTTTTTATCGTAAAATCTGAACCTGTGCTGATAACAGCTTCAGCGCCATAGTGCTCATGCCTGATCGAGGCCTGATGAGCACCGAACCATGAACGATAGAATCCCTTTTACAGAGATCGCCA  >  NC_002947/4014410‑4014678

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: ‑