Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,287 T→C 10.8% 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,2870TC10.8% 1765.5 / 148.7 577R88R (CGT→CGCyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (347/166);  new base C (55/7);  total (404/173)
Fisher's exact test for biased strand distribution p-value = 3.69e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 590 total aligned reads displayed.

TTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCA  >  CP009273/2762220‑2762354
                                                                   |                                                                   
ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                     <  1:7136961/68‑1 (MQ=255)
ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                     <  1:27686083/68‑1 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGGCGTGAGCATACTAAGCATATTCTGCGCCTACGTCgt                                                                  >  1:28236155/1‑67 (MQ=255)
     tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                                >  1:8810666/1‑68 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                             <  1:16760626/68‑1 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                             <  1:10951462/68‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                           <  1:1305292/66‑1 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        >  1:34312468/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                        >  1:14159798/1‑68 (MQ=255)
                   aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                   <  1:31792233/67‑1 (MQ=255)
                   aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                   <  1:20817694/67‑1 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                >  1:17447716/1‑68 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                >  1:17880873/1‑68 (MQ=255)
                      tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                               >  1:25265406/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                              >  1:29638909/1‑68 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                            >  1:32765835/1‑68 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                            >  1:5439384/1‑68 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                            >  1:12372924/1‑68 (MQ=255)
                          ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                           >  1:259095/1‑68 (MQ=255)
                                    gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                 <  1:4106869/68‑1 (MQ=255)
                                    gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                 <  1:28575796/68‑1 (MQ=255)
                                           tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATcc                          >  1:21410920/1‑68 (MQ=255)
                                            aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCt                         <  1:30015463/68‑1 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                        <  1:937904/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:16098467/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:14346317/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:13355384/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:31424205/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                       <  1:185853/68‑1 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:10034764/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:7819397/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:6589632/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:33968571/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:11374709/1‑68 (MQ=255)
                                                 acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                    >  1:27887648/2‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:6204669/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:8126553/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:26820894/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:17193476/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:29343515/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:31151976/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:11915352/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                   >  1:19175032/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                   >  1:15632952/1‑68 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:23101407/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:20468466/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:16921040/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:33366791/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:26675927/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:26360384/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:32431648/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:18470683/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:12836046/1‑67 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:29114777/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                    >  1:15175767/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                  <  1:28884586/68‑1 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                  <  1:27428937/68‑1 (MQ=255)
                                                    tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                 <  1:11133528/68‑1 (MQ=255)
                                                     atGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                <  1:26937690/68‑1 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGTGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct               <  1:11613902/68‑1 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct               <  1:15360727/68‑1 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:2660787/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:1655718/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:22767423/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                >  1:22431676/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCATGCTTAActc                >  1:1816039/1‑67 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTGCGGAAATTATCCTGCTTAACtctc              >  1:19667873/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct               <  1:24024384/67‑1 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcac              >  1:17707702/1‑66 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:26284667/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:585675/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:8940982/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:2213418/1‑68 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcgc              >  1:20975141/1‑65 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:9800992/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:7503084/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:5063764/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:5291392/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:25591834/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:14754283/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc              >  1:33541761/1‑67 (MQ=255)
                                                         tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa             <  1:19845283/67‑1 (MQ=255)
                                                         tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc            >  1:3030181/1‑68 (MQ=255)
                                                          gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg           <  1:23315365/68‑1 (MQ=255)
                                                            gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt         >  1:8420731/1‑68 (MQ=255)
                                                            gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAACTATCCTGCTTAACTCTCACGa          <  1:24937687/67‑1 (MQ=255)
                                                             ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg        >  1:21688594/1‑68 (MQ=255)
                                                              cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg       <  1:3073749/68‑1 (MQ=255)
                                                              cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg        >  1:34628954/1‑67 (MQ=255)
                                                              cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg        <  1:20158892/67‑1 (MQ=255)
                                                              cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg        <  1:31252045/67‑1 (MQ=255)
                                                               tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc      <  1:6909918/68‑1 (MQ=255)
                                                               tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc      <  1:4451475/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     <  1:10311697/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     <  1:10777752/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     <  1:3535485/68‑1 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     >  1:32074907/1‑67 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     >  1:13672812/1‑67 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa     >  1:2838349/1‑67 (MQ=255)
                                                                   tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa  >  1:24646520/1‑68 (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: