Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,278 G→A 14.7% 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,2780GA14.7% 1113.5 / 169.8 421L85L (CTG→CTA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (226/133);  new base A (55/7);  total (281/140)
Fisher's exact test for biased strand distribution p-value = 3.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

Only 100 of 559 total aligned reads displayed.

GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345
                                                                  |                                                                   
gCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    >  1:330859/1‑68 (MQ=255)
  cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    <  1:20748126/66‑1 (MQ=255)
     tatttGCCTGTCAGCCACGGGGAGGTGACCCCCGTCGTCGTGAGGATTCTAAGCATATGCTGCGCCTa                                                               <  1:32376936/66‑1 (MQ=25)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                >  1:5663189/1‑67 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:21915774/68‑1 (MQ=255)
       tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                             <  1:15376274/68‑1 (MQ=255)
        ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                            <  1:11497423/68‑1 (MQ=255)
        ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                            <  1:27686083/68‑1 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:25749429/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:10827154/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:598373/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:26831351/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:14430521/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:18223561/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGGCGTGAGCATACTAAGCATATTCTGCGCCTACGTCgt                                                         >  1:28236155/1‑67 (MQ=255)
            cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         <  1:26645076/67‑1 (MQ=255)
               tCAGACCCGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                     >  1:4911190/1‑68 (MQ=255)
               acAGACACGCGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAAGCTGCGCCTACCTCGGGAgg                                                     >  1:31170905/2‑68 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:3744753/68‑1 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:16760626/68‑1 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:1337079/68‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                  <  1:712007/66‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                  <  1:1305292/66‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAGGTCGGGAGGGGCaaa                                               >  1:27559602/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:16856789/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:22806186/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:28745745/68‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:28963075/68‑1 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:34312468/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:17447716/1‑68 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:2931226/1‑68 (MQ=255)
                               gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:31309916/1‑68 (MQ=255)
                                aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                     >  1:6522539/1‑67 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:16013827/1‑68 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:19743859/1‑68 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:33395152/68‑1 (MQ=255)
                                            gtgtGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:15178376/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:25997970/68‑1 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:30780246/68‑1 (MQ=255)
                                                gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                    >  1:32438826/1‑68 (MQ=255)
                                                  aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:25467685/1‑67 (MQ=255)
                                                   tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:12627807/1‑66 (MQ=255)
                                                   tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATcc                 >  1:21410920/1‑68 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGAAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:26302821/1‑65 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                >  1:18486239/1‑67 (MQ=255)
                                                     cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg               <  1:15622412/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:9280066/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:8916047/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:14772896/68‑1 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:31541674/1‑68 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:29572198/1‑68 (MQ=255)
                                                        aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt            >  1:12858648/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:10034764/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:17193476/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:32016544/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:14179275/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:6671434/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:270432/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:19468235/1‑68 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATGCTGCTTa           >  1:33993666/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:29719747/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:24078713/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:10979699/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:11171839/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:20468466/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:17849792/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:8035364/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:18678621/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:2856538/1‑67 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:4642141/68‑1 (MQ=255)
                                                            tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct        <  1:23101991/68‑1 (MQ=255)
                                                            tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         >  1:22163100/1‑67 (MQ=255)
                                                            tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         >  1:24728653/1‑67 (MQ=255)
                                                            taggctgcGCCTCCGTCGTGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAAct        <  1:10319220/65‑1 (MQ=255)
                                                             atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       <  1:13700081/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct      <  1:8373956/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:27926263/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:5406577/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:26591910/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:22431676/1‑67 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:11793163/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:18059418/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:7452361/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:34571840/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:26733988/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:12043985/1‑68 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcgc     >  1:20975141/1‑65 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:8759827/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:7815087/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:9800992/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:16428262/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:20022594/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:5291392/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:25142649/1‑67 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa    <  1:23606634/67‑1 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa    <  1:33629490/67‑1 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:1762603/1‑68 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:30048036/1‑68 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc   >  1:3030181/1‑68 (MQ=255)
                                                                  gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg  <  1:4153198/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: