Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,104,620 C→T 8.1% R387C (CGC→TGC)  PP_3612 → TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,104,6200CT 8.1% 137.6 / 2.4 49R387C (CGC→TGC) PP_3612TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TGCTGATCAAGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGC  >  NC_002947/4104474‑4104764
                                                                                                                                                  |                                                                                                                                                 
tgctgATCAAGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGAGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAgcgc                                                                                                                                                 >  1:85260/1‑149 (MQ=255)
        aaGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGa                                                                                                                                          <  6:134830/148‑1 (MQ=255)
         aGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGGAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGc                                                                                                                                        <  5:33256/149‑1 (MQ=255)
           cTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGTGCGTGCTGAGc                                                                                                                                        >  8:59183/1‑147 (MQ=255)
           cTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGTGCGTGCTGAGc                                                                                                                                        >  4:17304/1‑147 (MQ=255)
           cTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCa                                                                                                                                      <  6:246202/149‑1 (MQ=255)
            tCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGc                                                                                                                                        >  3:209637/1‑146 (MQ=255)
              aaCCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGAc                                                                                                                                   >  7:216717/1‑149 (MQ=255)
              aaCCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGAc                                                                                                                                   <  4:39601/149‑1 (MQ=255)
                  gcgGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAAc                                                                                                                               <  8:193685/149‑1 (MQ=255)
                   cgGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACg                                                                                                                              >  2:38222/1‑149 (MQ=255)
                       aaTGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGc                                                                                                                          >  7:130438/1‑149 (MQ=255)
                                 gcgATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCg                                                                                                                <  8:216717/149‑1 (MQ=255)
                                        ttACGGTTTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTAcc                                                                                                         <  1:343886/148‑1 (MQ=255)
                                         tACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCg                                                                                                        <  2:85260/149‑1 (MQ=255)
                                         tACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCg                                                                                                        >  5:63921/1‑149 (MQ=255)
                                          aCGGTATGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACAAACGTGGCCCAGGTACCGgtgaatc                                                                                                         >  2:343886/1‑143 (MQ=255)
                                            ggTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGa                                                                                                     >  6:306640/1‑149 (MQ=255)
                                            ggTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGa                                                                                                     >  7:338441/1‑149 (MQ=255)
                                                 gtTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCt                                                                                                >  3:323148/1‑149 (MQ=255)
                                                    cAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGTGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTg                                                                                              <  7:59183/148‑1 (MQ=255)
                                                    cAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGTGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTg                                                                                              <  3:17304/148‑1 (MQ=255)
                                                    cAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTg                                                                                              <  5:213319/148‑1 (MQ=255)
                                                         aGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTCCGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACAAACGTGGCCCAGGTACCGGTGTTCCGAGACGCCTTCGTGCCg                                                                                        >  4:170599/1‑149 (MQ=255)
                                                         aGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGc                                                                                          >  6:40342/1‑147 (MQ=255)
                                                         aGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGc                                                                                          <  5:40342/147‑1 (MQ=255)
                                                             ggAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGt                                                                                                                                               <  7:313917/90‑1 (MQ=255)
                                                             ggAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGt                                                                                                                                               >  8:313917/1‑90 (MQ=255)
                                                              gAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑cc                                                                                  <  6:125665/149‑1 (MQ=255)
                                                               aGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCg                                                                                 >  1:321863/1‑149 (MQ=255)
                                                                     aGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACAt                                                                           <  5:306640/149‑1 (MQ=255)
                                                                         acaaGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATgccg                                                                       <  4:209637/149‑1 (MQ=255)
                                                                             gcgcAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCaa                                                                   <  6:63921/149‑1 (MQ=255)
                                                                              cgcAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAAc                                                                  <  3:170599/149‑1 (MQ=255)
                                                                              cgcAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAAc                                                                  <  8:130438/149‑1 (MQ=255)
                                                                                           cATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAgg                                                     >  2:249238/1‑149 (MQ=255)
                                                                                              gcagcaCAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCa                                                  >  5:228854/1‑149 (MQ=255)
                                                                                               cagcaCAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCCGGAACATGCCGCCAACCTGTCTTCCAAGGGc                                                   <  1:38222/148‑1 (MQ=255)
                                                                                                         tgttgtACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCcag                                       >  7:129028/1‑149 (MQ=255)
                                                                                                           ttgtACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCcagca                                     <  4:323148/149‑1 (MQ=255)
                                                                                                                  gCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGaccacc                              >  1:132234/1‑149 (MQ=255)
                                                                                                                  gCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGaccacc                              <  8:129028/149‑1 (MQ=255)
                                                                                                                     tGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTgt                                                              >  4:135194/1‑114 (MQ=255)
                                                                                                                     tGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTgt                                                              <  3:135194/114‑1 (MQ=255)
                                                                                                                      ggAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACCTGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTc                          >  1:316987/1‑149 (MQ=255)
                                                                                                                             ggCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACGGTGACGCCTTGTTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCGGTGGTCTgctgga                   <  8:262277/149‑6 (MQ=255)
                                                                                                                                tttCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTg                  >  3:199184/1‑147 (MQ=255)
                                                                                                                                            tACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATc      >  4:104647/1‑147 (MQ=255)
                                                                                                                                             aCCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAg   >  3:117732/1‑149 (MQ=255)
                                                                                                                                              ccAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGc  >  6:267654/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGCTGATCAAGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTA‑CCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGC  >  NC_002947/4104474‑4104764

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