Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,312 C→G 10.2% Q97E (CAG→GAG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3120CG10.2% 2029.4 / 161.1 609Q97E (CAG→GAG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (295/252);  new base G (55/7);  total (350/259)
Fisher's exact test for biased strand distribution p-value = 4.36e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 616 total aligned reads displayed.

CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCAT  >  CP009273/2762245‑2762379
                                                                   |                                                                   
cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                     >  1:32765835/1‑68 (MQ=255)
  cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                   >  1:7870995/1‑68 (MQ=255)
   ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                                  <  1:34870453/68‑1 (MQ=255)
     gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                                                 >  1:19822110/1‑67 (MQ=255)
          cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                           >  1:11894630/1‑68 (MQ=255)
           gtgtGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                          <  1:15178376/68‑1 (MQ=255)
                  tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATcc                                                   >  1:21410920/1‑68 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTAAACAGGTTCCGGAAATTATCCt                                                  <  1:30015463/68‑1 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCg                                                             <  1:17910592/57‑1 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                     >  1:28690453/1‑65 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                  <  1:34325357/68‑1 (MQ=255)
                   aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGAAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                     >  1:26302821/1‑65 (MQ=255)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                  >  1:18486239/1‑67 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                <  1:13479558/68‑1 (MQ=255)
                      aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                               >  1:10303445/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:33968571/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:32385254/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:10034764/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:20130125/1‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                             >  1:16625383/2‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                             >  1:71282/2‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                             >  1:4017149/2‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                            >  1:19175032/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                            >  1:19571084/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                            >  1:30851418/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                            >  1:23690893/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                            >  1:1059315/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                                            >  1:15632952/1‑68 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:16921040/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:23101407/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:27745126/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:6737259/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:17119108/1‑67 (MQ=255)
                         cATATGCAGCGCCTACGTCGGGAGGCGCAAAATACCGGTAAACACGTTCCGCAAATTATACTGCATa                                             >  1:8476866/1‑67 (MQ=21)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                             >  1:29114777/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                           <  1:28447898/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                           <  1:23284781/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                           <  1:14584906/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                           <  1:16374527/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                          <  1:10850017/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                           >  1:24728653/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                         >  1:27926263/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                         >  1:626367/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                         >  1:5406577/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                         >  1:25720457/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                         >  1:2120389/1‑67 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCATGCTTAActc                                         >  1:1816039/1‑67 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:26649578/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:12185400/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:29672668/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:12680870/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:26733988/1‑68 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:22042334/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:8759827/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:11295086/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:7503084/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:25299862/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:14234876/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                       >  1:33541761/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCATAACtctc                                       >  1:13044381/1‑67 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                     >  1:3030181/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                     >  1:1762603/1‑68 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                   <  1:20961983/67‑1 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                   <  1:25082965/67‑1 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                  >  1:11542573/1‑68 (MQ=255)
                                    ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                 >  1:21688594/1‑68 (MQ=255)
                                    ccTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                 >  1:9649409/1‑68 (MQ=255)
                                     cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                 >  1:34628954/1‑67 (MQ=255)
                                     cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                 >  1:33335086/1‑67 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                               <  1:20444958/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                               <  1:26454329/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                               <  1:4451475/68‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                              <  1:28275756/68‑1 (MQ=255)
                                          tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa                           >  1:31354006/1‑68 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                          <  1:6407376/68‑1 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                          <  1:9305150/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:22744163/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:11094361/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:24184478/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:20690298/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:29345025/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                         <  1:15542069/68‑1 (MQ=255)
                                              gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGtt                        <  1:8136359/67‑1 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                       <  1:4972355/67‑1 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                      >  1:17239796/1‑68 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                      >  1:26847937/1‑68 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                     <  1:9117156/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                     <  1:26079664/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                    <  1:5897800/67‑1 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATGCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct           <  1:31499683/68‑1 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct           <  1:24903795/68‑1 (MQ=255)
                                                           ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg          <  1:27212296/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGTAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:31911466/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:22293120/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:18768511/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:30591622/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:31180691/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:524200/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:13746733/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg     <  1:14870994/68‑1 (MQ=255)
                                                                   |                                                                   
CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCAT  >  CP009273/2762245‑2762379

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: ‑