Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,926 G→C 10.6% G501G (GGC→GGG) ‡ recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9260GC10.6% 185.3 / 10.5 56G501G (GGC→GGG) ‡recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base C (2/4);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

AAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATA  >  NC_000913/2960781‑2961067
                                                                                                                                                 |                                                                                                                                             
aaGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCg                                                                                                                                            <  1:402968/149‑1 (MQ=255)
   gTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                         <  1:119066/149‑1 (MQ=255)
    tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGa                                                                                                                                          >  2:290327/1‑147 (MQ=255)
    tAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGa                                                                                                                                          <  1:290327/147‑1 (MQ=255)
      gagccGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCGAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                      >  4:318799/4‑149 (MQ=255)
      gAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                      >  4:16424/1‑149 (MQ=255)
       aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                     <  1:417511/149‑1 (MQ=255)
       aaCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGt                                                                                                                                     <  1:369192/149‑1 (MQ=255)
        aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGtt                                                                                                                                    <  1:116804/149‑1 (MQ=255)
         ccGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGGCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTg                                                                                                                                   <  2:268760/149‑1 (MQ=255)
         ccGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTg                                                                                                                                   <  4:448772/149‑1 (MQ=255)
          cGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgg                                                                                                                                  <  2:377446/149‑2 (MQ=255)
                 cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                                                                                                           <  2:293005/149‑1 (MQ=255)
                 cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggggcata                                                                                                                            <  1:327628/148‑8 (MQ=255)
                 cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTggggcata                                                                                                                            >  2:327628/1‑141 (MQ=255)
                  aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                          >  3:33846/1‑149 (MQ=255)
                     cGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGt                                                                                                                       <  2:185220/149‑1 (MQ=255)
                          tGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                         >  3:414037/1‑126 (MQ=255)
                          tGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                         <  4:414037/126‑1 (MQ=255)
                               gcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCa                                                                                                             >  2:324376/1‑149 (MQ=255)
                                 tctcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGc                                                                                                           >  1:183804/1‑149 (MQ=255)
                                   tcCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGccc                                                                                                         >  2:195780/1‑149 (MQ=255)
                                     cATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAAc                                                                                                                <  4:195424/140‑1 (MQ=255)
                                     cATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAAc                                                                                                                >  3:195424/1‑140 (MQ=255)
                                     cATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCctt                                                                                                       >  4:312860/1‑149 (MQ=255)
                                              gCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                              <  3:318799/149‑1 (MQ=255)
                                               ctcaacaaCATACGCGTCAGGCCAAATCGCGAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGAt                                                                                             <  3:388275/147‑1 (MQ=255)
                                                  aacaacATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                          >  1:21757/1‑149 (MQ=255)
                                                                             cAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATc                                                               <  1:324376/149‑1 (MQ=255)
                                                                               ggTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCa                                                             >  1:23598/1‑149 (MQ=255)
                                                                                    gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAAc                                                        <  3:62547/149‑1 (MQ=255)
                                                                                    gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAAc                                                        >  3:97820/1‑149 (MQ=255)
                                                                                       gtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGcc                                                     >  1:261989/1‑149 (MQ=255)
                                                                                       gtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGcc                                                     >  3:443574/1‑149 (MQ=255)
                                                                                          cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGc                                                  >  2:123409/1‑149 (MQ=255)
                                                                                            gTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                              <  3:356003/103‑1 (MQ=255)
                                                                                            gTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                              >  4:356003/1‑103 (MQ=255)
                                                                                               gCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                             >  3:196808/1‑149 (MQ=255)
                                                                                                    gAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTg                                        <  2:291291/149‑1 (MQ=255)
                                                                                                                 cgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                           <  3:312860/149‑1 (MQ=255)
                                                                                                                 cgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACt                           <  3:421115/149‑1 (MQ=255)
                                                                                                                               aTCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa             >  1:87702/1‑149 (MQ=255)
                                                                                                                                tCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATaa            <  2:5635/149‑1 (MQ=255)
                                                                                                                                  tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCctct                                         <  3:123228/118‑1 (MQ=255)
                                                                                                                                  tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCctct                                         >  4:123228/1‑118 (MQ=255)
                                                                                                                                  tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACa          >  2:192862/1‑149 (MQ=255)
                                                                                                                                  tATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACa          >  2:369318/1‑149 (MQ=255)
                                                                                                                                     gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGc       <  1:123409/149‑1 (MQ=255)
                                                                                                                                      ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt      <  4:355535/149‑1 (MQ=255)
                                                                                                                                        tcAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg        <  1:461708/144‑1 (MQ=255)
                                                                                                                                        tcAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAgg        >  2:461708/2‑145 (MQ=255)
                                                                                                                                        tcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa    >  3:27650/2‑149 (MQ=255)
                                                                                                                                        ccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt              <  2:12985/139‑1 (MQ=255)
                                                                                                                                        ccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt              >  1:12985/1‑139 (MQ=255)
                                                                                                                                         cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt   >  3:427897/1‑149 (MQ=255)
                                                                                                                                          aaCGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATa  >  2:143937/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATA  >  NC_000913/2960781‑2961067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑