Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,235,585 T→C 8.7% W307R (TGG→CGG)  fecA → outer membrane ferric citrate porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,235,5850TC 8.7% 125.3 / 2.2 46W307R (TGG→CGG) fecAouter membrane ferric citrate porin
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 = 1.00e+00

GCCACGTTCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGAT  >  NC_002947/5235438‑5235731
                                                                                                                                                   |                                                                                                                                                  
gCCACGTTCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTg                                                                                                                                                   >  1:305150/1‑149 (MQ=255)
gCCACGTTCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTAATCCGGCACCCGCGGCGGTGACTg                                                                                                                                                   >  5:260649/1‑149 (MQ=255)
gCCACGTTCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAACCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTg                                                                                                                                                   >  1:78885/1‑149 (MQ=255)
       tCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGaa                                                                                                                                            <  4:195938/149‑1 (MQ=255)
         aaGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACCCCGGCACCCGCGGGGGTGACTGGCGCGAaca                                                                                                                                          >  6:299930/1‑149 (MQ=255)
          aGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGTCACCCGCGGGTGTGACTGGCGCGAaca                                                                                                                                          >  4:296883/1‑148 (MQ=255)
                  aTGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGacac                                                                                                                                 <  8:259681/149‑1 (MQ=255)
                     cAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGc                                                                                                                               <  6:305301/148‑1 (MQ=255)
                                     cGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACCCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATc                                                                                                              <  5:299930/149‑1 (MQ=255)
                                      gTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCCACCTGCTGCTCGGTGGCACCAACGACAACGCTCAGCGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGCGTCACTGGCGCGAACACCGCGACACGCAGATCGACGACCTGATcc                                                                                                             >  7:131474/1‑149 (MQ=255)
                                           cGAGCCACGAAGGTGTCACCTCCAGCGCCATCCTGCGGGTCGGCGGCACCTACGACAACGGTCTGGGCGGCGCCGTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTcta                                                                                                         <  8:131474/148‑3 (MQ=255)
                                            gAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAg                                                                                                        >  2:316761/1‑148 (MQ=255)
                                            gAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAg                                                                                                        <  1:316761/148‑1 (MQ=255)
                                              gCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCgacgac                                                                                                                    >  5:326780/1‑134 (MQ=255)
                                              gCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCgacgac                                                                                                                    <  6:326780/134‑1 (MQ=255)
                                               ccACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGc                                                                                                     >  6:135370/1‑148 (MQ=255)
                                                         ttCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACt                                                                                          >  5:126083/1‑149 (MQ=255)
                                                            aaGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGgc                                                                                                                                                 <  5:253583/91‑1 (MQ=255)
                                                            aaGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGgc                                                                                                                                                 >  6:253583/1‑91 (MQ=255)
                                                            aaGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTgg                                                                                        <  5:135370/148‑1 (MQ=255)
                                                            aaGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTgg                                                                                        <  2:305150/148‑1 (MQ=255)
                                                            aaGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGa                                                                                       <  6:126083/149‑1 (MQ=255)
                                                                            tgctgGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAACGGCAAGCTGCAACTGGACGAAGCCAACAGCCTg                                                                       >  4:301472/1‑149 (MQ=255)
                                                                               tgGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAgg                                                                                                       >  6:339561/1‑114 (MQ=255)
                                                                               tgGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAgg                                                                                                       <  5:339561/114‑1 (MQ=255)
                                                                                    ggTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACCGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCAt                                                               >  2:343600/1‑149 (MQ=255)
                                                                                    ggTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACCGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCAt                                                               >  6:91870/1‑149 (MQ=255)
                                                                                               cGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGtact                                                    <  2:78885/149‑1 (MQ=255)
                                                                                                   aaCGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTTGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGa                                                <  3:301472/149‑1 (MQ=255)
                                                                                                     cGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAgg                                              >  7:299590/1‑149 (MQ=255)
                                                                                                             gcggcgCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGtactac                                                  >  6:105073/1‑137 (MQ=255)
                                                                                                             gcggcgCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGtactac                                                  <  5:105073/137‑1 (MQ=255)
                                                                                                                 cgcCCTGTTGTACTCCGGCACCCGCGGGGGTGACCGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACa                                  <  5:91870/149‑1 (MQ=255)
                                                                                                                 cgcCCTGTTGTACTCCGGCACCCGCGGGGGTGACCGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACa                                  <  1:343600/149‑1 (MQ=255)
                                                                                                                          gtACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCcggcg                         >  8:290138/1‑149 (MQ=255)
                                                                                                                             cTCCGGCACCCGCGGGGGTGACTGGCTCGAACACAGCGACACGCACATCTACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCggcggc                       <  6:260649/148‑1 (MQ=255)
                                                                                                                              tCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGc                                            >  6:20485/1‑126 (MQ=255)
                                                                                                                              tCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGc                                            <  5:20485/126‑1 (MQ=255)
                                                                                                                              tCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCt                     >  7:206930/1‑149 (MQ=255)
                                                                                                                              tCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTAATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCt                     >  1:271031/1‑149 (MQ=255)
                                                                                                                                cGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGACTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCt                     >  6:162857/1‑147 (MQ=255)
                                                                                                                                  gCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAg                  >  6:120838/1‑148 (MQ=255)
                                                                                                                                    aCCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCAc               >  4:342577/1‑149 (MQ=255)
                                                                                                                                           ggggTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGAt        <  5:120838/149‑1 (MQ=255)
                                                                                                                                              gTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTc     >  1:171871/1‑149 (MQ=255)
                                                                                                                                                 aCTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGAt  >  6:335431/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCACGTTCAAGGCCGCCATGCAAAACCAGATCAGCCCGTCTTCGAGCCACGAAGGTTTCAAGAACAGCGCCAACCTGCTGGTCGGTGGCACCAACGACAACGGTCTGGGCGGCGCCCTGTTGTACTCCGGCACCCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGAT  >  NC_002947/5235438‑5235731

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