Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,278 G→A 6.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,2780GA 6.1% 1341.6 / 55.3 410L85L (CTG→CTA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (243/142);  new base A (17/8);  total (260/150)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 576 total aligned reads displayed.

GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345
                                                                  |                                                                   
gCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    >  1:11864683/1‑68 (MQ=255)
  cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcg                                                                   <  1:10005801/67‑1 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                               <  1:10254139/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:429534/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:2836307/68‑1 (MQ=255)
       tttGCCAGTCAGACACGGGTACGTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                             <  1:20178447/68‑1 (MQ=255)
        ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                            <  1:25118929/68‑1 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTCCGTCgg                                                         >  1:13864014/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:18604197/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:22648605/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:8885145/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:10534563/1‑68 (MQ=255)
            cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCggg                                                        >  1:27638953/1‑68 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                  <  1:477922/66‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                  <  1:20786745/66‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:2959820/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:30353269/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:6189522/68‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:6825810/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:20906783/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:24514287/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:15830868/1‑68 (MQ=255)
                         gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                            >  1:9830526/1‑67 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:15101920/67‑1 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:23761846/67‑1 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:19877773/67‑1 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCggg                                       >  1:27856056/1‑67 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:16223879/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:9304331/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCAAATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       <  1:28672687/68‑1 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:14208451/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTAGGGAGGGGCAAATTACCGGTa                                      >  1:14634932/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:27296596/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:24834108/1‑68 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:24049204/1‑67 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:30899181/68‑1 (MQ=255)
                                     cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                               >  1:3629500/1‑68 (MQ=255)
                                      gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                              <  1:26719593/68‑1 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTTCCGGTAACCAGGTTCCGGa                         >  1:12944678/1‑68 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:25862134/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:3454037/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:574117/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:15546186/68‑1 (MQ=255)
                                               aGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                     >  1:17639768/1‑68 (MQ=255)
                                               aGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                     >  1:7034101/1‑68 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:3347828/1‑68 (MQ=255)
                                                  aTACTAATCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:4665783/1‑67 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:31334998/1‑65 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                <  1:31776279/68‑1 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:19025282/1‑67 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:13754102/1‑67 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:25670969/1‑67 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGTAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:20520001/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:25493962/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:545840/68‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:231771/64‑1 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:598549/1‑68 (MQ=255)
                                                        aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt            >  1:21934528/1‑68 (MQ=255)
                                                        aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt            >  1:8564377/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:8155729/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTa           >  1:14140250/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:14974770/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:29833905/1‑68 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:11209283/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:10920874/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:1658027/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:18406350/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:20290117/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:21396301/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:22925856/1‑67 (MQ=255)
                                                          cATATGATGCGCCTACGTCGGGACGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:3891466/1‑67 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaaa         >  1:31965895/1‑67 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa          <  1:22206148/67‑1 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:18081244/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:26512019/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:1605856/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct      <  1:5202421/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:16749900/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:7602667/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:11668997/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:25366935/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:27950892/1‑67 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:26252407/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:10745803/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:17823536/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:2923695/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:738687/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:13344640/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:29368616/1‑68 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAACTTACCGGTACACAGGTTCCGGAAATTATCCTGCTGAACtctc     >  1:12595038/1‑67 (MQ=25)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:6532716/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:28126751/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:17028622/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:19535656/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:23530618/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTAACtctc     >  1:19402173/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCTGAAATTATCGTGCTTAACactc     >  1:26938502/1‑67 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa    <  1:2238201/67‑1 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:3725209/1‑68 (MQ=255)
                                                                  gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg  <  1:19788299/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: