Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,281 C→T 11.5% 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,2810CT11.5% 1577.8 / 162.4 538R86R (CGC→CGTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (334/142);  new base T (55/7);  total (389/149)
Fisher's exact test for biased strand distribution p-value = 1.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 562 total aligned reads displayed.

CATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATG  >  CP009273/2762214‑2762348
                                                                   |                                                                   
   tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                   >  1:4770125/1‑67 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACTTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:32351958/68‑1 (MQ=255)
    cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                 <  1:27683362/68‑1 (MQ=255)
     tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                                <  1:30248620/68‑1 (MQ=255)
      ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                               <  1:11497423/68‑1 (MQ=255)
        gCCTGTCAGACACGGGTACCTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                             >  1:5757109/1‑68 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:22359485/1‑68 (MQ=255)
         ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                            >  1:10827154/1‑68 (MQ=255)
           tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                          <  1:2660911/68‑1 (MQ=255)
           tGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGa                                                          >  1:8810666/1‑68 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                     <  1:1305292/66‑1 (MQ=255)
                  cacGGGTACGTGACCCCCGTCGGCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                     <  1:1196793/66‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:34312468/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  <  1:26956558/68‑1 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  >  1:8145295/1‑68 (MQ=255)
                   acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                  <  1:33422042/68‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:28203190/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:6528016/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:20817694/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                             <  1:31868090/67‑1 (MQ=255)
                         aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCg                                            >  1:15306778/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATGCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:3243333/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:17598364/1‑68 (MQ=255)
                           gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                          >  1:29393144/1‑68 (MQ=255)
                             gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                        >  1:13679133/1‑68 (MQ=255)
                               cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                      >  1:3066911/1‑68 (MQ=255)
                               cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                      >  1:12372924/1‑68 (MQ=255)
                               cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                      >  1:19743859/1‑68 (MQ=255)
                                ccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACa                                     >  1:259095/1‑68 (MQ=255)
                                 cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                    >  1:31143727/1‑68 (MQ=255)
                                              gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                       >  1:16796289/1‑68 (MQ=255)
                                                 tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:1514343/1‑66 (MQ=255)
                                                 tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATcc                    >  1:21410920/1‑68 (MQ=255)
                                                  aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                      >  1:31240921/1‑65 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:18486239/1‑67 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                   >  1:534547/1‑67 (MQ=255)
                                                   cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                  <  1:937904/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:14346317/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:13355384/68‑1 (MQ=255)
                                                    tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                 <  1:22161951/68‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                 <  1:2183009/64‑1 (MQ=255)
                                                     aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                >  1:29134721/4‑68 (MQ=255)
                                                     aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                >  1:31541674/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:7819397/1‑68 (MQ=255)
                                                       gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:10034764/1‑68 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa              >  1:71282/2‑68 (MQ=255)
                                                       acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCACGAAGTTCCTGAAATTATCCTGCTTa              >  1:5137769/2‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:1781363/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:6671434/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:19468235/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:25073811/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:28674242/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:26820894/1‑68 (MQ=255)
                                                        cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa             >  1:17976275/1‑68 (MQ=255)
                                                        cATATGCTGCGCCTAGGTCGGGAGGGGCAATTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:11712770/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATGCTGCTTa              >  1:33993666/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:28941115/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:11171839/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:10979699/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:24078713/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:27807798/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:26675927/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:26360384/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:25267689/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:12836046/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:20468466/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:18678621/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:17254874/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:17119108/1‑67 (MQ=255)
                                                        cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:16690053/1‑67 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:14172434/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa              >  1:18288468/1‑66 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            <  1:27428937/68‑1 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            <  1:16374527/68‑1 (MQ=255)
                                                         atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            <  1:3295778/68‑1 (MQ=255)
                                                          tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct           <  1:23101991/68‑1 (MQ=255)
                                                          tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc            >  1:24728653/1‑67 (MQ=255)
                                                          taggctgcGCCTCCGTCGTGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAAct           <  1:10319220/65‑1 (MQ=255)
                                                            tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc          >  1:22767423/1‑67 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:26284667/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:8940982/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:34571840/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:15996758/1‑68 (MQ=255)
                                                             gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:29672668/1‑68 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGTTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa       >  1:3791757/1‑68 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActcgc        >  1:20975141/1‑65 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:25591834/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:9800992/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:14754283/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:15621734/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:20022594/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:7503084/1‑67 (MQ=255)
                                                              ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc        >  1:5566545/1‑67 (MQ=255)
                                                               tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa       <  1:33629490/67‑1 (MQ=255)
                                                               tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc      >  1:2991751/1‑68 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:23315365/68‑1 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:6094453/68‑1 (MQ=255)
                                                                gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg     <  1:4153198/68‑1 (MQ=255)
                                                                  gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt   >  1:28437337/1‑68 (MQ=255)
                                                                  gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt   >  1:8420731/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: