Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,281 C→T 10.8% R86R (CGC→CGT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2810CT10.8% 1580.4 / 149.4 536R86R (CGC→CGTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (337/140);  new base T (48/10);  total (386/150)
Fisher's exact test for biased strand distribution p-value = 6.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 560 total aligned reads displayed.

CATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATG  >  CP009273/2762214‑2762348
                                                                   |                                                                   
   tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                   >  1:1025966/1‑67 (MQ=255)
   tCTTTGCCTGTCAGACACGGGTACGTGACCACCGTCGTCGTGAGCATACTAAGCATATGCTgcgcat                                                                   >  1:8257410/1‑65 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:36440337/68‑1 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:33342917/68‑1 (MQ=255)
      ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                               <  1:561563/68‑1 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:8624348/1‑68 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:25182252/1‑68 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:33136400/1‑68 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCCTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:22181683/1‑68 (MQ=255)
              cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                       <  1:13874118/68‑1 (MQ=255)
              cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                       <  1:32708937/68‑1 (MQ=255)
               aGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                       <  1:6558230/67‑1 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                      <  1:24603180/65‑1 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                     <  1:23266246/66‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:36203553/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  <  1:10078144/68‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  <  1:17962607/68‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAACCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:17331746/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTAGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:24001203/1‑68 (MQ=255)
                      ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                >  1:965517/1‑67 (MQ=255)
                      ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                >  1:12700037/1‑67 (MQ=255)
                       gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                               >  1:20657949/1‑67 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACTTCGGGAGGGGCAAATTAcc                                             <  1:26750128/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:29982088/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:31589269/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                            >  1:33699883/1‑68 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                            >  1:19743522/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          <  1:16905441/68‑1 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          <  1:4450918/68‑1 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:21099327/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:1496588/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          <  1:10883701/68‑1 (MQ=255)
                            tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                         >  1:21518317/1‑68 (MQ=255)
                            tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                         >  1:13190941/1‑68 (MQ=255)
                             gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:31900032/1‑68 (MQ=255)
                             gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:12997664/1‑68 (MQ=255)
                                 cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                    >  1:3229626/1‑68 (MQ=255)
                                   cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                  >  1:36980699/1‑68 (MQ=255)
                                     tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                 >  1:31167317/1‑67 (MQ=255)
                                      cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                               >  1:28300957/1‑68 (MQ=255)
                                          gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                           <  1:340390/68‑1 (MQ=255)
                                          gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                           <  1:20849463/68‑1 (MQ=255)
                                          gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                           <  1:29357111/68‑1 (MQ=255)
                                                 tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:21810085/1‑66 (MQ=255)
                                                  aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:20346590/1‑65 (MQ=255)
                                                  aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:27070261/1‑65 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:36712695/1‑67 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:18930148/1‑67 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:1580181/1‑66 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:507932/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:12321554/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:26352666/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:7039055/68‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                 <  1:10650072/64‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                >  1:25478932/4‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:19186196/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:24805616/1‑68 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa              >  1:34765993/2‑68 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa              <  1:28646105/67‑1 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa              >  1:35069360/2‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:34381639/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:22511800/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:19405137/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:35655160/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:29167153/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:27672331/1‑68 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTATCCTGCTTa              >  1:7473879/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:16032928/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:4021867/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:22258793/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:13399908/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:2349853/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACACGTTCCGGAAATTATCCTGCTTa              >  1:36031626/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGCAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtaa              >  1:3862718/1‑65 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCTAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:30666162/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:23718973/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            <  1:16403554/68‑1 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            <  1:16789092/68‑1 (MQ=255)
                                                          tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            >  1:24957783/1‑67 (MQ=255)
                                                          tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAAc            >  1:5814034/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct         <  1:27316215/68‑1 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:35539236/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:24356355/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:4834525/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:25823595/1‑67 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:938151/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:1467280/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:18261734/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:11084528/1‑68 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:15379430/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:229569/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:8428517/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:6850811/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:29571947/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:33518845/1‑67 (MQ=255)
                                                               tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc      >  1:9028214/1‑68 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:18499026/68‑1 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:14465818/68‑1 (MQ=255)
                                                                   ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg  >  1:6149790/1‑68 (MQ=255)
                                                                   ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg  >  1:20001341/1‑68 (MQ=255)
                                                                   |                                                                   
CATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATG  >  CP009273/2762214‑2762348

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: