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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2870TC10.4% 1720.2 / 138.1 560R88R (CGT→CGCyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (352/149);  new base C (48/10);  total (401/159)
Fisher's exact test for biased strand distribution p-value = 4.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

Only 100 of 563 total aligned reads displayed.

TTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACC  >  CP009273/2762220‑2762353
                                                                   |                                                                  
  gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                                  >  1:7630686/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                 >  1:26143098/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                 >  1:10994489/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                 >  1:33620384/1‑68 (MQ=255)
   ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                 >  1:17740624/1‑68 (MQ=255)
    cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                 <  1:27901665/67‑1 (MQ=255)
       tCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                             >  1:30980424/1‑68 (MQ=255)
        cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                            <  1:32708937/68‑1 (MQ=255)
         aGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                            <  1:23949670/67‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                           <  1:24603180/65‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgggg                                                           <  1:12994709/65‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                          <  1:23266246/66‑1 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                        >  1:21497031/1‑68 (MQ=255)
            cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                        >  1:16670257/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTGGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                       >  1:35253383/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                       >  1:23630164/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                       >  1:36203553/1‑68 (MQ=255)
             acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                       >  1:8967265/1‑68 (MQ=255)
                ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGTGAGGGGCAAAtt                                                     >  1:5388421/1‑67 (MQ=255)
                ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                     >  1:965517/1‑67 (MQ=255)
                   aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                  <  1:31589269/67‑1 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                               >  1:19987496/1‑68 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                               >  1:14388452/1‑68 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                               <  1:4450918/68‑1 (MQ=255)
                     gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCGACGTCGGGAGGGGCAAATTACCGGt                                               >  1:24989719/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                             >  1:6828903/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                             >  1:6554061/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                             >  1:6096671/1‑68 (MQ=255)
                       gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                             >  1:20164607/1‑68 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                           >  1:31848446/1‑68 (MQ=255)
                         cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                           >  1:17279396/1‑68 (MQ=255)
                             cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                       >  1:33353133/1‑68 (MQ=255)
                             cgtcgtcgtGAGCATACTAAGCATATGCTGCCCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                       <  1:32973487/68‑1 (MQ=255)
                                    gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                <  1:340390/68‑1 (MQ=255)
                                     tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCATATTACCGGTAAACAGGTTCCGGaa                                <  1:13601350/67‑1 (MQ=255)
                                        tcATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                             <  1:12559011/66‑1 (MQ=255)
                                        gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                             <  1:21897116/67‑1 (MQ=255)
                                        gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGATCCGGAAATTa                            >  1:19106866/1‑68 (MQ=255)
                                           tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                           >  1:29687609/1‑66 (MQ=255)
                                            aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                           >  1:32282536/1‑65 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                        >  1:20804241/1‑67 (MQ=255)
                                             cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                        >  1:20552107/1‑67 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                        >  1:1580181/1‑66 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                      <  1:8587115/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                      <  1:27071715/68‑1 (MQ=255)
                                              tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                      <  1:26802737/68‑1 (MQ=255)
                                               aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                      <  1:503113/64‑1 (MQ=255)
                                               aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                     >  1:25478932/4‑68 (MQ=255)
                                               aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                     >  1:22257882/1‑68 (MQ=255)
                                                 gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:24805616/1‑68 (MQ=255)
                                                 acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                   >  1:34765993/2‑68 (MQ=255)
                                                 acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                   >  1:35069360/2‑68 (MQ=255)
                                                 acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                   >  1:28461694/2‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAATTTATCCTGCTTaa                  >  1:29310161/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                  >  1:35655160/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                  >  1:34381639/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                  >  1:10358584/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                  >  1:14660618/1‑68 (MQ=255)
                                                  cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                  >  1:16845306/1‑68 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:14947450/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:13176436/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:30470906/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:4021867/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACACGTTCCGGAAATTATCCTGCTTa                   >  1:36031626/1‑67 (MQ=255)
                                                  cATATGCTGCGCCTACGTCGGCAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtaa                   >  1:3862718/1‑65 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:16356409/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                   >  1:12184708/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                   >  1:7090567/1‑66 (MQ=255)
                                                   atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCTTa                   >  1:1963288/1‑66 (MQ=255)
                                                    tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                 >  1:24252180/1‑67 (MQ=255)
                                                     atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc               <  1:10749554/68‑1 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc               >  1:23438089/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc               >  1:4834525/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc               >  1:8346460/1‑67 (MQ=255)
                                                      tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc               >  1:18205399/1‑67 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:22469968/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:11547232/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:36441355/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:19404525/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:21789976/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:26430725/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:938151/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:29460949/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:2741438/1‑68 (MQ=255)
                                                       gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac             >  1:5735920/1‑66 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:15379430/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:36717156/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:229569/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:33518845/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:20905274/1‑67 (MQ=255)
                                                        ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc             >  1:18852042/1‑67 (MQ=255)
                                                         tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa            <  1:13356642/67‑1 (MQ=255)
                                                          gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg          <  1:18499026/68‑1 (MQ=255)
                                                              cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg       >  1:25784640/1‑67 (MQ=255)
                                                               tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc     <  1:10176499/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa    <  1:37037686/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa    <  1:28922981/68‑1 (MQ=255)
                                                                aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa    <  1:1003428/68‑1 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa    >  1:16020129/1‑67 (MQ=255)
                                                                 cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCCGAAATAATCCTGCTTAACTCTCACGATGGCa    >  1:252142/1‑67 (MQ=255)
                                                                   |                                                                  
TTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACC  >  CP009273/2762220‑2762353

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: