Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,288 C→A 10.5% R89R (CGG→AGG)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2880CA10.5% 1904.8 / 164.7 584R89R (CGG→AGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (354/166);  new base A (55/7);  total (411/173)
Fisher's exact test for biased strand distribution p-value = 5.96e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

Only 100 of 591 total aligned reads displayed.

GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355
                                                                  |                                                                   
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:14430521/1‑68 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:3744753/68‑1 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:10951462/68‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGGCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                            <  1:1196793/66‑1 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:24166278/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:8145295/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:34312468/1‑68 (MQ=255)
               gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                                      >  1:5884050/1‑67 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                    <  1:6528016/67‑1 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                    <  1:15215867/67‑1 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                    <  1:21455720/67‑1 (MQ=255)
                   gTGACCCCCGTCGTCGTGAGCATGCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 >  1:3243333/1‑68 (MQ=255)
                   gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 >  1:32079546/1‑68 (MQ=255)
                   gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 >  1:14776168/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:6746867/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:25265406/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:29638909/1‑68 (MQ=255)
                       cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                             >  1:19743859/1‑68 (MQ=255)
                        ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                            >  1:259095/1‑68 (MQ=255)
                         cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                           >  1:7870995/1‑68 (MQ=255)
                          ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                          <  1:33395152/68‑1 (MQ=255)
                             tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                        >  1:17210782/1‑67 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:17726481/68‑1 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAttt                              >  1:8820723/1‑67 (MQ=255)
                                        aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             <  1:20947479/67‑1 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCg                                     <  1:17910592/57‑1 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:534547/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:19373950/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:18486239/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                         <  1:937904/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:31443047/1‑66 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:13479558/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:22161951/68‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                       >  1:15552139/4‑68 (MQ=255)
                                              aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGcat                      >  1:23163474/1‑66 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:11374709/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:33968571/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:10034764/1‑68 (MQ=255)
                                               acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                     >  1:8947264/2‑68 (MQ=255)
                                               acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                     >  1:16625383/2‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:26820894/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:31151976/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:28426930/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:28674242/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:29343515/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:7544965/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:23690893/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:13127551/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:14179275/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:16362126/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:17000598/1‑68 (MQ=255)
                                                cATATGCTGCGCCTAGGTCGGGAGGGGCAATTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:11712770/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26360384/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26675927/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:27745126/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:18678621/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:7206717/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCCGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:15945767/1‑67 (MQ=255)
                                                cATATGCAGCGCCTACGTCGGGAGGCGCAAAATACCGGTAAACACGTTCCGCAAATTATACTGCATa                     >  1:8476866/1‑67 (MQ=21)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:18288468/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:16794792/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:29114777/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:21226390/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:22592360/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:27428937/68‑1 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:4642141/68‑1 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:3295778/68‑1 (MQ=255)
                                                  tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                  <  1:11133528/68‑1 (MQ=255)
                                                   atGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 <  1:26937690/68‑1 (MQ=255)
                                                   atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 <  1:13700081/68‑1 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGTTTCAGGAAATTATCCTGCTTAActc                 >  1:20404366/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:626367/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:2660787/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:22859665/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:5566545/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:2628323/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:9800992/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:20022594/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:17506257/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:12849174/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:25591834/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:12499770/1‑67 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa              <  1:33629490/67‑1 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc             >  1:30048036/1‑68 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc             >  1:3030181/1‑68 (MQ=255)
                                                        gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg            <  1:4153198/68‑1 (MQ=255)
                                                          gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAACTATCCTGCTTAACTCTCACGa           <  1:24937687/67‑1 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg        <  1:3073749/68‑1 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:34628954/1‑67 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         <  1:31252045/67‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:28257133/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:10777752/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:27899363/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:10311697/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:20527154/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:21921880/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:3185012/68‑1 (MQ=255)
                                                                 tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa   >  1:24646520/1‑68 (MQ=255)
                                                                 tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa   >  1:506796/1‑68 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:28915567/68‑1 (MQ=255)
                                                                  |                                                                   
GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355

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: