Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,311 A→G 10.5% K96K (AAA→AAG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3110AG10.5% 1713.5 / 130.7 550K96K (AAA→AAG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (296/195);  new base G (48/10);  total (344/206)
Fisher's exact test for biased strand distribution p-value = 8.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 575 total aligned reads displayed.

CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378
                                                                  |                                                                   
cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                    >  1:17279396/1‑68 (MQ=255)
  cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                  >  1:7621365/1‑68 (MQ=255)
          cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                          >  1:25731211/1‑68 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:20849463/68‑1 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:14589935/68‑1 (MQ=255)
           gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                         <  1:3814604/68‑1 (MQ=255)
            tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCATATTACCGGTAAACAGGTTCCGGaa                                                         <  1:13601350/67‑1 (MQ=255)
               gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                     >  1:29090370/1‑68 (MQ=255)
                 aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                    <  1:19811091/67‑1 (MQ=255)
                  tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                    >  1:4417212/1‑66 (MQ=255)
                    cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:20552107/1‑67 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:33565807/1‑66 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:13334363/1‑66 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:32462362/1‑66 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                 >  1:1580181/1‑66 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:32255920/68‑1 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:8587115/68‑1 (MQ=255)
                     tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                               <  1:18519850/68‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                               <  1:503113/64‑1 (MQ=255)
                      aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                              >  1:33436514/4‑68 (MQ=255)
                      aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                              >  1:25439470/1‑68 (MQ=255)
                       aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                             >  1:10268429/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:32782061/1‑68 (MQ=255)
                        gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:28770668/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAATTTATCCTGCTTaa                                           >  1:29310161/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:13153286/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:8358526/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:24846290/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:22511800/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:10807458/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:28529388/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:16845306/1‑68 (MQ=255)
                         cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                           >  1:27593488/1‑68 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:26310310/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:25865058/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:6699370/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:2181347/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:19249605/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:22258793/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:23403951/1‑67 (MQ=255)
                         cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:14947450/1‑67 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:35955362/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:23718973/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:21701616/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:16356409/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                            >  1:1908363/1‑66 (MQ=255)
                          atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                          <  1:18093018/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                         <  1:23057566/68‑1 (MQ=255)
                           tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                         <  1:30890128/68‑1 (MQ=255)
                            atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        <  1:11056360/68‑1 (MQ=255)
                             tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                        >  1:24356355/1‑67 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:29460949/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:18293567/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:36060846/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:16569908/1‑68 (MQ=255)
                              gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac                                      >  1:5735920/1‑66 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:15379430/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:15945357/1‑67 (MQ=255)
                               ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                      >  1:29572030/1‑67 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:24008604/1‑68 (MQ=255)
                                tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                    >  1:31785292/1‑68 (MQ=255)
                                 gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                   <  1:14465818/68‑1 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                 >  1:26866613/1‑68 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                 >  1:34966877/1‑68 (MQ=255)
                                   gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                 >  1:33196554/1‑68 (MQ=255)
                                    ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                >  1:6149790/1‑68 (MQ=255)
                                    ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTACCGGAAATTATCCTGCTTAACTCTCACGATg                                >  1:1947982/1‑68 (MQ=255)
                                     cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATgg                               <  1:33928722/68‑1 (MQ=255)
                                     cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                <  1:30227307/67‑1 (MQ=255)
                                     cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                <  1:11468439/67‑1 (MQ=255)
                                       aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             <  1:1003428/68‑1 (MQ=255)
                                        cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                             >  1:26531393/1‑67 (MQ=255)
                                           cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                         >  1:35156928/1‑68 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:34300834/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:20141685/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:27155470/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:12440422/68‑1 (MQ=255)
                                            gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACACTCACGATGGCACCAGt                        <  1:9757184/68‑1 (MQ=255)
                                            gggACGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        <  1:34691368/68‑1 (MQ=255)
                                             ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                        >  1:3953232/1‑67 (MQ=255)
                                               aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                     >  1:4719760/1‑68 (MQ=255)
                                                ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    >  1:9313480/1‑68 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:36262659/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:3146571/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:12693931/66‑1 (MQ=255)
                                                  ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                    <  1:35580214/66‑1 (MQ=255)
                                                      aaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGa              >  1:27064246/1‑68 (MQ=255)
                                                       aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt             >  1:7060932/1‑68 (MQ=255)
                                                         tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt             >  1:5438967/1‑66 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct          <  1:24604743/68‑1 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct          <  1:36488501/68‑1 (MQ=255)
                                                          aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct          <  1:2806791/68‑1 (MQ=255)
                                                               tAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCgg     >  1:18878215/1‑68 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:22131470/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:25123663/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:11878761/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:10482520/68‑1 (MQ=255)
                                                                aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg    <  1:30493125/68‑1 (MQ=255)
                                                                  aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa  >  1:36901439/1‑68 (MQ=255)
                                                                  aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa  <  1:8876275/68‑1 (MQ=255)
                                                                  |                                                                   
CCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCA  >  CP009273/2762245‑2762378

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: