Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,311 A→G 8.8% K96K (AAA→AAG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3110AG 8.8% 1965.2 / 116.4 603K96K (AAA→AAG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (374/175);  new base G (38/15);  total (413/190)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 634 total aligned reads displayed.

CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378
                                                                  |                                                                   
cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                    >  1:15756071/1‑68 (MQ=255)
   ccGGCGTCGTGTGCATACTAAGCATATGCTGCGCCTACGTCGGTAGGGGCAAATTACCGGTAAACAgg                                                                 <  1:20202122/68‑1 (MQ=25)
      tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTATACAGGtt                                                               >  1:25411700/1‑67 (MQ=255)
       cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                                             >  1:15492339/1‑68 (MQ=255)
         tcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                                 <  1:32805014/62‑1 (MQ=255)
          cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                          >  1:34836713/1‑68 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:31917439/68‑1 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:3836424/68‑1 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTa                                                     >  1:34242273/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:28536765/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:21694550/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:13225213/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:4846397/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCCGAAAtt                                                      <  1:21490720/67‑1 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                    >  1:8155601/1‑65 (MQ=255)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:23684233/1‑67 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:24884327/68‑1 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:11504983/68‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:38657030/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:12469876/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:9204068/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:35663952/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                              >  1:17929572/4‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:36395695/1‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                            >  1:35035228/2‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:23952795/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:27032757/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:12154951/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:29734959/1‑68 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:19604682/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:16312676/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:15092909/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:38131500/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:24595391/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCAGGAAATTATCATGCTTa                                            >  1:38893700/1‑67 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:36974872/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:1597119/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:14419037/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTTTCCTGCTTAAc                                          >  1:35352536/1‑67 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          >  1:20691830/1‑67 (MQ=255)
                            atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        <  1:26505972/68‑1 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:18247064/1‑67 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:31371538/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:32142443/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:10161998/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:28109961/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:3744360/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:9627290/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:10907848/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:14801249/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:13606898/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:1746369/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:12786450/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:18458888/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:19308863/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:20890270/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:21254811/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTAACtctc                                      >  1:11729119/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:29352645/1‑68 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:16818229/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:34515959/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:39489898/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:3261895/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:7693984/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:2995183/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:29201857/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:24357782/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:8369284/1‑67 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:33794281/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:30775219/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:2782278/1‑68 (MQ=255)
                                 gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                   <  1:2255698/68‑1 (MQ=255)
                                    ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                >  1:18941359/1‑68 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:6975211/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:5937421/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:11171146/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:6493028/68‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             <  1:36052752/68‑1 (MQ=255)
                                        cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                             >  1:10014583/1‑66 (MQ=255)
                                          gcGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAGCCTGCTTAACTCTCACGAGGGCACCa                          >  1:32434524/2‑68 (MQ=21)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         <  1:30482580/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:17004663/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:8610570/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:8525858/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:7428150/68‑1 (MQ=255)
                                                ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    >  1:19871886/1‑68 (MQ=255)
                                                ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    >  1:22279399/1‑68 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:26230630/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:27224908/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:5390210/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:34024234/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:4467520/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                   <  1:177286/67‑1 (MQ=255)
                                                       aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt             <  1:4147321/68‑1 (MQ=255)
                                                            cGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctgc        <  1:37879755/68‑1 (MQ=255)
                                                             ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc       >  1:33295246/1‑68 (MQ=255)
                                                             ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc       <  1:36510943/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:23049513/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:25763108/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:13007525/68‑1 (MQ=255)
                                                                  |                                                                   
CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378

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: ‑