Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,407 G→C 7.6% R22P (CGC→CCC)  uvrD → DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,991,4070GC 7.6% 2317.4 / 24.1 735R22P (CGC→CCC) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (230/435);  new base C (67/2);  total (297/438)
Fisher's exact test for biased strand distribution p-value = 3.36e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 793 total aligned reads displayed.

CAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGC  >  CP009273/3991341‑3991473
                                                                  |                                                                  
cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACg                                                                    <  1:39073874/67‑1 (MQ=255)
cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACg                                                                    <  1:37837458/67‑1 (MQ=255)
cAATGGACGGTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGCGGCCGCGCCACg                                                                    <  1:42271430/67‑1 (MQ=255)
  aTGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgcag                                                                 <  1:42953244/68‑1 (MQ=255)
   tGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgcagc                                                                >  1:38654231/1‑68 (MQ=255)
     gACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACccagcc                                                               >  1:5333693/1‑66 (MQ=255)
      aCGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCaa                                                              >  1:10022570/1‑67 (MQ=255)
       cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAcc                                                            <  1:29581333/68‑1 (MQ=255)
       cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                             <  1:31511866/67‑1 (MQ=255)
        gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                           <  1:7214807/68‑1 (MQ=255)
        gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                           <  1:7539343/68‑1 (MQ=255)
        gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                           <  1:31900428/68‑1 (MQ=255)
           tCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCt                                                        <  1:35846013/68‑1 (MQ=255)
              tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                     <  1:12155468/68‑1 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACATCGCGAAGCGGTGGCCGCGCCACCCACCAAACCTCAGGTGc                                                   >  1:4833678/1‑68 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTgg                                                   >  1:38471979/1‑67 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCGGGt                                                     >  1:26538902/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCCACCTTCTGGt                                                     >  1:1746090/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCTTGTggg                                                     >  1:22695512/1‑65 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCTTGTGGt                                                     >  1:42608450/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                     >  1:5631300/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                     >  1:14917425/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                     >  1:29928777/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTGTGGTGc                                                   >  1:34723895/1‑68 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCCACCTTCAGGt                                                     >  1:25702892/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCACCCACAAACCATCTGGTGc                                                   >  1:41931565/1‑68 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCCCGCCACGCAGCAACCTTCAGGt                                                     >  1:30203011/1‑66 (MQ=255)
                ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCCCCCCACGCCGCCACCGTCTGGTGc                                                   >  1:5969479/1‑68 (MQ=255)
                 cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCATATGGt                                                     >  1:40255126/1‑65 (MQ=255)
                  tGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTg                                                 >  1:41012114/1‑68 (MQ=255)
                        aCAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGggc                                           <  1:32372264‑M2/68‑4 (MQ=255)
                        aCAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAACTTCTCGTGGTGGCGggc                                           <  1:15201116‑M2/68‑4 (MQ=255)
                             cTTAATGACAAACCGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGggcg                                          >  1:38027476‑M2/1‑60 (MQ=255)
                             cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGggcg                                          >  1:35618755‑M2/1‑60 (MQ=255)
                             cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGGGGTGGCGggcgcgg                                       <  1:13221608‑M2/67‑8 (MQ=255)
                               tAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCa                                    <  1:36862823/68‑1 (MQ=255)
                               tAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCa                                    <  1:16546131/68‑1 (MQ=255)
                                 aTGCCACACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATCTGGAGCTGGCTGGCGCGGGCAGt                                  <  1:2020153/68‑1 (MQ=255)
                                 aTGACAACCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:8201022/68‑1 (MQ=255)
                                 aTGACAAACCGCCCGACGCGGTGGCCTCGCAACGCCGCAACCTGCTGGCGCTGGCGGGCGCGGGCAGt                                  <  1:10495098/68‑1 (MQ=25)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:17879056/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:18152156/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:19064416/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:5086423/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:39518202/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:11042628/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:33284560/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:21998388/68‑1 (MQ=255)
                                 aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCATCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:2723550/68‑1 (MQ=255)
                                 aTGACAAACACCGCGACGCGGGGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                  <  1:23273622/68‑1 (MQ=38)
                                  tGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCCGTg                                 >  1:23774546/1‑68 (MQ=38)
                                  tGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCACCAAACTTCTTGTGCTGGCGGGCGCGGGCAGTg                                 >  1:28499895/1‑68 (MQ=255)
                                   gACAACCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                <  1:15459041/68‑1 (MQ=255)
                                   gACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCCGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                >  1:24303350/1‑68 (MQ=255)
                                   gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                <  1:13036136/68‑1 (MQ=255)
                                   gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                <  1:11498687/68‑1 (MQ=255)
                                   gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                <  1:1593436/68‑1 (MQ=255)
                                      aaaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTaa                             >  1:33811167/1‑68 (MQ=255)
                                       aaCAGCGCGAAGCGGTGGCCGCGCCCCGCAGGAACCTTCTGGTGCGGGCGGGCGCGGGCAGTGGTTAg                            >  1:21048882/1‑68 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCCTGTGGTGCTGGCGGGCGCGGGGAGTGGTAAg                            >  1:21401776/1‑67 (MQ=38)
                                        aCAGCGCGAAGCGGTGGCCGCGCCCCCCCGCAACCGTGTGGTGGTGGCGGGCGCGGGCCGTGGTAAg                            >  1:7927623/1‑67 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                            >  1:13689282/1‑67 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                            >  1:28042139/1‑67 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                            >  1:27502144/1‑67 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCGCCCCACCCACCAACCGTCTGGTGCTGGCGGGCGCCGGCAGTGGTAAg                            >  1:26915413/1‑67 (MQ=255)
                                        aCAGCGCGAAGCGGTGGCCCCGCCACCCAGCAACCTTCTGGTGGTGGCGGGCGCGGGCAGTGGGAAGa                           >  1:14407376/1‑68 (MQ=255)
                                         cAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                          >  1:6695279/1‑68 (MQ=255)
                                          cgcgcGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                         <  1:24742048/67‑1 (MQ=255)
                                          aGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                          >  1:9313682/1‑67 (MQ=255)
                                            cgcgAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgc                        >  1:26087114/1‑67 (MQ=255)
                                             gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                      <  1:19890455/68‑1 (MQ=255)
                                             gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                      <  1:36233556/68‑1 (MQ=255)
                                             gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                      <  1:25365815/68‑1 (MQ=255)
                                             gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgc                        <  1:6849585/66‑1 (MQ=255)
                                             gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                         <  1:20604359/65‑1 (MQ=255)
                                               gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                    <  1:9670016/68‑1 (MQ=255)
                                               gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                    <  1:10813322/68‑1 (MQ=255)
                                               gAAGCGGGGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                    <  1:40631541/68‑1 (MQ=255)
                                                 aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                   <  1:3281770/67‑1 (MQ=255)
                                                  gCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                    >  1:18413862/1‑65 (MQ=255)
                                                   cactgccCGCGCCACGCAGCAACCTTCTGGTGCTGGCGggcgcggccaccgcttcgcgctgtttgtca                >  1:39836484‑M2/4‑38 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                 <  1:12493586/66‑1 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                 <  1:8974152/66‑1 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                 <  1:30646997/66‑1 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                 <  1:18835381/66‑1 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg               <  1:34219072/68‑1 (MQ=255)
                                                    ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg               <  1:1408603/68‑1 (MQ=255)
                                                      tgcccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt              >  1:8611659/4‑67 (MQ=255)
                                                      tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt              <  1:33480556/67‑1 (MQ=255)
                                                       ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc            <  1:31200695/68‑1 (MQ=255)
                                                          cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATc         <  1:22915708/68‑1 (MQ=255)
                                                           gcgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCg        <  1:29111734/68‑1 (MQ=255)
                                                           gcgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCg        <  1:35115635/68‑1 (MQ=255)
                                                           gcgcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCg        <  1:41384/68‑1 (MQ=255)
                                                             gcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt       <  1:19475594/67‑1 (MQ=255)
                                                             gcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt       <  1:6326637/67‑1 (MQ=255)
                                                             gcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt       <  1:16951027/67‑1 (MQ=255)
                                                             gcCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt       <  1:4592937/67‑1 (MQ=255)
                                                                 cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGc  <  1:433669/68‑1 (MQ=255)
                                                                 cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGc  <  1:37422439/68‑1 (MQ=255)
                                                                  |                                                                  
CAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGC  >  CP009273/3991341‑3991473

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