Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,278 G→A 13.1% 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,2780GA13.1% 1214.5 / 154.7 442L85L (CTG→CTA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (246/138);  new base A (48/10);  total (294/148)
Fisher's exact test for biased strand distribution p-value = 4.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

Only 100 of 551 total aligned reads displayed.

GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345
                                                                  |                                                                   
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                >  1:1025966/1‑67 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                               <  1:31479535/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:27449612/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:21916314/68‑1 (MQ=255)
       tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                             <  1:29197408/68‑1 (MQ=255)
       tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCCCCTACg                                                             >  1:20940317/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:15218558/1‑68 (MQ=255)
            cTGTCAGACACGGGTACGTGACCCCCGTCCTCGTGAGCATACTAAGCATATCCTGCGCCTACGTCgg                                                         <  1:34738044/67‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                   <  1:12994709/65‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                >  1:21497031/1‑68 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                >  1:16670257/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:10078144/68‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:8967265/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:23352254/1‑68 (MQ=255)
                        ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                             >  1:965517/1‑67 (MQ=255)
                        ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                             >  1:12700037/1‑67 (MQ=255)
                         gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                            >  1:28919407/1‑67 (MQ=255)
                          tACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          >  1:20827545/1‑68 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:35248502/67‑1 (MQ=255)
                            cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                        <  1:7762291/68‑1 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       <  1:10883701/68‑1 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:14388452/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:19987496/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:36864570/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:35601570/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:6096671/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:7464413/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:6828903/1‑68 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGGTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:24818130/1‑67 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:7003291/1‑67 (MQ=255)
                                   cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                 >  1:29984880/1‑68 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:20593609/68‑1 (MQ=255)
                                       tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                              >  1:28027358/1‑67 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                         <  1:25568601/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:25292721/68‑1 (MQ=255)
                                             tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCATATTACCGGTAAACAGGTTCCGGaa                        <  1:13601350/67‑1 (MQ=255)
                                                tcATACTAAGCATATGCTGCGCCTACGTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                     <  1:36388259/66‑1 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                     <  1:25854650/67‑1 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:16827318/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:33867895/1‑68 (MQ=255)
                                                  aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   <  1:36616487/67‑1 (MQ=255)
                                                   tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:4417212/1‑66 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGTGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:24138420/1‑68 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:36081081/1‑65 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:17186273/1‑67 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:13334363/1‑66 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:32255920/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:34167215/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:26802737/68‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:10650072/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:3939147/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:23548057/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt             >  1:33436514/4‑68 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt             >  1:18146492/4‑68 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:22257882/1‑68 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:30853483/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:32782061/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:27072766/1‑68 (MQ=255)
                                                         acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa           >  1:28461694/2‑68 (MQ=255)
                                                         acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa           >  1:5042458/2‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:4716601/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:33582676/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:35979772/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:14660618/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:2501108/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:18397998/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:1628436/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:8358526/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:32363425/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:244936/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTCaa          >  1:33178475/1‑68 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:14947450/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:13176436/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:37121378/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:19400578/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:8560515/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:31718595/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:22427002/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:23086310/1‑67 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:12184708/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:34992402/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:1908363/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCTTa           >  1:1963288/1‑66 (MQ=255)
                                                            aatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTCCTTAAc         >  1:29423744/2‑67 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:938151/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:11084528/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:26430725/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:23863279/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:555491/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac     >  1:5735920/1‑66 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:15945357/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:15596212/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:1772216/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:2461882/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:18852042/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:29572030/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:22851683/1‑67 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:6315253/1‑68 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:2018962/1‑68 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:21800946/1‑68 (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: