Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,955,163 G→A 10.6% Q281* (CAG→TAG)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,955,1630GA10.6% 131.7 / 4.8 46Q281* (CAG→TAG) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (16/25);  new base A (2/3);  total (18/28)
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

CCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACC  >  NC_000913/2955025‑2955311
                                                                                                                                          |                                                                                                                                                    
ccAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                                                                                                                            <  2:563392/149‑1 (MQ=255)
ccAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                                                                                                                            <  4:226001/149‑1 (MQ=255)
ccAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGttt                                                                                                                                            <  4:465659/149‑1 (MQ=255)
                 aaTGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACCGTGGATGTCGCGGGGTTTCCCCCCCGGCCCTCGTGCGTTCTTCTAGGAAACGGTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTCCCCa                                                                                                                           >  1:330695/1‑149 (MQ=255)
                 aaTGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCa                                                                                                                           >  4:324279/1‑149 (MQ=255)
                    ggTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGgc                                                                                                                        >  2:466751/1‑149 (MQ=255)
                         tGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGa                                                                                                                   <  3:310917/149‑1 (MQ=255)
                         tGCAAGCAGTGGATCGATCTCCTCAAACATTGGATGTCGCGGGCTTGCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTAATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGa                                                                                                                   <  2:6831/149‑1 (MQ=255)
                            aaGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCt                                                                                                                >  3:126765/1‑149 (MQ=255)
                               cAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGt                                                                                                             <  4:126765/149‑1 (MQ=255)
                               cAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGt                                                                                                             <  2:102778/149‑1 (MQ=255)
                                   tGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGAt                                                                                                         <  4:57094/149‑1 (MQ=255)
                                   tGATCAATCGCTCCAAACAGTGGATGTCCCGGGGTTT‑CCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATc                                                                                                        <  3:413476/149‑1 (MQ=255)
                                    gatcgatcGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATc                                                                                                        <  2:278547/149‑1 (MQ=255)
                                    gatcgatcGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATc                                                                                                        <  1:354944/149‑1 (MQ=255)
                                      tcgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCa                                                                                                      >  2:280531/1‑149 (MQ=255)
                                        gatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                                                                                                                       >  4:168341/1‑114 (MQ=255)
                                        gatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                                                                                                                       <  3:168341/114‑1 (MQ=255)
                                        gatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCAtt                                                                                                    <  3:423051/149‑1 (MQ=255)
                                              cTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTAATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                                                                                              >  4:13437/1‑149 (MQ=255)
                                              cTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTAATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCt                                                                                              >  4:108679/1‑149 (MQ=255)
                                                  aaCAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGAtt                                                                                          >  3:158169/1‑149 (MQ=255)
                                                       tGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTAc                                                                                     <  1:466751/149‑1 (MQ=255)
                                                       tGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTAc                                                                                     <  4:568168/149‑1 (MQ=255)
                                                              cgcgGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTaaa                                                                              <  1:247180/149‑1 (MQ=255)
                                                               gcgGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATc                                                                                                                 >  1:347033/1‑113 (MQ=255)
                                                               gcgGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATc                                                                                                                 <  2:347033/113‑1 (MQ=255)
                                                                           ccccGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCa                                                                 >  1:446612/1‑149 (MQ=255)
                                                                              cGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTAATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATa                                                              <  3:108679/149‑1 (MQ=255)
                                                                              cGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTAATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATa                                                              <  3:13437/149‑1 (MQ=255)
                                                                                 ccTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCa                                                           <  2:97950/149‑1 (MQ=255)
                                                                                  cTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCag                                                          <  3:56850/149‑1 (MQ=255)
                                                                                              ttctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCa                                              >  4:167786/1‑149 (MQ=255)
                                                                                                ctAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAgc                                            <  2:276135/149‑1 (MQ=255)
                                                                                                   aGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCa                                         >  2:74755/1‑149 (MQ=255)
                                                                                                     aaaCGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATc                                       >  2:224622/1‑149 (MQ=255)
                                                                                                        cGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAg                                    <  1:131523/149‑1 (MQ=255)
                                                                                                              gagaATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCAcc                              >  4:289078/1‑149 (MQ=255)
                                                                                                               agaATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCAccc                             <  1:490758/149‑1 (MQ=255)
                                                                                                                   tttttCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTg                         >  2:311520/1‑149 (MQ=255)
                                                                                                                   tttttCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTg                         >  3:204994/1‑149 (MQ=255)
                                                                                                                       tCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTc                     <  1:224622/149‑1 (MQ=255)
                                                                                                                                   ggCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCActgctg         >  2:231160/1‑149 (MQ=255)
                                                                                                                                   ggCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCActgctg         <  1:246339/149‑1 (MQ=255)
                                                                                                                                    gCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGt        <  2:494474/149‑1 (MQ=255)
                                                                                                                                     cAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGtt       <  4:158169/149‑1 (MQ=255)
                                                                                                                                          gATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTCGCCTTATTGCTACGGTTAAACTTACGTCGCTCAATACCAGAAGATTCGATCAGCGCATCCAGGTCACCCACTGCGCCGCGCCACTGCTGTTTTAcc  >  2:387408/1‑149 (MQ=255)
                                                                                                                                          gATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTAcc  >  4:169737/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGCCTGATTGCTACGGTTAAACTTGCGTCGATCAATACCAGAAGATTCGATCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACC  >  NC_000913/2955025‑2955311

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