Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,287 T→C 8.5% R88R (CGT→CGC yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2870TC 8.5% 2022.2 / 118.9 626R88R (CGT→CGCyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (428/143);  new base C (38/15);  total (467/159)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

Only 100 of 633 total aligned reads displayed.

TTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCA  >  CP009273/2762220‑2762354
                                                                   |                                                                   
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                  >  1:22851757/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                  >  1:19630948/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                  >  1:18463036/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                  >  1:25685813/1‑68 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                             <  1:31675346/68‑1 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                             <  1:28293891/68‑1 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                             <  1:15140049/68‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                         >  1:5473092/1‑68 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                         >  1:21968735/1‑68 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                         >  1:38174597/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        >  1:3217552/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        <  1:7637894/68‑1 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        >  1:13134951/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        >  1:34853143/1‑68 (MQ=255)
                ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                      >  1:400852/1‑67 (MQ=255)
                   aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                   <  1:9134562/67‑1 (MQ=255)
                   aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                                  >  1:28588709/1‑68 (MQ=255)
                    cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                 <  1:7988529/68‑1 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                >  1:19328847/1‑68 (MQ=255)
                      tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                               >  1:23257494/1‑68 (MQ=255)
                      tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                               >  1:15468172/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                              >  1:17826409/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                              >  1:35400497/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                               <  1:21185745/67‑1 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                            >  1:15756071/1‑68 (MQ=255)
                          ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                           >  1:33170984/1‑68 (MQ=255)
                             cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                        <  1:27642730/68‑1 (MQ=255)
                               tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTATACAGGtt                                       >  1:25411700/1‑67 (MQ=255)
                                  tcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                         <  1:32756381/62‑1 (MQ=255)
                                   cgtTAGAATACTAATCATATTATGCTCCTACGTGGGGAGTGGCAAATTACCGGTAAACAGTTTCCGGa                                  <  1:11926348/68‑1 (MQ=255)
                                   cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                  <  1:20364579/68‑1 (MQ=255)
                                    gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                 <  1:4998228/68‑1 (MQ=255)
                                    gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                 <  1:31130885/68‑1 (MQ=255)
                                        gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                             >  1:21694550/1‑68 (MQ=255)
                                           tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                            >  1:125266/1‑66 (MQ=255)
                                            aCTAAGCTTATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                         <  1:424228/68‑1 (MQ=255)
                                            aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                            >  1:36466560/1‑65 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                         >  1:23684233/1‑67 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                         >  1:10123082/1‑67 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                        <  1:17212806/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                         >  1:11012758/1‑66 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:36057537/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:18287022/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:24884327/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:11592911/68‑1 (MQ=255)
                                               aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                       <  1:37939870/64‑1 (MQ=255)
                                               aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                       <  1:1346438/64‑1 (MQ=255)
                                               aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                      >  1:17929572/4‑68 (MQ=255)
                                               aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                      >  1:13932463/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:7230402/1‑68 (MQ=255)
                                                 acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                    >  1:35035228/2‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:21478570/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:29734959/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:11459902/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:12311248/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:23952795/1‑68 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:6107103/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:35798874/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:9624550/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:27844137/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:18100352/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:24595391/1‑67 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAATTTATCCTGCTTa                    >  1:10812979/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:8279610/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:5816887/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:19091486/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:36974872/1‑66 (MQ=255)
                                                    tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                 <  1:26729317/68‑1 (MQ=255)
                                                    tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                  >  1:12823875/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct               <  1:30205234/68‑1 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:30647426/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:27043654/1‑67 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:3744360/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:32095835/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:39019129/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:14801249/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:20786163/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:34174338/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:8744137/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTAACtctc              >  1:29352645/1‑68 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                             >  1:15902398/1‑36 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:14309690/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:32604596/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:29201857/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:6819477/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:16256160/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:16818229/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:17467332/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:34515959/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:38664988/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:38453663/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:33828871/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:24357782/1‑67 (MQ=255)
                                                         tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                 >  1:8472195/1‑63 (MQ=255)
                                                            gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt         >  1:22494353/1‑68 (MQ=255)
                                                            gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt         >  1:2606300/1‑68 (MQ=255)
                                                             ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg        >  1:4553976/1‑68 (MQ=255)
                                                               tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc      <  1:26455181/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     <  1:19892947/68‑1 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc     >  1:10014583/1‑66 (MQ=255)
                                                                   |                                                                   
TTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCA  >  CP009273/2762220‑2762354

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: