Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,312 C→G 8.4% Q97E (CAG→GAG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3120CG 8.4% 2174.2 / 134.7 625Q97E (CAG→GAG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (374/196);  new base G (38/15);  total (414/211)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 632 total aligned reads displayed.

CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378
                                                                   |                                                                  
cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                    >  1:15756071/1‑68 (MQ=255)
 ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                                                   >  1:33170984/1‑68 (MQ=255)
   ccGGCGTCGTGTGCATACTAAGCATATGCTGCGCCTACGTCGGTAGGGGCAAATTACCGGTAAACAgg                                                                 <  1:20202122/68‑1 (MQ=25)
      tcgtcgtGAGCATACTCAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGGAAACAGGtt                                                               >  1:35018895/1‑67 (MQ=255)
      tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                               >  1:34179812/1‑67 (MQ=255)
         tcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                                 <  1:32805014/62‑1 (MQ=255)
          cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                          >  1:26559188/1‑68 (MQ=255)
          cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                          >  1:34836713/1‑68 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:31917439/68‑1 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                      <  1:17083761/67‑1 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:28536765/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:4846397/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:13225213/1‑68 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCCGAAAtt                                                      <  1:21490720/67‑1 (MQ=255)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:35972012/1‑67 (MQ=255)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:23684233/1‑67 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:11504983/68‑1 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:24884327/68‑1 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:25412480/68‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:12469876/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:38657030/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:9204068/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:35663952/64‑1 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:7061655/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:16904944/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:36395695/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:23056633/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:23952795/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:27032757/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:29734959/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:37400364/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:16008889/1‑68 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:19604682/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:11138135/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:24595391/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:17850023/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:15092909/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:16312676/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCAGGAAATTATCATGCTTa                                            >  1:38893700/1‑67 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:5690446/1‑66 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTTTCCTGCTTAAc                                          >  1:35352536/1‑67 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                         <  1:21014522/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          >  1:20691830/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:11986069/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:15396334/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:32420721/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:18247064/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:31321880/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:36767365/1‑67 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:38355018/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:32142443/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:10161998/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:9627290/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:14801249/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:28109961/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:10907848/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:1746369/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:13606898/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:19308863/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:21254811/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:21703356/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:12786450/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTAACtctc                                      >  1:11729119/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:29352645/1‑68 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:3261895/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:39489898/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:34515959/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:7693984/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:18983887/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:14309690/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:8369284/1‑67 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:33794281/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:2782278/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:30775219/1‑68 (MQ=255)
                                 gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                   <  1:2255698/68‑1 (MQ=255)
                                 gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                   <  1:17964910/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:5937421/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:6493028/68‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             <  1:19892947/68‑1 (MQ=255)
                                        cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                             >  1:10014583/1‑66 (MQ=255)
                                        cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             >  1:4152931/1‑67 (MQ=255)
                                         gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc                           >  1:29936908/1‑68 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         <  1:30482580/68‑1 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         <  1:13000159/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:8525858/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:7428150/68‑1 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                      <  1:24371422/67‑1 (MQ=255)
                                                ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    >  1:22279399/1‑68 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:34024234/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:8228079/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:27224908/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:26230630/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:4467520/66‑1 (MQ=255)
                                                     aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg               <  1:8691163/68‑1 (MQ=255)
                                                        ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt             >  1:18524223/1‑67 (MQ=255)
                                                            cGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctgc        <  1:37879755/68‑1 (MQ=255)
                                                             ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc       >  1:2915426/1‑68 (MQ=255)
                                                             ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc       <  1:36510943/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:25763108/68‑1 (MQ=255)
                                                                 aaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGc   >  1:38043915/1‑68 (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: ‑