Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,288 C→A 10.1% 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.1% 1842.4 / 146.4 561R89R (CGG→AGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (352/149);  new base A (48/10);  total (401/160)
Fisher's exact test for biased strand distribution p-value = 4.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 571 total aligned reads displayed.

GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355
                                                                  |                                                                   
gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                                    >  1:7630686/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:32908015/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:2462930/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:10994489/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:7096778/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:8624348/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGACGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:25448381/1‑68 (MQ=255)
  cTGTCAGACACGGGTACGTGACCCCCGTCCTCGTGAGCATACTAAGCATATCCTGCGCCTACGTCgg                                                                   <  1:34738044/67‑1 (MQ=255)
   tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                                 >  1:29731596/1‑68 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:22375149/68‑1 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:21130112/68‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                             <  1:12994709/65‑1 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCata                                                         >  1:26600750/1‑66 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:35234467/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:8967265/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:23352254/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAACCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:17331746/1‑68 (MQ=255)
              ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                       >  1:965517/1‑67 (MQ=255)
               gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                                      >  1:20657949/1‑67 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                                   >  1:19743522/1‑68 (MQ=255)
                   gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 >  1:1496588/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:13190941/1‑68 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:35601570/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:6554061/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:6096671/1‑68 (MQ=255)
                       cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                             >  1:31848446/1‑68 (MQ=255)
                           cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                         >  1:36980699/1‑68 (MQ=255)
                           cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGCAGGCGCAAATTACCGGTTAACAGGt                                         >  1:17967374/1‑68 (MQ=25)
                              cgtcgtGAGCATACTAAGCATATGCGGAGCCTACGTCGGGAGGGGGAAATTACCGGTAAACAGGTTcc                                      >  1:34260853/1‑68 (MQ=25)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:29357111/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:13737249/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:5026265/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:20849463/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:3814604/68‑1 (MQ=255)
                                   tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGaaa                                 >  1:19433918/1‑68 (MQ=255)
                                      tcATACTAAGCATATGCTGCGCCTACGTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                               <  1:36388259/66‑1 (MQ=255)
                                      tcATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                               <  1:12559011/66‑1 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                              >  1:29090370/1‑68 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                              >  1:18273587/1‑68 (MQ=255)
                                        aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACACGTTCCGGAAATTa                              <  1:36143109/66‑1 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:21810085/1‑66 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:4417212/1‑66 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:16042058/1‑66 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:36712695/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:533892/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:18930148/1‑67 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:32462362/1‑66 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:32255920/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:12321554/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:18519850/68‑1 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:20018192/68‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:23775217/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:3939147/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:23548057/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:22596190/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                       >  1:33436514/4‑68 (MQ=255)
                                              aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                      >  1:33607625/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:2815601/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:19186196/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26891669/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:35979772/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:27593488/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:20438081/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:4716601/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:8358526/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:10807458/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:25038051/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:24846290/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:10358584/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:244936/1‑68 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:25865058/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:13399908/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:30470906/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26310310/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:16989675/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:23086310/1‑67 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:34992402/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:25617469/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:21701616/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                   <  1:16403554/68‑1 (MQ=255)
                                                  tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                  <  1:30890128/68‑1 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:11547232/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:1467280/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:15458161/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:15816364/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:938151/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:29538897/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:31585970/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac               >  1:5735920/1‑66 (MQ=255)
                                                     gctgcGACTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATAATCCTGCTTAACtctc               >  1:16794994/1‑68 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:28612384/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:6850811/1‑67 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc             >  1:24008604/1‑68 (MQ=255)
                                                        gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg            <  1:22203588/68‑1 (MQ=255)
                                                        gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg            <  1:14465818/68‑1 (MQ=255)
                                                          gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt          >  1:27087110/1‑68 (MQ=255)
                                                          gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt          >  1:33196554/1‑68 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg        <  1:33928722/68‑1 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:10176499/68‑1 (MQ=255)
                                                              aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa      <  1:1003428/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: