Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,103,353 T→A 10.6% intergenic (‑108/‑109) PP_3611 ← / → PP_3612 hypothetical protein/TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,103,3530TA10.6% 123.9 / 5.9 47intergenic (‑108/‑109)PP_3611/PP_3612hypothetical protein/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base A (3/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCATTTCCATGC  >  NC_002947/4103206‑4103501
                                                                                                                                                   |                                                                                                                                                    
cAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCAtt                                                                                                                                                     <  3:259811/149‑1 (MQ=255)
         cAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGAAATGCGTCGcca                                                                                                                                            <  2:194503/149‑2 (MQ=255)
                 gcgagAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACt                                                                                                                                    <  3:329154/149‑1 (MQ=255)
                    agagTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCACCGAAATCACt                                                                                                                                    <  7:328798/146‑1 (MQ=255)
                    agagTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGACTTGCCGGCCCTCGCGGCTCGCGACGCATTTCGTCACCGAAATCACt                                                                                                                                    >  8:328798/1‑146 (MQ=255)
                       gTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTtata                                                                                                                              <  6:145145/149‑1 (MQ=255)
                                  gTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGc                                                                                                                            >  5:228304/1‑140 (MQ=255)
                                  gTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGc                                                                                                                            <  6:228304/140‑1 (MQ=255)
                                     cATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAAtt                                                                                                                 <  3:195236/148‑1 (MQ=255)
                                        aGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCtt                                                                                                              >  7:239157/1‑148 (MQ=255)
                                         ggCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTa                                                                                                            <  8:239157/149‑1 (MQ=255)
                                                  cccTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGc                                                                                                     <  6:147898/147‑1 (MQ=255)
                                                       tgATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTAc                                                                                                <  1:299157/147‑1 (MQ=255)
                                                       tgATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACa                                                                                               <  5:83413/148‑1 (MQ=255)
                                                         aTCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACAttt                                                                                            <  6:51213/149‑1 (MQ=255)
                                                          tCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTa                                                                                                            >  6:278506/1‑132 (MQ=255)
                                                          tCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTa                                                                                                            <  5:278506/132‑1 (MQ=255)
                                                             aGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCa                                                                                         >  1:299039/1‑148 (MQ=255)
                                                               tCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTtc                                                                                                                                                   >  6:210904/1‑88 (MQ=255)
                                                               tCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTtc                                                                                                                                                   <  5:210904/88‑1 (MQ=255)
                                                                   aaTAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCAt                                                                                   >  3:186128/1‑148 (MQ=255)
                                                                     tAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTAc                                                                                <  2:299039/149‑1 (MQ=255)
                                                                       ggCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTa                                                                              >  3:166067/1‑149 (MQ=255)
                                                                          gcgTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTc                                                                           >  8:110982/1‑149 (MQ=255)
                                                                            gTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGt                                                                         >  3:43789/1‑149 (MQ=255)
                                                                                  ggCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTTCCGGCACTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTTCACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTcc                                                                   <  7:173610/149‑1 (MQ=255)
                                                                                   gCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTa                                                                                 <  5:145471/134‑1 (MQ=255)
                                                                                   gCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTa                                                                                 >  6:145471/1‑134 (MQ=255)
                                                                                       aCGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCAATTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGa                                                                <  5:151402/147‑1 (MQ=255)
                                                                                       aCGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCAATTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGa                                                                >  6:151402/1‑147 (MQ=255)
                                                                                         gcgcTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTGTAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTc                                                                           <  7:221425/134‑1 (MQ=255)
                                                                                         gcgcTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTGTAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTc                                                                           >  8:221425/1‑134 (MQ=255)
                                                                                           gcTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATaacaa                                                          >  1:337774/1‑149 (MQ=255)
                                                                                            cTGACAACCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATaacaac                                                         <  2:70838/149‑1 (MQ=255)
                                                                                               aCAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTCTTATAGCTAATACGAATTCTTTACACTTCCTTTACATTTGCAACTCATTACTATTCGTGACTCCCGATAACAACaaa                                                      <  4:185475/149‑1 (MQ=255)
                                                                                                    cATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTc                                                 <  2:337774/149‑1 (MQ=255)
                                                                                                          cATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCttt                                                                                                  >  5:88844/1‑94 (MQ=255)
                                                                                                          cATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCttt                                                                                                  <  6:88844/94‑1 (MQ=255)
                                                                                                           aTTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGt                                           >  7:98357/1‑148 (MQ=255)
                                                                                                                cGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGc                                     <  4:43789/149‑1 (MQ=255)
                                                                                                                  gccgccTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCa                                    >  2:168414/1‑148 (MQ=255)
                                                                                                                       cTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAgg                              >  4:55074/1‑149 (MQ=255)
                                                                                                                                ccTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGcc                     >  3:224129/1‑149 (MQ=255)
                                                                                                                                ccTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGcc                     <  4:166067/149‑1 (MQ=255)
                                                                                                                                         cggCGACGAAATGCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATc             >  1:170351/2‑148 (MQ=255)
                                                                                                                                          ggCGACGAAATGCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCa           >  1:104021/1‑149 (MQ=255)
                                                                                                                                           gCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCAt          >  7:205854/1‑149 (MQ=255)
                                                                                                                                                 cATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTATTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCATTTCCAt    <  5:66418/149‑1 (MQ=255)
                                                                                                                                                   tttCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGAATCCCGATAACAACAAAGAGTCGACTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCATTTCCATGc  >  5:136707/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAAGGGTTACAGCACGAGCGAGAGTACGAATGGGGTTCATAGGCGTGCTTCCCTGTGATCAAGTCCGAATAGGCGCGTGAAAGGCTCACGCGCTGACAAGCATAGCCATTACCGGCCGCCTTGCCGGCCCTCTCGGCTGGCGACGCATTTCGTCGCCGAAATCACTGTTATAGCTAATACGAATTCTTTACACTTGCTTTACATTTGCAAATCATTACTATTCGTGACTCCCGATAACAACAAAGAGTCGCCTGTAATGGCAAAAAGGTCTGTTGCCCCAGCATCCATTTCCATGC  >  NC_002947/4103206‑4103501

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