Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,951,877 G→A 6.8% A928V (GCT→GTT)  recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,951,8770GA 6.8% 233.8 / 4.9 73A928V (GCT→GTT) recBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base A (2/3);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTA  >  AP009048/2951808‑2951944
                                                                     |                                                                   
ggTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAAAGACGCTGGCAACGCCTgcagc                                                                    <  1:3226302/71‑1 (MQ=255)
 gTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcagca                                                                   <  1:1288309/71‑1 (MQ=255)
 gTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTgcaaca                                                                   <  1:3232997/71‑1 (MQ=255)
   gACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATc                                                                 >  1:1513721/1‑71 (MQ=255)
        aCCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt                                                            >  1:1673003/1‑71 (MQ=255)
           gcgcgGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATc                                                           >  1:2512268/1‑69 (MQ=255)
               ggAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGcc                                                      <  1:2305506/70‑1 (MQ=255)
                 aaaCTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAg                                                   >  1:3045005/1‑71 (MQ=255)
                     tGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACAt                                                            >  1:1837598/1‑58 (MQ=255)
                       aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc                                              >  1:1940694/1‑70 (MQ=255)
                       aTGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATc                                              <  1:4128131/70‑1 (MQ=255)
                        tgtgGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCa                                             >  1:1236620/1‑70 (MQ=255)
                             tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                       >  1:1100537/1‑71 (MQ=255)
                             tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                       >  1:3063817/1‑71 (MQ=255)
                             tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                       <  1:921802/71‑1 (MQ=255)
                             tgtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCt                                       <  1:597560/71‑1 (MQ=255)
                              gtTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAg                                                   <  1:949905/58‑1 (MQ=255)
                               ttAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCa                                             <  1:2104561/63‑1 (MQ=255)
                                 aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa                                                    >  1:1998254/1‑54 (MQ=255)
                                 aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGc                                   >  1:862568/1‑71 (MQ=255)
                                 aaCGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCACCCGAGGCATCAAATCCTGGGc                                   >  1:410971/1‑71 (MQ=255)
                                   cGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGa                                 >  1:4310832/1‑71 (MQ=255)
                                    gTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGAt                                <  1:3744279/71‑1 (MQ=255)
                                      cGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAc                              >  1:3837762/1‑71 (MQ=255)
                                       ggTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCaa                                            >  1:2732000/1‑56 (MQ=255)
                                       ggTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAcc                             >  1:1485063/1‑71 (MQ=255)
                                       ggTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTg                                      <  1:3268711/62‑1 (MQ=255)
                                        gTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg                            <  1:3844001/71‑1 (MQ=255)
                                        gTTCTTCAACGACGCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATc                                         >  1:3065334/1‑58 (MQ=255)
                                         ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATc                                         <  1:2008777/57‑1 (MQ=255)
                                         ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           <  1:4419994/71‑1 (MQ=255)
                                         ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           <  1:933213/71‑1 (MQ=255)
                                         ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           <  1:4434778/71‑1 (MQ=255)
                                         ttcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           <  1:3484423/71‑1 (MQ=255)
                                          tcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           >  1:2783076/1‑70 (MQ=255)
                                          tcttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGt                           >  1:4468824/1‑70 (MQ=255)
                                           cttcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGa                                                    >  1:448918/1‑44 (MQ=255)
                                             tcAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc                        >  1:1719323/1‑70 (MQ=255)
                                              cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTg                                      <  1:3695344/55‑1 (MQ=255)
                                              cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAcc                       <  1:1336076/70‑1 (MQ=255)
                                              cAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGAc                        <  1:159247/69‑1 (MQ=255)
                                              cAACGACGCTGGCAACGCCTGCAGCATCGACATACAGCCGAGGCATCAAATCCTg                                      <  1:1686209/55‑1 (MQ=255)
                                               aaCGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCt                                       <  1:2845138/53‑1 (MQ=255)
                                                acgacgCTGGCAACGCCTGCAGCATCGACATCAAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACg                    >  1:1390546/1‑71 (MQ=255)
                                                 cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGa                                 >  1:842361/1‑57 (MQ=255)
                                                 cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc                   >  1:942023/1‑71 (MQ=255)
                                                 cgacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGc                   <  1:238887/71‑1 (MQ=255)
                                                  gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAgga                                                 >  1:3227710/1‑39 (MQ=255)
                                                  gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCt                  <  1:2813839/71‑1 (MQ=255)
                                                  gacgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCt                  <  1:996566/71‑1 (MQ=255)
                                                   acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGc                                                 >  1:4468753/1‑39 (MQ=255)
                                                   acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg                 <  1:868808/71‑1 (MQ=255)
                                                   acgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg                 <  1:1143069/71‑1 (MQ=255)
                                                   acgCTGGCAACGCCTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTg                 <  1:3541470/71‑1 (MQ=255)
                                                    cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt                <  1:3578991/71‑1 (MQ=255)
                                                    cgCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt                <  1:3195975/71‑1 (MQ=255)
                                                     gCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATcc                                        >  1:3954324/1‑46 (MQ=255)
                                                      cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg              <  1:1542637/71‑1 (MQ=255)
                                                      cTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTg              >  1:1864046/1‑71 (MQ=255)
                                                       tGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg                            >  1:1888447/1‑56 (MQ=255)
                                                         gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCg                            <  1:4530258/54‑1 (MQ=255)
                                                         gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa            >  1:3955084/1‑70 (MQ=255)
                                                         gCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa            >  1:4075723/1‑70 (MQ=255)
                                                           aaCGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTg                          >  1:1020327/1‑54 (MQ=255)
                                                             cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCaaa          >  1:220708/1‑68 (MQ=255)
                                                             cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCa       <  1:4314943/71‑1 (MQ=255)
                                                             cGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCa       <  1:2621859/71‑1 (MQ=255)
                                                              gCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCag      <  1:1941353/71‑1 (MQ=255)
                                                               ccTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGt                >  1:1719946/1‑60 (MQ=255)
                                                               ccTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCaga     <  1:4300510/71‑1 (MQ=255)
                                                                cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATAcc                             <  1:2821657/46‑1 (MQ=255)
                                                                cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa            >  1:634707/1‑63 (MQ=255)
                                                                cTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCagag    >  1:3638036/1‑71 (MQ=255)
                                                                cTGCAACATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCaga     >  1:1336670/1‑70 (MQ=255)
                                                                 tGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt               >  1:2102890/1‑59 (MQ=255)
                                                                 tGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCa            <  1:3730080/62‑1 (MQ=255)
                                                                  gcagcaTCGACATCCAGCCAAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCaaa          <  1:1521639/63‑1 (MQ=255)
                                                                  gcaacaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGtt               >  1:2908658/1‑58 (MQ=255)
                                                                   cagcaTCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTa  <  1:688587/70‑1 (MQ=255)
                                                                     |                                                                   
GGTGACGCACCGCGCGGAAACTGATGTGGTGTTAACGTCGGTTCTTCAACGACGCTGGCAACGCCTGCAGCATCGACATCCAGCCGAGGCATCAAATCCTGGGCGATACCGTGACCACGCTGTTGCAAACCAGAGTA  >  AP009048/2951808‑2951944

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