Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,416 T→G 8.6% L25R (CTT→CGT)  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,4160TG 8.6% 2435.3 / 35.3 932L25R (CTT→CGT) uvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base T (282/522);  new base G (83/8);  total (401/531)
Fisher's exact test for biased strand distribution p-value = 2.76e-26
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 967 total aligned reads displayed.

GTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGAT  >  CP009273/3991349‑3991482
                                                                   |                                                                  
gTTTCTTGCCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                                    <  1:27263533/68‑1 (MQ=255)
gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                                    <  1:16223218/68‑1 (MQ=255)
      tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCAGCCACCATGTGGt                                                              >  1:12044256/1‑68 (MQ=255)
      tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              <  1:5150293/68‑1 (MQ=255)
      tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              <  1:20112683/68‑1 (MQ=255)
      tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              <  1:31370502/68‑1 (MQ=255)
      tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              <  1:33324919/68‑1 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGt                                                              >  1:25213411/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCCCAACCGTCTGGTgg                                                            >  1:21675886/1‑67 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              >  1:50815009/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              >  1:50572013/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              >  1:13571558/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              >  1:29200054/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                              >  1:16895953/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTATGGt                                                              >  1:46459546/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCACCCACCAACCGTCTGGt                                                              >  1:35034574/1‑66 (MQ=255)
        ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCCCGCCACGCAGCAACCGTCGGGTGc                                                            >  1:49462578/1‑68 (MQ=255)
        acTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCACCAACCGTCTGGTgg                                                            >  1:19580516/2‑67 (MQ=255)
          tGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTg                                                          <  1:32271174/68‑1 (MQ=255)
           gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTgg                                                         <  1:11113287/68‑1 (MQ=255)
               gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCAGCCACCTTCTGGTGCTggcgg                                                      >  1:10017268/1‑67 (MQ=39)
               gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                                     <  1:35790330/68‑1 (MQ=255)
                  aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGGTGgcgggcgc                                                  >  1:47547128/1‑68 (MQ=255)
                  aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCGGGTGCTGgcgggcgc                                                  >  1:24728913/1‑68 (MQ=255)
                     tgcaccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgg                                                >  1:9683623‑M1/29‑67 (MQ=255)
                     cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTCGTGCTGgcgggcgcgg                                                <  1:29645040/67‑1 (MQ=255)
                       tAATGACAATCCCCGCGTAGCGGTGGCCGCGCCACGCAGCAACCTGCTGGTGCTGGCGGGCGCGGGCa                                             <  1:29913726/68‑1 (MQ=37)
                       tAATGACAAACACCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCa                                             <  1:23392651/68‑1 (MQ=255)
                         aTGACACACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:18056166/68‑1 (MQ=255)
                         aTGACAAACAGCTCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:32564169/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:18847884/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:43623377/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:28825905/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:15198366/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:26428834/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:13193519/68‑1 (MQ=255)
                         aTGACAAACAGCGCGAAGCGGGGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                           <  1:20318292/68‑1 (MQ=255)
                          tGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTg                                          >  1:50177239/1‑68 (MQ=39)
                          tGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCATCTGGGTGTGGCGGGCGCGGGCAggg                                          >  1:27814888/1‑66 (MQ=255)
                          tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTg                                          >  1:49680557/1‑68 (MQ=255)
                          tGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCGTGTGGGGCTGGCGGGCGCGGGCAGTg                                          >  1:2553994/1‑68 (MQ=255)
                           gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                         <  1:41431202/68‑1 (MQ=255)
                           gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                         <  1:14573345/68‑1 (MQ=255)
                           gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                                         <  1:8936738/68‑1 (MQ=255)
                            aCAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGt                                        >  1:39132072/1‑68 (MQ=255)
                              aaaCAGCGCGAAGCGGTGGCCGCCCCACCCAGCAACCGTCTGGTGCTGGCGGGGGCGGGCAGTGGTaa                                      >  1:37215841/1‑68 (MQ=255)
                               aaCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGGTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:30457214/1‑68 (MQ=255)
                                acacCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGa                                    <  1:7823107/68‑1 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCCCGCCGCCACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:48362462/1‑67 (MQ=38)
                                aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGGAGTGGTAAg                                     >  1:42325050/1‑67 (MQ=38)
                                aCAGCGCGAAGCGGTGGCCGCGCCCCCCAGCCACCGTGTGGGGCTGGCGGGCGCGGGCCGTGGTaat                                     >  1:1559612/1‑66 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGGTGGCGGGCGCGGGCAGTGGTTAg                                     >  1:34416354/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:21377550/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:35534502/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTATGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:38070681/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:10577055/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                                     >  1:31807703/1‑67 (MQ=255)
                                aCAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTaat                                     >  1:40971113/1‑66 (MQ=255)
                                  aGCGCGAAGCGGTGGCCGCGCCCCCCAGCAACCATCTGGTGCTGGCGGGCGCGGGGAGTGGTAAGAc                                   >  1:17573905/1‑67 (MQ=38)
                                  aGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                                   >  1:2600881/1‑67 (MQ=255)
                                    cgcgAAGCGGTGGCCGCGCCACGCCGCAACCGTCTGGGGCTGGCGGGCGCGGGCCGTGGTAAGAcgc                                 >  1:44000780/1‑67 (MQ=255)
                                     gcgcAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                               <  1:44676224/68‑1 (MQ=255)
                                     gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                                  <  1:22361132/65‑1 (MQ=255)
                                     gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                                  <  1:1178835/65‑1 (MQ=255)
                                     gcgAAGCGCTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                               <  1:47889218/68‑1 (MQ=255)
                                       gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgcg                              <  1:37495413/67‑1 (MQ=255)
                                       gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                             <  1:48908369/68‑1 (MQ=255)
                                       gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                             <  1:43231453/68‑1 (MQ=255)
                                        aaGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                            <  1:22590867/68‑1 (MQ=255)
                                         aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                            <  1:41718435/67‑1 (MQ=255)
                                         aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                            <  1:6698279/67‑1 (MQ=255)
                                          gCGGTGGCCGCGCCCCCCCGCAACCTGCTGGTGGTGGCGGGCGCGGGGAGTGGTAAGAcgcgcgg                             >  1:6117261/1‑64 (MQ=255)
                                          ccactgcCCGCGCCACGCAGCACCCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                             <  1:21025004‑M1/58‑1 (MQ=255)
                                            actgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAGACGCGCGTGCTgg                        <  1:33969807‑M1/63‑1 (MQ=255)
                                              tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt                       <  1:5781152/67‑1 (MQ=255)
                                              tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt                       <  1:19106115/67‑1 (MQ=255)
                                               ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc                     <  1:46016057/68‑1 (MQ=255)
                                                 CCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCa                    <  1:7248393‑M1/67‑1 (MQ=255)
                                                        aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCg            <  1:46995772/68‑1 (MQ=255)
                                                         cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGc           <  1:40481262/68‑1 (MQ=255)
                                                          gcagcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCg            >  1:12433190/1‑66 (MQ=255)
                                                             gcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTgg       >  1:26977514/1‑68 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTTAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:24318931/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTg        <  1:37804537/65‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTg        <  1:38434835/65‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGtt     >  1:36448025/1‑68 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:20772061/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:28255683/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:45729727/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:44406337/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:42807040/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:30963705/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:15940201/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:2386222/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:36717607/67‑1 (MQ=255)
                                                               aaCCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGt      <  1:36001664/67‑1 (MQ=255)
                                                                 ccTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTtga   >  1:45109949/1‑68 (MQ=255)
                                                                  cTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTtgat  >  1:14156974/1‑68 (MQ=255)
                                                                  cTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTtgat  >  1:39780484/1‑68 (MQ=255)
                                                                  cTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCAGGTtgat  >  1:8473425/1‑68 (MQ=255)
                                                                   |                                                                  
GTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGAT  >  CP009273/3991349‑3991482

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: