Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,281 C→T 8.9% 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,2810CT 8.9% 1878.1 / 129.9 601R86R (CGC→CGTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (410/136);  new base T (38/15);  total (450/151)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 631 total aligned reads displayed.

CATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATG  >  CP009273/2762214‑2762348
                                                                   |                                                                   
cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcgc                                                                     <  1:8293046/68‑1 (MQ=255)
cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcgc                                                                     <  1:33390075/68‑1 (MQ=255)
   tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                   >  1:24649696/1‑67 (MQ=255)
   tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                                  <  1:24769964/68‑1 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:25480666/68‑1 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:34898260/68‑1 (MQ=255)
        gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                             <  1:14111529/68‑1 (MQ=255)
              cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                       <  1:2929307/68‑1 (MQ=255)
              cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                        <  1:19767802/67‑1 (MQ=255)
              cAGACACGGGTACGTGACCCCCGTCGGCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                        <  1:4610347/67‑1 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                     <  1:8791617/66‑1 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                   >  1:21968735/1‑68 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                   >  1:38174597/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  <  1:10039887/68‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:25154065/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:17976188/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:15207833/1‑68 (MQ=255)
                      ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                >  1:400852/1‑67 (MQ=255)
                      ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                >  1:17553310/1‑67 (MQ=255)
                      ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGGCCGGAGGGGCAAAtt                                                >  1:25935395/1‑67 (MQ=255)
                          cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                           <  1:17226924/68‑1 (MQ=255)
                          cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                           <  1:5544763/68‑1 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:24386948/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:26310230/1‑68 (MQ=255)
                           gTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          <  1:2779193/68‑1 (MQ=255)
                            tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                         >  1:15468172/1‑68 (MQ=255)
                            tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                         >  1:23257494/1‑68 (MQ=255)
                             gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                         <  1:21185745/67‑1 (MQ=255)
                              aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:18160099/1‑67 (MQ=255)
                              aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:38454370/1‑67 (MQ=255)
                              aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:4243609/1‑67 (MQ=255)
                               cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                      >  1:30288299/1‑68 (MQ=255)
                               cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                      >  1:29534834/1‑68 (MQ=255)
                                     tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                 >  1:34179812/1‑67 (MQ=255)
                                     tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                 >  1:13170270/1‑67 (MQ=255)
                                         cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                            >  1:32141972/1‑68 (MQ=255)
                                         cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                            >  1:26559188/1‑68 (MQ=255)
                                         cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                            >  1:6837267/1‑68 (MQ=255)
                                         cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                            >  1:32967020/1‑68 (MQ=255)
                                          gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                           <  1:32656537/68‑1 (MQ=255)
                                           tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                          >  1:22638424/1‑68 (MQ=255)
                                              gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                       >  1:21694550/1‑68 (MQ=255)
                                              gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                       >  1:18940874/1‑68 (MQ=255)
                                                aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      <  1:17888413/67‑1 (MQ=255)
                                                  aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:11611314/1‑65 (MQ=255)
                                                  aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:8155601/1‑65 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:20192164/1‑67 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:23684233/1‑67 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTAGCCt                   >  1:27040984/1‑67 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:11012758/1‑66 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:31866954/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:11504983/68‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                 <  1:37939870/64‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                 <  1:9204068/64‑1 (MQ=255)
                                                      aGCATATGCTTCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt               >  1:12386781/1‑68 (MQ=255)
                                                      aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt               >  1:15795880/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:9573318/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:28972316/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa              >  1:18432790/1‑68 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa              >  1:35035228/2‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:21478570/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:23952795/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:12154951/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:16008889/1‑68 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:8481953/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:32300869/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:12947739/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:16996040/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:37543677/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:19604682/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:28434766/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:37052456/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:14524379/1‑67 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAATTTATCCTGCTTa              >  1:10812979/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:35773061/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:5816887/1‑66 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:16280174/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:7468090/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:33856206/1‑67 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:31321880/1‑67 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:6224096/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:13606898/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:14973123/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:10161998/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:39019129/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:36556589/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:34505055/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:20786163/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:28109961/1‑68 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:38664988/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:3607579/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:35363828/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:7774288/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:27376120/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:5031035/1‑67 (MQ=255)
                                                               tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc      >  1:30775219/1‑68 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:19227813/68‑1 (MQ=255)
                                                                  gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa    <  1:12653297/67‑1 (MQ=255)
                                                                  gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt   >  1:22494353/1‑68 (MQ=255)
                                                                   ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg  >  1:29830575/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: