Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,527,591 C→T Q375* (CAG→TAG)  yrfF → putative RcsCDB‑response attenuator, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,527,5910CT97.8% 122.0 / NA 47Q375* (CAG→TAG) yrfFputative RcsCDB‑response attenuator, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (23/22);  minor base G (0/1);  total (23/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 = 9.76e-01

AAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGCAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTT  >  NC_000913/3527451‑3527728
                                                                                                                                            |                                                                                                                                          
aaGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGACTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAg                                                                                                                                          <  2:79688/143‑1 (MQ=255)
      aacaCTTTCCGTTACAGCACTGTCTGCTCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCAtt                                                                                                                                    <  2:259284/140‑1 (MQ=255)
       gagcTTTCCGGGACAGCACTGGCGGCGCAGTACGATTATCGCTGCGGGCTCGCTGATGGTGCTGTTTATGCTGGTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTg                                                                                                                                   <  1:406901/140‑1 (MQ=255)
       aagcTTTCCGTTGCGGCACGGGGTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTg                                                                                                                                   <  2:554275/140‑1 (MQ=255)
       aaaCTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTGATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCACTCATCGATGAAAGGCGCGTAGACCATTg                                                                                                                                   >  2:482771/1‑143 (MQ=255)
          cTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCAGTTATTCTGGGTCCCGCTGGATATGCCGCAGCAATCCACTCACTCATGGATGAAAGGCGAGTAGACCACTGAAg                                                                                                                                >  2:278177/1‑143 (MQ=255)
              ccGTTAGAGCACTGGCTGCGCAGTACGATTACCGCTGCGGGCTCGCTGCGGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCCCTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGccac                                                                                                                            <  1:677567/143‑1 (MQ=255)
              ccGTTACTGCCCTGTCTTCACAGTACGATTATCGCGGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGccac                                                                                                                            <  2:424932/143‑1 (MQ=255)
                       cACTGGCTGCGCAGTACGCTTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTaaa                                                                                                                   <  2:20522/143‑1 (MQ=255)
                          tGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCCATGAAGCCCCCAACGTAAaaca                                                                                                                >  2:602614/1‑143 (MQ=255)
                          tGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAaaca                                                                                                                <  1:281308/143‑1 (MQ=255)
                          tGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCACGCTGGATATGCCGCTGCAATTCGCTCTCTCATGGGTGAAGGGCGCGTAGACCATTGAAGCCGCCAGGGTAAaaca                                                                                                                >  1:341517/1‑143 (MQ=255)
                           ggcggcGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCGACCAGCGTAAaacaa                                                                                                               <  2:605066/143‑1 (MQ=255)
                            gCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAaacaac                                                                                                              <  2:712078/143‑1 (MQ=255)
                             ctccgcAGTTCGATTCTCGCTTCGGGCTCGCTACTGGTTCTGTTTATGCTGTTATTCTGGATCCCGCTGAATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACt                                                                                                             <  2:487992/140‑1 (MQ=255)
                                        gactATCGCTGCGGGCTCGCCGCTGGTGCTGTGTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTg                                                                                                  <  1:174328/140‑1 (MQ=255)
                                        gATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATCTCAGCCACCAGCGTAAAACAACCGGCTGAAGCTg                                                                                                  >  2:415312/1‑143 (MQ=255)
                                                cTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGAGCCCGCTGGATATGCCGCTGGAGTTCACTCTCTCATGGATGAAAGGCGCGTAGACCGTTGAAGCGAGCAGCGGAAGACACCGGGGTGGCGGTGGCGTGAgg                                                                                          >  2:2186/1‑143 (MQ=255)
                                                  gcgcgcTCGCGGCTTCTGCTGTTTATGCTGTTAGTCTGGATCCCGCTGGATATGCCGCTGAAATTCACGCCCTCCGGGATGCAACGCGCGTCGACCATTGAAGCCGCCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGt                                                                                        <  2:38253/143‑1 (MQ=255)
                                                      gCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTTTCATCGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGGCAAACAACTGGCTGTCGCTGCCGTGCGGCTCGGc                                                                                    >  1:66056/1‑143 (MQ=255)
                                                       cTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAACGGCGCGTAGACCATTGAAGCAACCAGAGAAAAACAACTGGATAACGCTGGCGTGTGGGTAGGCg                                                                                   >  1:228013/1‑143 (MQ=255)
                                                       cTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTCAACCATTGAAGCCCCCAGCGTACCACAACTGGCTGACGCTGGCGTGCGGGCAGGcc                                                                                   >  2:568174/1‑142 (MQ=255)
                                                        tCGCGGCTGGTGCTGTTTATGCTGTTATTCTGGCTCCCGCTGTATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGa                                                                                  <  2:43554/143‑1 (MQ=255)
                                                                 gTGCTGTTTATGCTGTTATACTGGAGCGCACTGGGTTTGGCTCTCAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCg                                                                        <  1:518530/143‑1 (MQ=255)
                                                                  tGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGt                                                                       >  2:149696/1‑143 (MQ=255)
                                                                   gCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTa                                                                      >  1:691716/1‑143 (MQ=255)
                                                                      gTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAATGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTAGCTGAAGCTGGCGTGCGGGTAGGCGATACATTGCGTATTa                                                                   >  1:10711/1‑143 (MQ=255)
                                                                       tttATGCTGTTATTCTAGATCCCGCTGGAGATGCCGCTGAAATTCACTCTCTCATGGATGAACGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCggcgtag                                                                  >  1:690442/1‑137 (MQ=255)
                                                                          acgCTGTTACTCTGGATCCCGCTCGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGCCCCTCGACGCCACTAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCgg                                                               <  1:268699/141‑1 (MQ=255)
                                                                              tGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGGACg                                                           >  1:33341/1‑143 (MQ=255)
                                                                               gTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACgg                                                          <  1:87698/143‑1 (MQ=255)
                                                                               gTTATTCTGGACCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCCTTGAAGCCACCAGCGTAAAACAACTGGCGGACGCTGGCGTGCGGGAAGGCGATAGACTGCCGTATAGCGGTACgg                                                          >  2:248837/1‑143 (MQ=255)
                                                                                     cTGGATACCGTTGGATATGCCGCTGAAATTCACGCTCTCGTGGCTGGAAGGCGCAGAGACGGTTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAAtgt                                                    <  1:300536/143‑1 (MQ=255)
                                                                                      tGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAAtgtg                                                   >  1:154119/1‑143 (MQ=255)
                                                                                                   ggacgCCGCTGGTATTCACTCTCTCGTGGATGCGGGGCGGGTGGCCCATTGAAGCCACCAGCGTAAGACAACTGGCTGGGGCGGGTTTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACt                                      <  2:471713/139‑1 (MQ=255)
                                                                                                         cGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAAATGGCTGACGCTGGCGCGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATACGAACTTCCGGa                                >  2:307215/1‑143 (MQ=255)
                                                                                                                 ttCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTCAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGCACTTCCGGAACCTGGAg                        >  2:267063/1‑143 (MQ=255)
                                                                                                                  tCACTCTCTCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAgc                       >  2:673149/1‑143 (MQ=255)
                                                                                                                  tCACTCCCTAATGGATTAAAGGCGCGTAGCCCATGGAGGCCACCAGCGTAAAACAACTGGCTGACGCGGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAgc                       <  2:113735/143‑1 (MQ=255)
                                                                                                                     ctctctCATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAACATTCGAAGTTCCGCAACCTGGAGcgcg                    >  1:697296/1‑143 (MQ=255)
                                                                                                                      tctctcATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAAGGTGAAATATTAGAACTTCCGGAACCAAGAgcgcaa                   >  2:50963/1‑141 (MQ=255)
                                                                                                                       ctctcATGGATGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGaa                  <  1:369434/143‑1 (MQ=255)
                                                                                                                         ccaaTGGCTGAACGGCGCGTAGGCCATTGGCGCCACCAGCGTAAAACAACTGGCGTACGCTGGCGTGCGGGTAGGCGATACATTGCGTATGAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGaaaa                <  2:29401/140‑1 (MQ=255)
                                                                                                                               gATGAAGGGCGCGTAGCGCATTGAAGCCACCAGCGTACAGGAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAAtt          <  1:227471/143‑1 (MQ=255)
                                                                                                                               catacacggCGCGCAGCCCATGAAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAAtt          <  2:610751/136‑1 (MQ=255)
                                                                                                                                aTGAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAGTATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTc         >  1:447239/1‑143 (MQ=255)
                                                                                                                                  gAAAGGCGCGTAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGAGCATTAGTCCCACTTGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCACTTccc       >  1:163435/1‑141 (MQ=255)
                                                                                                                                       gcgcgtAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttt  <  1:149696/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
AAGTAAAAAACTTTCCGTTACAGCACTGGCTGCGCAGTACGATTATCGCTGCGGGCTCGCTGCTGGTGCTGTTTATGCTGTTATTCTGGATCCCGCTGGATATGCCGCTGAAATTCACTCTCTCATGGATGAAAGGCGCGCAGACCATTGAAGCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTT  >  NC_000913/3527451‑3527728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: