Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,991,405 A→C 13.0% 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,4050AC13.0% 1383.3 / 67.1 545P21P (CCA→CCCuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base A (179/287);  new base C (75/4);  total (254/291)
Fisher's exact test for biased strand distribution p-value = 3.37e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 555 total aligned reads displayed.

CCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCG  >  CP009273/3991340‑3991472
                                                                 |                                                                   
 cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgc                                                                  >  1:5912308/1‑68 (MQ=255)
 cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACg                                                                   <  1:2976715/67‑1 (MQ=255)
 cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACg                                                                   <  1:4865859/67‑1 (MQ=255)
 cAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACg                                                                   <  1:5324736/67‑1 (MQ=255)
  aaTGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgca                                                                 <  1:22043243/68‑1 (MQ=255)
    tGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCgcagc                                                               >  1:10604865/1‑68 (MQ=255)
    tGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgcagc                                                               <  1:24516083/68‑1 (MQ=255)
    tGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACgcagc                                                               <  1:22948399/68‑1 (MQ=255)
      gACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCccagca                                                              >  1:29500937/1‑67 (MQ=255)
        cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                            <  1:26578957/67‑1 (MQ=255)
        cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                            <  1:28233163/67‑1 (MQ=255)
        cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                            <  1:12056986/67‑1 (MQ=255)
        cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAAc                                                            <  1:16284572/67‑1 (MQ=255)
        cGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCCCGCCGCAAcc                                                           >  1:18476211/1‑68 (MQ=255)
         gTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCt                                                          <  1:19130956/68‑1 (MQ=255)
            tCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCAGCACCCt                                                          >  1:27347761/1‑65 (MQ=255)
               tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:25846045/68‑1 (MQ=255)
               tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:11651163/68‑1 (MQ=255)
               tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:13461109/68‑1 (MQ=255)
               tACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    <  1:22193090/68‑1 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCAGCAACCATCTGGTgg                                                  >  1:17819817/1‑67 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCAACCTTCTGGt                                                    >  1:20518880/1‑66 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTgtggtgg                                                  >  1:9039207/1‑67 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:18282975/1‑66 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:29217165/1‑66 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGt                                                    >  1:16793965/1‑66 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGc                                                  >  1:24894619/1‑68 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGc                                                  >  1:6182888/1‑68 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTgg                                                  >  1:13351366/1‑67 (MQ=255)
                 ccTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCTTCTggg                                                    >  1:778525/1‑65 (MQ=255)
                    ggcgatacgatgcaccagtacgcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTCTGGTGCTg                                                <  1:13787223‑M1/29‑1 (MQ=255)
                        gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggcggg                                           <  1:9636259/68‑1 (MQ=255)
                           aGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgc                                        >  1:16550025/1‑68 (MQ=255)
                              tgcaccagtacgcgcgtcttaccactgcCCGCGCCACGCACCAACCTTCTGGTGCTGgcgggcgcgg                                      >  1:25486671‑M1/29‑67 (MQ=255)
                              cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcggg                                     <  1:29341195/68‑1 (MQ=255)
                              cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcggg                                     <  1:9337409/68‑1 (MQ=255)
                              cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGgcgggcgcggg                                     <  1:6951799/68‑1 (MQ=255)
                              cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCGGgcgggcccggg                                     >  1:24723357/1‑68 (MQ=255)
                                  atcaCAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:19583326/65‑1 (MQ=255)
                                  aTGACAACCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:23348947/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCTGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:18761160/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:26976499/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:11024774/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:23972655/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:15851304/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:7169488/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:17458874/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:21121283/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:28388273/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:18942333/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:7520597/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:3155966/68‑1 (MQ=255)
                                  aTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGt                                 <  1:19233812/68‑1 (MQ=255)
                                   tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:14446488/1‑68 (MQ=255)
                                   tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                <  1:4335053/68‑1 (MQ=255)
                                   tGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:3683755/1‑68 (MQ=255)
                                   tGACAAACAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTg                                >  1:15020840/1‑68 (MQ=255)
                                    gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTgg                               <  1:11202972/68‑1 (MQ=255)
                                     aCAAACAGCGCGAAGCGGTGGCCGCGCCACGCCGCCACCGGGGGGTGGTGGGGGGCGGGGGGggtggt                              >  1:23707385/1‑68 (MQ=255)
                                       aaaCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTaa                            >  1:21869413/1‑68 (MQ=255)
                                         gcgcgtcttaccactgcCCGCGCCACGCAGCAACCTTGTGGTGCTGGCGGGCGGGGGCAGTGGGAAg                           <  1:12879288‑M1/50‑1 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:15676604/1‑67 (MQ=39)
                                         aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCTTCTTGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:20188760/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCCCGCAGCAACCGTCTGGTGCGGGCGGGCGCGGGGAGTGGTTAg                           >  1:18060212/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCCCCCCGCCACCGTCTGGTGGTGGGGGGGGGGGGGCGTGGTTAg                           >  1:10003612/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCCCCCAGCAACCTTATGGTGCTGGCGGGCGCGGGGCGTGGGAAg                           >  1:22632705/1‑67 (MQ=37)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACGCCCCCACCGTGTGGTGGTGGCGGGCGCGGGCAGTGGTTAg                           >  1:10705158/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGGTGGCGGGCGCGGGCAGTGGTAAg                           >  1:8783090/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTggc                                              >  1:6705985/1‑48 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:21530679/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTGATGGTGCTGGCGGGCGCGGGCAGTGGTAAg                           >  1:14063589/1‑67 (MQ=255)
                                         aCAGCGCGAAGCGGTGGCCGCGCCACCCAGCAACCTTGGGGTGCTGGCCGGCGCGGGCAGTGGTAAg                           >  1:3425336/1‑67 (MQ=255)
                                          cAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAc                         >  1:7347045/1‑68 (MQ=255)
                                           aGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAgcc                         >  1:12522540/1‑65 (MQ=255)
                                              gcgATGCGGTGGCCGCGCCACGCTGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:5489956/65‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     <  1:20848117/68‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     <  1:14554881/68‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     <  1:15221622/68‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:11511777/65‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGAcg                        <  1:25042518/65‑1 (MQ=255)
                                              gcgAAGCGGTGGCCGCGCCACGCACCAACCGTCTGGTGGTGGCGGGCGCGGGCAGTGGTAAGAcgcgc                     >  1:809967/1‑68 (MQ=255)
                                                gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGt                   >  1:12629448/1‑68 (MQ=255)
                                                   gCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:19877366/68‑1 (MQ=255)
                                                     gggggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:19366210/63‑1 (MQ=255)
                                                     ggTGGCCGCGCCCCGCAGCAACCGTCTGGTGGGGGCGGGCGCGGGCAGTGGTAAGACGCGCGGACTgg              >  1:26184722/1‑68 (MQ=255)
                                                     ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACt                <  1:14222688/66‑1 (MQ=255)
                                                     ggTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTgg              <  1:8329745/68‑1 (MQ=255)
                                                     actgcCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTTAGACGCGCGTTCTgg              <  1:5771122‑M1/63‑1 (MQ=255)
                                                       tGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGt             <  1:3874086/67‑1 (MQ=255)
                                                        ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:15460322/68‑1 (MQ=255)
                                                        ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:17230760/68‑1 (MQ=255)
                                                        ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGc           <  1:20999369/68‑1 (MQ=255)
                                                         cCCGCGCCACGCAGCAACCGTCTGGTGCTGCCGGGCGCGGGCAGTGGTAAGACGCGCGGACTGgtgga          <  1:6523199‑M1/67‑3 (MQ=255)
                                                              GCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt      >  1:27897993‑M1/1‑67 (MQ=255)
                                                              GCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt      >  1:13058787‑M1/1‑67 (MQ=255)
                                                              GCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt      >  1:29017142‑M1/1‑67 (MQ=255)
                                                              GCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGt      >  1:27567898‑M1/1‑67 (MQ=255)
                                                               CCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTa     <  1:10353876‑M1/67‑1 (MQ=255)
                                                                 aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATc   >  1:25925909/1‑67 (MQ=255)
                                                                 aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCg  <  1:7696456/68‑1 (MQ=255)
                                                                 |                                                                   
CCAATGGACGTTTCTTACCTGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCG  >  CP009273/3991340‑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: