Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,311 A→G 10.4% K96K (AAA→AAG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3110AG10.4% 1829.5 / 134.0 581K96K (AAA→AAG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (296/222);  new base G (53/7);  total (352/229)
Fisher's exact test for biased strand distribution p-value = 1.08e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 616 total aligned reads displayed.

CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378
                                                                  |                                                                   
cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                                     <  1:3242139/67‑1 (MQ=255)
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)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:31355089/1‑67 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:13479558/68‑1 (MQ=255)
                      aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                              >  1:31541674/1‑68 (MQ=255)
                      aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                              >  1:10303445/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:10034764/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:20130125/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:33968571/1‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                            >  1:71282/2‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                            >  1:4017149/2‑68 (MQ=255)
                        acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                            >  1:16625383/2‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:1059315/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:32016544/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:23690893/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:19175032/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:19571084/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:30905236/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGATCCTGAAATTATCCTGCTTaa                                           >  1:15632952/1‑68 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:6737259/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:27745126/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:23101407/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:16921040/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:16374527/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:23284781/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:28447898/68‑1 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:14584906/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:29672668/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:12185400/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:26649578/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:26733988/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:12680870/1‑68 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:25299862/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:31203378/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:22042334/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:33541761/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:14234876/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:7503084/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:11295086/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:8759827/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:4451475/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:26454329/68‑1 (MQ=255)
                                      tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                              <  1:20444958/68‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             <  1:28275756/68‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             <  1:30599487/68‑1 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         <  1:9305150/68‑1 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         <  1:6407376/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:29345025/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:22744163/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:26847937/1‑68 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                     >  1:17239796/1‑68 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:26079664/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:9117156/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                   <  1:5897800/67‑1 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct          <  1:24903795/68‑1 (MQ=255)
                                                           ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg         <  1:27212296/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:524200/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:13746733/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:18768511/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:14870994/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:22293120/68‑1 (MQ=255)
                                                                  |                                                                   
CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378

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