Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,308 T→G 8.5% G95G (GGT→GGG yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3080TG 8.5% 2155.2 / 137.2 628G95G (GGT→GGGyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (395/179);  new base G (39/15);  total (434/194)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 647 total aligned reads displayed.

GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGG  >  CP009273/2762241‑2762375
                                                                   |                                                                   
gTGACCCCCGTCGTCGTGAGCATACTTAGCATATCCTGCGCCTACGTCGGGAGGGCCAAAATACCGGt                                                                     >  1:20118287/1‑68 (MQ=255)
gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                     >  1:6274670/1‑68 (MQ=255)
gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                     >  1:19146793/1‑68 (MQ=255)
gTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                     <  1:2779193/68‑1 (MQ=255)
 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                    >  1:15468172/1‑68 (MQ=255)
  gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                   >  1:13820386/1‑68 (MQ=255)
  gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                    <  1:21185745/67‑1 (MQ=255)
   aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                   >  1:23634805/1‑67 (MQ=255)
    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                 >  1:23915625/1‑68 (MQ=255)
    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                 >  1:19756477/1‑68 (MQ=255)
      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                               >  1:31906057/1‑68 (MQ=255)
      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                               >  1:26723228/1‑68 (MQ=255)
       ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                              <  1:22445930/68‑1 (MQ=255)
              cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                       >  1:32967020/1‑68 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                      <  1:17237376/68‑1 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                      <  1:36729434/68‑1 (MQ=255)
                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                  >  1:35192212/1‑68 (MQ=255)
                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                  >  1:22868885/1‑68 (MQ=255)
                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                  >  1:13188643/1‑68 (MQ=255)
                     aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                 <  1:37953113/67‑1 (MQ=255)
                     aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCCCTAAACAGGTTCCGGAAATTAt                                                 <  1:18806805/67‑1 (MQ=255)
                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                 >  1:24732509/1‑66 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                 >  1:32634548/1‑65 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                 >  1:8155601/1‑65 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                 >  1:36466560/1‑65 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                              <  1:5308271/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                            <  1:11504983/68‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                            <  1:12469876/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                            <  1:34642707/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                            <  1:35842517/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                            <  1:10872406/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                            <  1:7676932/64‑1 (MQ=255)
                           aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                      <  1:27295295/56‑1 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:19394956/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:18352508/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:1279666/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:34112183/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:18116929/1‑68 (MQ=255)
                            acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                         >  1:38528130/2‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:26048485/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:4063398/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:3541254/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:37777499/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:21290694/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:21753957/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:24309621/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                        >  1:12154951/1‑68 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:16996040/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:30665976/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:28434766/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:11758459/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:15092909/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:26344559/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:17902881/1‑67 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:5816887/1‑66 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                         >  1:2158446/1‑66 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                     >  1:27043654/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                     >  1:25271490/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                     >  1:31321880/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                     >  1:33393002/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                     >  1:32420721/1‑67 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:34965355/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:38355018/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:37076615/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:4389475/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:9627290/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:32142443/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:10161998/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:20844515/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:14801249/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:28109961/1‑68 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                                  >  1:15902398/1‑36 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTATCCTGCTGAACtctc                                   >  1:16227905/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:3607579/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:24542710/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:14309690/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:6819477/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:8369284/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                   >  1:33828871/1‑67 (MQ=255)
                                         cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                             >  1:15695252/1‑67 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                           <  1:11171146/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                          <  1:20401150/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                          <  1:17703128/68‑1 (MQ=255)
                                            cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                          >  1:10014583/1‑66 (MQ=255)
                                             gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc                        >  1:7245991/1‑68 (MQ=255)
                                             gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc                        >  1:29936908/1‑68 (MQ=255)
                                               cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                      <  1:2559142/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:38817270/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:13061926/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:34264651/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:8525858/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:8610570/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     <  1:30396738/68‑1 (MQ=255)
                                                 ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     >  1:39451025/1‑67 (MQ=255)
                                                 ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                     >  1:29650323/1‑67 (MQ=255)
                                                   aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                   <  1:9053158/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                 <  1:29308904/66‑1 (MQ=255)
                                                             tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgc        <  1:4782553/68‑1 (MQ=255)
                                                                 ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc    >  1:2915426/1‑68 (MQ=255)
                                                                   tAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCgg  >  1:16513033/1‑68 (MQ=255)
                                                                   |                                                                   
GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGG  >  CP009273/2762241‑2762375

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: