Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,276 C→T 14.6% L85L (CTG→TTG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2760CT14.6% 1142.1 / 169.1 425L85L (CTG→TTG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (226/137);  new base T (55/7);  total (281/144)
Fisher's exact test for biased strand distribution p-value = 1.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

Only 100 of 545 total aligned reads displayed.

CAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCA  >  CP009273/2762210‑2762343
                                                                  |                                                                   
cAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgct                                                                    <  1:17841463/68‑1 (MQ=255)
  gCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                  >  1:330859/1‑68 (MQ=255)
    cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                  <  1:20748126/66‑1 (MQ=255)
    cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                  <  1:24081337/66‑1 (MQ=255)
     aTTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGcc                                                               <  1:18214708/68‑1 (MQ=255)
        cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                            <  1:21915774/68‑1 (MQ=255)
         tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                           <  1:15376274/68‑1 (MQ=255)
            gCCTGTCAGACACGGGTACCTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                        >  1:5757109/1‑68 (MQ=255)
             ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                       >  1:25749429/1‑68 (MQ=255)
             ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                       >  1:14430521/1‑68 (MQ=255)
             ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                       >  1:32497407/1‑68 (MQ=255)
             ccTGTCAGACACGGGTACGTGACACCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                       >  1:34036685/1‑68 (MQ=255)
               tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                     >  1:8810666/1‑68 (MQ=255)
                  cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                  <  1:16760626/68‑1 (MQ=255)
                      cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                <  1:1305292/66‑1 (MQ=255)
                       acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                             <  1:28963075/68‑1 (MQ=255)
                       acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                             <  1:33422042/68‑1 (MQ=255)
                       acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                             >  1:16856789/1‑68 (MQ=255)
                             aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                        <  1:15215867/67‑1 (MQ=255)
                             aCGTGACCCCCGTCGTCGTCAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                       >  1:26134181/1‑68 (MQ=255)
                               gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                     >  1:32079546/1‑68 (MQ=255)
                               gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                     >  1:28776015/1‑68 (MQ=255)
                               gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                     >  1:2931226/1‑68 (MQ=255)
                                 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                   <  1:27934808/68‑1 (MQ=255)
                                 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                   >  1:13679133/1‑68 (MQ=255)
                                 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                   >  1:3363042/1‑68 (MQ=255)
                                   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                 >  1:16013827/1‑68 (MQ=255)
                                   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                 >  1:5439384/1‑68 (MQ=255)
                                   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                 >  1:19743859/1‑68 (MQ=255)
                                     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                               >  1:6076309/1‑68 (MQ=255)
                                              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                      <  1:28575796/68‑1 (MQ=255)
                                              gtGAGCATACAAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                      <  1:11604433/68‑1 (MQ=255)
                                                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                  >  1:4775688/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)
                                                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                 >  1:30180890/1‑65 (MQ=255)
                                                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt              <  1:34325357/68‑1 (MQ=255)
                                                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt              >  1:31355089/1‑67 (MQ=255)
                                                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt              >  1:534547/1‑67 (MQ=255)
                                                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt              >  1:19373950/1‑67 (MQ=255)
                                                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg             <  1:937904/68‑1 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt              >  1:10800885/1‑66 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc            <  1:8916047/68‑1 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc            <  1:13355384/68‑1 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc            <  1:14772896/68‑1 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc            <  1:22161951/68‑1 (MQ=255)
                                                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc            <  1:26856458/68‑1 (MQ=255)
                                                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc            <  1:11031599/64‑1 (MQ=255)
                                                         aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt           >  1:29572198/1‑68 (MQ=255)
                                                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt          >  1:12858648/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:6589632/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:32385254/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:27074642/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:10034764/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:34808626/1‑68 (MQ=255)
                                                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:7819397/1‑68 (MQ=255)
                                                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa         >  1:71282/2‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGGAATTATCCTGCTTaa        >  1:24820268/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:17193476/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:4188290/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:30905236/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:8126553/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:23690893/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa        >  1:31771890/1‑68 (MQ=255)
                                                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa        >  1:15632952/1‑68 (MQ=255)
                                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:11171839/1‑67 (MQ=255)
                                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:29719747/1‑67 (MQ=255)
                                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:27745126/1‑67 (MQ=255)
                                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:6737259/1‑67 (MQ=255)
                                                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:20468466/1‑67 (MQ=255)
                                                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa         >  1:14172434/1‑66 (MQ=255)
                                                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct      <  1:23101991/68‑1 (MQ=255)
                                                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc       >  1:28292367/1‑67 (MQ=255)
                                                              taggctgcGCCTCCGTCGTGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAAct      <  1:10319220/65‑1 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct    <  1:8373956/68‑1 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:31237878/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:26591910/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:22767423/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:22431676/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:18477752/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:17310677/1‑67 (MQ=255)
                                                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc     >  1:1655718/1‑67 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcac   >  1:17707702/1‑66 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:26345048/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:30667107/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:26733988/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:26649578/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:18617881/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:18059418/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:12043985/1‑68 (MQ=255)
                                                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:11793163/1‑68 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGTTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa  >  1:3791757/1‑68 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcgc   >  1:20975141/1‑65 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:5063764/1‑67 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:27586736/1‑67 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:25142649/1‑67 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:7503084/1‑67 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:20022594/1‑67 (MQ=255)
                                                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc   >  1:9800992/1‑67 (MQ=255)
                                                                  |                                                                   
CAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCA  >  CP009273/2762210‑2762343

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: