Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,104,627 T→C 8.7% L389P (CTG→CCG)  PP_3612 → TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,104,6270TC 8.7% 123.3 / 2.8 46L389P (CTG→CCG) PP_3612TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (24/18);  new base C (2/2);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGCTGAAC  >  NC_002947/4104483‑4104769
                                                                                                                                                |                                                                                                                                              
aGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGc                                                                                                                                            <  3:207845/149‑1 (MQ=255)
 gCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACGAGCGCGTGCTGAGc                                                                                                                                            >  4:71807/1‑148 (MQ=255)
         gcgGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAAc                                                                                                                                   >  2:249941/1‑149 (MQ=255)
         gcgGCAATGTCGCCCGAGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAATCAACAAGCGCAACTGGCCGGCATGCAGCTCAACCTGTTGTACGGATTGGAAGCAGGCTTTCAGGACAAGTACCAGAGCGTGATGAGCAAGACCAAc                                                                                                                                   >  6:309036/1‑149 (MQ=255)
         gcgGCAATCACGCCCGCCAAGAGTACGGAGTCTTCAACCAGACGCAGCTCAAGCAACAAGCGCAACTCGCCCGCATGCAGCACAACCAGATGTACGCCGAGGAAGCAGGCTTTCAGGTCAAGTACCAGCGCCTCCTGCGCCAGACCAAc                                                                                                                                   >  6:313527/1‑149 (MQ=255)
            gCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTg                                                                                                                                >  1:223524/1‑149 (MQ=255)
                      ccGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTAc                                                                                                                      >  8:23862/1‑149 (MQ=255)
                               gTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTAcc                                                                                                             >  2:248804/1‑149 (MQ=255)
                                                    ggAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTAc                                                                                        >  4:35112/1‑149 (MQ=255)
                                                           aaGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACTGGTGTACCGTGACGCCTTGGTGCCGGTACCGGaa                                                                                   <  1:249941/147‑1 (MQ=255)
                                                           aaGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCCGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACa                                                                                 >  6:178425/1‑149 (MQ=255)
                                                           aaGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCCGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACa                                                                                 >  6:174063/1‑149 (MQ=255)
                                                             gCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGCCGCCTTGGTGCCGGTACCGGAACATg                                                                               >  8:136292/1‑149 (MQ=255)
                                                              caacaaGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATgc                                                                              <  3:35112/149‑1 (MQ=255)
                                                                acaaGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACTAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATgccg                                                                            <  1:248804/149‑1 (MQ=255)
                                                                acaaGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCg                                                                                            >  1:196454/1‑133 (MQ=255)
                                                                acaaGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCg                                                                                            <  2:196454/133‑1 (MQ=255)
                                                                   aGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATgccgcc                                                                          >  2:183062/1‑148 (MQ=255)
                                                                    gcgcAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCaa                                                                        <  2:223524/149‑1 (MQ=255)
                                                                       cAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATgc                                                                              <  8:230308/140‑1 (MQ=255)
                                                                       cAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATgc                                                                              >  7:230308/1‑140 (MQ=255)
                                                                          cTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCg                                                                                            <  3:5029/123‑1 (MQ=255)
                                                                          cTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCg                                                                                            >  4:5029/1‑123 (MQ=255)
                                                                              ccGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTcc                                                              <  6:70986/149‑1 (MQ=255)
                                                                                     gcagcaCAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCaa                                                            >  5:98526/1‑144 (MQ=255)
                                                                                     gcagcaCAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCaa                                                            <  6:98526/144‑1 (MQ=255)
                                                                                      cagcaCAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCa                                                       <  8:121499/148‑1 (MQ=255)
                                                                                         cacaACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTAcc                                                                                                             <  1:201024/91‑1 (MQ=255)
                                                                                         cacaACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTAcc                                                                                                             >  2:201024/1‑91 (MQ=255)
                                                                                            aaCCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAAc                                                                       <  5:120820/126‑1 (MQ=255)
                                                                                            aaCCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAAc                                                                       >  6:120820/1‑126 (MQ=255)
                                                                                                     tACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGac                                       >  7:42695/1‑149 (MQ=255)
                                                                                                           gTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCgg                                 >  4:66562/1‑149 (MQ=255)
                                                                                                             ggAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTc                               <  1:183062/149‑1 (MQ=255)
                                                                                                             ggAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTc                               >  3:244884/1‑149 (MQ=255)
                                                                                                              gAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCCGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTc                               <  5:174063/148‑1 (MQ=255)
                                                                                                              gAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCCGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTc                               <  5:178425/148‑1 (MQ=255)
                                                                                                                   aGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAAc                                                  >  1:67518/1‑124 (MQ=255)
                                                                                                                   aGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAAc                                                  <  2:67518/124‑1 (MQ=255)
                                                                                                                             gACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACAt               >  5:178437/1‑149 (MQ=255)
                                                                                                                             gACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGAACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGAGCAAGGGCACCAACTTCCAGCAGACCACCGGACTGTACGTGCAGCAcct               >  7:48245/1‑147 (MQ=255)
                                                                                                                             gACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTATACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACGTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCGGGACAt               >  7:46941/1‑149 (MQ=255)
                                                                                                                                aaGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGAt            <  7:136292/149‑1 (MQ=255)
                                                                                                                                       aGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCcag                                            <  6:306203/110‑1 (MQ=255)
                                                                                                                                       aGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCcag                                            >  5:306203/1‑110 (MQ=255)
                                                                                                                                          catctcTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGCTGAAc  <  5:212860/144‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
AGCTCAACCGCGGCAATGTCGCCCGCGATGAGTACGGTGTGTTCAACCAGACGGAGCTCAAGCAACAAGCGCAACTGGCCGGCATGCAGCACAACCTGTTGTACGGCGTGGAAGCAGGCTTTCAGGACAAGTACCAGCGCGTGCTGAGCCAGACCAACGTGGCCCAGGTACCGGTGTACCGTGACGCCTTGGTGCCGGTACCGGAACATGCCGCCAACCTGTCGTCCAAGGGCACCAACTTCCAGCAGACCACCGGTCTGTACGTGCAGGACATGATCGAGCTGAAC  >  NC_002947/4104483‑4104769

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