Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,955,062 A→T D314E (GAT→GAA 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,0620AT98.3% 207.0 / NA 60D314E (GAT→GAArecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (33/26);  minor base C (0/1);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGG  >  NC_000913/2954965‑2955157
                                                                                                 |                                                                                               
cACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTgttc                                                                                               <  1:1038114/100‑1 (MQ=255)
   ggcggcGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTGgttcgat                                                                                            <  2:840341/100‑1 (MQ=255)
       gcgTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGcc                                                                                        <  2:142665/100‑1 (MQ=255)
       gcgTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGcc                                                                                        >  1:1303565/1‑100 (MQ=255)
        cgTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCt                                                                                       >  1:355089/1‑100 (MQ=255)
             ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAAc                                                                                  >  1:352755/1‑100 (MQ=255)
             ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAAc                                                                                  <  2:1174322/100‑1 (MQ=255)
                 cgcgCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTg                                                                              >  1:1425094/1‑100 (MQ=255)
                 cgcgCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTg                                                                              >  2:249/1‑100 (MQ=255)
                   cgcTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGa                                                                            >  2:413443/1‑100 (MQ=255)
                       aCTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTc                                                                        >  2:1383847/1‑100 (MQ=255)
                        cTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTcg                                                                       >  2:315711/1‑100 (MQ=255)
                             tCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCgggg                                                                  >  1:1329379/1‑100 (MQ=255)
                              cgcgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGTAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGt                                                                 <  1:984886/100‑1 (MQ=255)
                              cgcgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGt                                                                 >  2:509377/1‑100 (MQ=255)
                                cgGATCTCAGCCAATGCGCGGGTGATCGCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGAGGGGttt                                                               <  1:521508/100‑1 (MQ=255)
                                 ggATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTc                                                              <  2:497206/100‑1 (MQ=255)
                                 ggATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTc                                                              >  2:844394/1‑100 (MQ=255)
                                  gATCTCAGCCAATGCTCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTcc                                                             <  2:1578084/100‑1 (MQ=255)
                                  gATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTcc                                                             >  2:1047154/1‑100 (MQ=255)
                                  gATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTcc                                                             <  1:1383833/100‑1 (MQ=255)
                                   aTCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccc                                                            >  2:1063021/1‑100 (MQ=255)
                                    tctcAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTcccc                                                           <  2:607466/100‑1 (MQ=255)
                                      tcAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTcccccc                                                         >  2:540479/1‑100 (MQ=255)
                                       cAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccccccc                                                        >  1:221116/1‑100 (MQ=255)
                                       cAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccccccc                                                        >  1:1154775/1‑100 (MQ=255)
                                       cAGCCAATGCGCGGGTGATCACCAGACCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccccccc                                                        >  1:266133/1‑100 (MQ=255)
                                        aGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCg                                                       >  2:1147230/1‑100 (MQ=255)
                                        aGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCg                                                       <  2:1550957/100‑1 (MQ=255)
                                           cAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGAGGGGTTTCCCCCCCGGcc                                                    <  1:179685/100‑1 (MQ=255)
                                              tGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAATCAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTc                                                 <  2:222265/100‑1 (MQ=255)
                                               gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACCGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCg                                                >  1:391184/1‑100 (MQ=255)
                                                cgcgGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGt                                               >  1:180438/1‑100 (MQ=255)
                                                 gcgGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTg                                              >  2:482734/1‑100 (MQ=255)
                                                   gggTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCg                                            <  2:1372856/100‑1 (MQ=255)
                                                   gggTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCg                                            <  1:1145361/100‑1 (MQ=255)
                                                    ggTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGa                                           >  2:1223311/1‑100 (MQ=255)
                                                     gTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAt                                          >  2:1456177/1‑100 (MQ=255)
                                                          cACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttct                                     >  1:669221/1‑100 (MQ=255)
                                                           aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTa                                    <  1:482730/100‑1 (MQ=255)
                                                           aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTa                                    <  2:409299/100‑1 (MQ=255)
                                                             cAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAg                                  <  2:727749/100‑1 (MQ=255)
                                                                aTCGCGGATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGaaa                               >  2:49901/1‑100 (MQ=255)
                                                                      gATCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTgg                         >  2:1197410/1‑100 (MQ=255)
                                                                       aTCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTggg                        <  2:1010667/100‑1 (MQ=255)
                                                                       aTCGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTggg                        <  1:413439/100‑1 (MQ=255)
                                                                         cGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGgag                      >  1:1044034/1‑100 (MQ=255)
                                                                         cGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGgag                      <  1:844387/100‑1 (MQ=255)
                                                                         cGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGgag                      >  1:1543197/1‑100 (MQ=255)
                                                                         cGACAATGGTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGgag                      >  1:1238324/1‑100 (MQ=255)
                                                                                ggTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt               <  2:1425108/100‑1 (MQ=255)
                                                                                ggTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt               <  2:526644/100‑1 (MQ=255)
                                                                                 gTTCTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTc              <  2:894690/100‑1 (MQ=255)
                                                                                  ttCTGCAAGCAGTTGCTCGATCGCCTCAAACAGTGGATGTAGCGGGGGTTCCCCCCCGGCCTTTGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTcc             <  1:152867/100‑1 (MQ=255)
                                                                                  tccTGCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTcc             >  2:1082578/3‑100 (MQ=255)
                                                                                      gCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCg         <  2:1303579/100‑1 (MQ=255)
                                                                                      gCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCg         >  1:996391/1‑100 (MQ=255)
                                                                                      gCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCg         <  1:1545143/100‑1 (MQ=255)
                                                                                      gCAAGCAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCg         <  1:1063010/100‑1 (MQ=255)
                                                                                           cAGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTcc    <  2:1212974/100‑1 (MQ=255)
                                                                                            aGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGCATTTTTCCCGCGACTCCg   >  1:1300197/1‑100 (MQ=255)
                                                                                            aGTTGTTCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCg   >  2:499035/1‑100 (MQ=255)
                                                                                             gttgttCGATCGCCTCAAACAGTGGAGGTCGCGGGGTTTCCCCCCCGCCCTTCGTGCGATCTTCGAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCgg  <  1:1448006/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGG  >  NC_000913/2954965‑2955157

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