Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,164 G→T 22.4% intergenic (‑128/‑140) PP_5175 ← / → PP_5176 HlyD family secretion protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,1640GT22.4% 43.5 / 11.8 31intergenic (‑128/‑140)PP_5175/PP_5176HlyD family secretion protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base T (3/4);  total (12/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 = 5.16e-01

CGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATG  >  NC_002947/5903019‑5903306
                                                                                                                                                 |                                                                                                                                              
cGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGc                                                                                                                                             >  3:20431/1‑149 (MQ=255)
     gCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcgcg                                                                                                                                        <  4:20431/149‑1 (MQ=255)
      caccAACGCCAGGCCGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGGTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGa                                                                                                                                       <  8:248506/146‑1 (MQ=255)
       aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGaa                                                                                                                                      <  3:169770/149‑1 (MQ=255)
       aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCcgcg                                                                                                                                        >  2:217884/1‑147 (MQ=255)
       aTCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGGTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGaa                                                                                                                                      <  4:83418/149‑1 (MQ=255)
          aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAaga                                                                                                                                    <  5:313389/148‑1 (MQ=255)
          aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAagag                                                                                                                                   <  1:217884/149‑1 (MQ=255)
           aCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAgg                                                                                                                                  >  7:11480/1‑149 (MQ=255)
            cGCCATGCTGGCTCCTTGCTTGTGCACAAATGTTTTGTGGGATCCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGc                                                                                                                                 <  7:69452/149‑1 (MQ=255)
                          cttgcttgTGCACTAATGCTTTTTGTGCTGCTACTGTGAGCGTAGTGCTGGGTTCGATATACATCACTGAGGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCGTCTTCGGGGGCGAGCCCGCGAAGAGGCCTCCACAGCCAt                                                                                                                     <  6:334753/147‑1 (MQ=255)
                           ttgcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTa                                                                                                                  >  1:243624/1‑149 (MQ=255)
                               ttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAg                                                                                                              >  8:61689/1‑149 (MQ=255)
                               ccgatctcaAATGATTTGGGGCATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGCCGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTGCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa                                                                                                               <  5:205259/141‑1 (MQ=255)
                                   gCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCACATAGCGGGCCACCCCGCGAAGACGCCTCCACACCCATAGATGAGTAt                                                                                                           >  4:111759/1‑148 (MQ=255)
                                       agaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGGGTGCCCGCGCCGAGGGCGCCCCAGCCATGTATAAGTAtgtggc                                                                                                      >  4:192243/3‑147 (MQ=255)
                                               tgtgGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa                                                                                               <  6:121487/148‑1 (MQ=255)
                                                   ggATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGaaa                                                                                           <  7:61689/148‑1 (MQ=255)
                                                      tGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCAt                                                                                       >  1:42312/1‑149 (MQ=255)
                                                        cGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCa                                                                                     <  2:243624/149‑1 (MQ=255)
                                                                      gTTCAGGGTTCGATAATCTTCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATc                                                                        <  3:111759/148‑1 (MQ=255)
                                                                                      aCATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTAc                                                       <  2:42312/149‑1 (MQ=255)
                                                                                        aTCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCGTCTGCGCGGGCTAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCt                                                     <  3:162086/149‑1 (MQ=255)
                                                                                                  ggggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGa                                           <  8:11480/149‑1 (MQ=255)
                                                                                                               cATTTTTGGTGCCGGCGGGGGCCTCTGCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCaa                              <  3:192243/149‑1 (MQ=255)
                                                                                                                                gggtCCTCTTCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCg             <  3:186058/149‑1 (MQ=255)
                                                                                                                                 aggCCTCTTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGt            >  1:221023/2‑149 (MQ=255)
                                                                                                                                       ctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc       >  5:43446/1‑148 (MQ=255)
                                                                                                                                        ttCGCGGGCTCGCCCGCGAAGAGGCCCCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCg      >  2:71389/1‑148 (MQ=255)
                                                                                                                                           gcgGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCATGCATCTGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAAt   <  7:236068/148‑1 (MQ=255)
                                                                                                                                           gcgGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATg  >  8:295853/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAATG  >  NC_002947/5903019‑5903306

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: