Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,405 A→C 7.7% P21P (CCA→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,4050AC 7.7% 2041.8 / 24.9 712P21P (CCA→CCCuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base A (243/398);  new base C (66/4);  total (310/402)
Fisher's exact test for biased strand distribution p-value = 2.33e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 722 total aligned reads displayed.

GCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCG  >  CP009273/3991339‑3991472
                                                                  |                                                                   
 ccAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCAc                                                                    <  1:4893016/67‑1 (MQ=255)
  cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgc                                                                  <  1:10365803/68‑1 (MQ=255)
   aaTGGACGTTTCTTACCTGCTCGACAGCCTTAATGGCAAACAGCGCGAAGCGGTGGCCGCGCCACgca                                                                 <  1:14153059/68‑1 (MQ=255)
     tGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACccagc                                                               >  1:27310657/1‑68 (MQ=255)
         cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAcc                                                           >  1:23391312/1‑68 (MQ=255)
         cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                            <  1:26086709/67‑1 (MQ=255)
          gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                          <  1:12529935/68‑1 (MQ=255)
            ttCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCtt                                                         <  1:17772029/67‑1 (MQ=255)
             tCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                          >  1:19552310/1‑65 (MQ=255)
                tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:13975236/68‑1 (MQ=255)
                tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:33117462/68‑1 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCCACCTTCGGGTGc                                                  >  1:28626533/1‑68 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTgg                                                  >  1:23887402/1‑67 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTGGGGTGc                                                  >  1:26222420/1‑68 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCAGCAACCTTCTGGTGc                                                  >  1:7531411/1‑68 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCGTCTGGTGc                                                  >  1:11802117/1‑68 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCGTCGGGt                                                    >  1:19861242/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCATCGGGt                                                    >  1:17482262/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCCACCTTCTGGTGc                                                  >  1:1810556/1‑68 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:34204949/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:10797897/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:3061223/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:32329317/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:25835328/1‑66 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCGTCTGGTgg                                                  >  1:3529868/1‑67 (MQ=255)
                  ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCACGCAGCAACCTTCTGGt                                                    >  1:24420280/1‑66 (MQ=255)
                   cTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCCACCTTCTggg                                                    >  1:22932608/1‑64 (MQ=255)
                     gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAACGTCTGGTGCggg                                               >  1:7329544/1‑68 (MQ=255)
                         gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                           <  1:5380997/68‑1 (MQ=255)
                         gACAGCCTTAACGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                           <  1:11604724/68‑1 (MQ=255)
                         gACACCCTTAATGCCATAGAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                           <  1:15078723/68‑1 (MQ=255)
                          aCAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCACGCAGCAACCTTCTGGTGCTggcgggc                                          >  1:21154077/1‑68 (MQ=255)
                           cAGGCTTAATGCCACACAGCGCGAAGCGGTGGCGGAGCCACGCAGCAACCTGCTGGTGCTggcgggcg                                         <  1:32162647/68‑1 (MQ=255)
                            aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCCACCATCTGGTGCTGgcgggcgc                                        >  1:15279029/1‑68 (MQ=255)
                            aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCTGgcgggcgc                                        >  1:21541572/1‑68 (MQ=255)
                             gCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcg                                       >  1:25934480/1‑68 (MQ=255)
                              ccTTAATGACCAACAGCGCCACCCGGGGCCCGCGCCACGCCGCAACCTTCTGGGGCTGgcgggcgcgg                                      <  1:4464861/68‑1 (MQ=255)
                              ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgg                                      >  1:20131609/1‑68 (MQ=255)
                               cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcgg                                      <  1:29052285/67‑1 (MQ=255)
                                 tAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCa                                   <  1:33682321/68‑1 (MQ=255)
                                 tAATGACAAACAGCCCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCa                                   <  1:13616516/68‑1 (MQ=255)
                                   atcaCAAACAGCGCGACGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:2922708/65‑1 (MQ=255)
                                   aagACAAACCGCGCGACGCGGAGGCCGCGCCACGCAGCAACCTGCTGGTGCTTGCGGGCGCGGGCAGt                                 <  1:3859519/66‑1 (MQ=255)
                                   aTGACAAACCGCGCGACGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:1897326/68‑1 (MQ=255)
                                   aTGACAAACCGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:26550073/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:6496022/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:31432645/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:961795/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:25427313/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:12925432/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:24732014/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:25135132/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:23103587/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:16308333/68‑1 (MQ=255)
                                   aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:1866836/68‑1 (MQ=255)
                                   aTGACAAACACCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:4409694/68‑1 (MQ=255)
                                   aTGACAAACACCGCGAAGCGGTGGCCGCCCCACGCAGCAACCTTCTGGCGCTGGCGGGCGCGGGCAGt                                 <  1:10694520/68‑1 (MQ=39)
                                    tGACAAACCGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                <  1:6275447/68‑1 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:17004263/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:21769668/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:14644866/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:16768459/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTg                                >  1:1933981/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:32817037/1‑68 (MQ=255)
                                    tGACAAACAGCGCGAAGCGGTGGCCCCCCCACGCAGCCACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:6844907/1‑68 (MQ=38)
                                     gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                               <  1:9142891/68‑1 (MQ=255)
                                          aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:29780871/1‑67 (MQ=255)
                                          aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGGGGTAAg                           >  1:5164948/1‑67 (MQ=38)
                                          aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:22383575/1‑67 (MQ=255)
                                          aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:15521203/1‑67 (MQ=255)
                                          aCAGCGCGAAGCGGTGGCCGCCCCACCCCGCAACCCTCTGGTGCAGGCGGGCGCGGGCAGTGGttac                           >  1:20370019/1‑64 (MQ=255)
                                               gcgACGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     <  1:8928682/68‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     <  1:24927502/68‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:31720862/65‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:33378277/65‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:30154628/65‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:8108739/65‑1 (MQ=255)
                                               gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGAcgcg                      <  1:26645283/67‑1 (MQ=255)
                                                   cgCGGTTGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTTGTAAGACGCGCGTa                  <  1:845836/66‑1 (MQ=39)
                                                   aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                  <  1:5607428/67‑1 (MQ=255)
                                                   aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                  <  1:1113440/67‑1 (MQ=255)
                                                   aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTa                  <  1:10025452/67‑1 (MQ=255)
                                                    gCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgcgg                   >  1:3674011/1‑64 (MQ=255)
                                                    gCGGTGGCCGCCCCCCCCAGCAACCTTCTTGTGCTGGGGGGCGCGGGGAGTGGGTAGTCGCGCGt                   >  1:3105678/1‑65 (MQ=255)
                                                    ccactgcCCGCGCCCCGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGGCGTGGTAAGACGCGCGt                   <  1:6124590‑M1/58‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:4072575/66‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:13832361/66‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:15815504/66‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:16524595/66‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:31250829/66‑1 (MQ=255)
                                                      ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:7904702/66‑1 (MQ=255)
                                                      actgcCCGCGCCCCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg              <  1:18274511‑M1/63‑1 (MQ=255)
                                                        tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt             <  1:31934953/67‑1 (MQ=255)
                                                        tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt             <  1:28225418/67‑1 (MQ=255)
                                                         ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:12015609/68‑1 (MQ=255)
                                                         ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:22655439/68‑1 (MQ=255)
                                                          gCCGCGCCACGCAGCAACCTGCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:23990630/67‑1 (MQ=255)
                                                                 CACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATc   >  1:27749215‑M1/1‑68 (MQ=255)
                                                                  aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCg  <  1:20794255/68‑1 (MQ=255)
                                                                  aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCg  <  1:14805358/68‑1 (MQ=255)
                                                                  |                                                                   
GCCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCG  >  CP009273/3991339‑3991472

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: