Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,314 G→A 8.6% Q97Q (CAG→CAA) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3140GA 8.6% 1882.3 / 130.6 615Q97Q (CAG→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (371/191);  new base A (38/15);  total (409/206)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 638 total aligned reads displayed.

CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381
                                                                   |                                                                   
cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                     >  1:27461473/1‑68 (MQ=255)
    tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                                  >  1:27939580/1‑67 (MQ=255)
    tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                                  >  1:13170270/1‑67 (MQ=255)
       tcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                                    <  1:32756381/62‑1 (MQ=255)
        cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                             >  1:32141972/1‑68 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:3499921/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:20368943/68‑1 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                         <  1:17083761/67‑1 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:9716805/1‑68 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                        >  1:21694550/1‑68 (MQ=255)
             gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCCGAAAtt                                                         <  1:21490720/67‑1 (MQ=255)
                tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:9355795/1‑66 (MQ=255)
                 aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:11611314/1‑65 (MQ=255)
                  cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:35972012/1‑67 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:2926060/1‑66 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:31866954/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:24884327/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:18287022/68‑1 (MQ=255)
                    aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                  <  1:38657030/64‑1 (MQ=255)
                    aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                                 >  1:25296695/4‑68 (MQ=255)
                     aGCATATGCTTCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                >  1:12386781/1‑68 (MQ=255)
                     aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                >  1:35461940/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:1358708/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:28972316/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:7061655/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:16904944/1‑68 (MQ=255)
                       cttatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:15053780/3‑67 (MQ=255)
                       cATGTGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:28403635/1‑67 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAGCGGGCAGGAAGTTCCTGAAATTATCCTGCGTaa                                              >  1:27071063/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:11459902/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:37400364/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:26048485/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:16008889/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:22011970/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:23952795/1‑68 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:30665976/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:11138135/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:26344559/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:16312676/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:17850023/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCAGGAAATTATCATGCTTa                                               >  1:38893700/1‑67 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:5690446/1‑66 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                            <  1:21014522/68‑1 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                             >  1:20691830/1‑67 (MQ=255)
                          atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           <  1:2366073/68‑1 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                                          <  1:8536071/68‑1 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:11986069/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:36767365/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:15394670/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:19311738/1‑67 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:33112417/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:34174338/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:10161998/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:4389475/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:38355018/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:34505055/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:29536983/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:20033362/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:21254811/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:1746369/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:12786450/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:837377/1‑68 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:18983887/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:12959965/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:7774288/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:24357782/1‑67 (MQ=255)
                               gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                      <  1:17964910/68‑1 (MQ=255)
                               gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                      <  1:32377299/68‑1 (MQ=255)
                               gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                      <  1:2255698/68‑1 (MQ=255)
                                   cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:15695252/1‑67 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:5937421/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:34827809/68‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:26815094/68‑1 (MQ=255)
                                      cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                                >  1:10014583/1‑66 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:2559142/68‑1 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            >  1:36316041/1‑68 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            >  1:36498219/1‑68 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:19627252/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:7428150/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:14646281/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:30396738/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:35255055/68‑1 (MQ=255)
                                           ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           >  1:39451025/1‑67 (MQ=255)
                                           ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           >  1:6609700/1‑67 (MQ=255)
                                            gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                         >  1:4958167/1‑68 (MQ=255)
                                             aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                        >  1:33775378/1‑68 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:31314422/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:24588039/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:8228079/66‑1 (MQ=255)
                                                   aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg                  <  1:8691163/68‑1 (MQ=255)
                                                     aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                <  1:4147321/68‑1 (MQ=255)
                                                      ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                >  1:18524223/1‑67 (MQ=255)
                                                         ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg            <  1:29927868/68‑1 (MQ=255)
                                                          cGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctgc           <  1:37879755/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:10873994/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:23049513/68‑1 (MQ=255)
                                                               aaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGc      >  1:38043915/1‑68 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     >  1:33911287/1‑68 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     >  1:32591472/1‑68 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     >  1:14188203/1‑68 (MQ=255)
                                                                   |                                                                   
CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381

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: ‑