Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,278 G→A 11.6% L85L (CTG→CTA) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2780GA11.6% 1313.3 / 136.6 458L85L (CTG→CTA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (270/134);  new base A (38/15);  total (309/149)
Fisher's exact test for biased strand distribution p-value = 5.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 628 total aligned reads displayed.

GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345
                                                                  |                                                                   
 ccATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    >  1:19777306/1‑67 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                >  1:24649696/1‑67 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                               <  1:24769964/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:34898260/68‑1 (MQ=255)
          gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                          <  1:14111529/68‑1 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGGCTACGTCgg                                                         >  1:18168187/1‑68 (MQ=255)
             tGTCAGACCCGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgggc                                                       >  1:5109140/1‑67 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:28293891/68‑1 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:23085146/68‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                  <  1:9099362/66‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                >  1:21968735/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGTAGGGGCaaa                                               >  1:3903720/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:13134951/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:17976188/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:10039887/68‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:1961600/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:15762944/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:15207833/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:25154065/1‑68 (MQ=255)
                        ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                             >  1:400852/1‑67 (MQ=255)
                        ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                             >  1:17553310/1‑67 (MQ=255)
                        ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGGCCGGAGGGGCAAAtt                                             >  1:25935395/1‑67 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:38211686/67‑1 (MQ=255)
                            cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                        <  1:17226924/68‑1 (MQ=255)
                            cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                        <  1:5544763/68‑1 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:10176700/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:6274670/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:15468172/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:35400497/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      <  1:21185745/67‑1 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:38454370/1‑67 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:4243609/1‑67 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:29534834/1‑68 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:30288299/1‑68 (MQ=255)
                                   cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                 >  1:26723228/1‑68 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:22445930/68‑1 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:35191997/68‑1 (MQ=255)
                                     cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                               <  1:27642730/68‑1 (MQ=255)
                                           cgtTAGAATACTAATCATATTATGCTCCTACGTGGGGAGTGGCAAATTACCGGTAAACAGTTTCCGGa                         <  1:11926348/68‑1 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                         >  1:6837267/1‑68 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                         >  1:32141972/1‑68 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                         >  1:32967020/1‑68 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:32656537/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:8120639/68‑1 (MQ=255)
                                             tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                       >  1:22638424/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTa                    >  1:34242273/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:21694550/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:28716282/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:18940874/1‑68 (MQ=255)
                                                  aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   <  1:37953113/67‑1 (MQ=255)
                                                   tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:9355795/1‑66 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:11611314/1‑65 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:20192164/1‑67 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:31866954/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:11504983/68‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:10872406/64‑1 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:25459953/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:1358708/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa           >  1:18432790/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:8282571/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:16008889/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:21478570/1‑68 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:873338/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:12880185/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:14954896/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:11138135/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:16996040/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:17902881/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:27844137/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:27237331/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:26344559/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:23761492/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:32300869/1‑67 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:5816887/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:35773061/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:36974872/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:4748301/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:7468090/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:33889850/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:33393002/1‑67 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:3744360/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:29148810/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:36556589/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:31371538/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:20786163/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:12584660/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:19308863/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:12234533/1‑68 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:29298502/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:24357782/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:7774288/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:14428568/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:2995183/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:38664988/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:3607579/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:34515959/1‑67 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct        >  1:8472195/1‑63 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:24012834/1‑68 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:30775219/1‑68 (MQ=255)
                                                                  gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg  <  1:16307567/68‑1 (MQ=255)
                                                                  |                                                                   
GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345

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