Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,299 A→T 8.3% Q92H (CAA→CAT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2990AT 8.3% 2083.8 / 129.7 639Q92H (CAA→CATyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (408/178);  new base T (38/15);  total (446/193)
Fisher's exact test for biased strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 671 total aligned reads displayed.

CACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAG  >  CP009273/2762232‑2762365
                                                                   |                                                                  
 acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                   >  1:23844093/1‑68 (MQ=255)
    ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                                 >  1:23529508/1‑67 (MQ=255)
    ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                                 >  1:17553310/1‑67 (MQ=255)
        cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                            <  1:5258266/68‑1 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                           >  1:5762144/1‑68 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                           >  1:26310230/1‑68 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                         >  1:35400497/1‑68 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                         <  1:2590097/68‑1 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                         >  1:17826409/1‑68 (MQ=255)
                ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                    <  1:35191997/68‑1 (MQ=255)
                   tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                  >  1:13170270/1‑67 (MQ=255)
                       cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                             >  1:22610643/1‑68 (MQ=255)
                        gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                            <  1:8120639/68‑1 (MQ=255)
                        gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                            <  1:3836424/68‑1 (MQ=255)
                         tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                           >  1:32702849/1‑68 (MQ=255)
                            gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                        >  1:1696891/1‑68 (MQ=255)
                            gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                        >  1:4846397/1‑68 (MQ=255)
                                aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                       >  1:36466560/1‑65 (MQ=255)
                                 cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                    >  1:20192164/1‑67 (MQ=255)
                                 cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                    >  1:24162564/1‑67 (MQ=255)
                                 cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                    >  1:10123082/1‑67 (MQ=255)
                                 cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                   <  1:224120/68‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                  <  1:2450127/64‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                  <  1:38657030/64‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                  <  1:34642707/64‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                  <  1:10872406/64‑1 (MQ=255)
                                    aGCATATGCTTCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                >  1:12386781/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:34112183/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:2795211/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:28972316/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:31670131/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:7061655/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:18116929/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:1358708/1‑68 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                               >  1:18432790/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:11459902/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:21290694/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:21753957/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:37400364/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:12154951/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:25042180/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                              >  1:16008889/1‑68 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:19604682/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:28409686/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:27741480/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:6107103/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:24689239/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:16312676/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:35798874/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCAGGAAATTATCATGCTTa                               >  1:38893700/1‑67 (MQ=255)
                                       atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                               >  1:17226443/1‑66 (MQ=255)
                                        tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                            <  1:21014522/68‑1 (MQ=255)
                                         atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                           <  1:26505972/68‑1 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                          <  1:30205234/68‑1 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                           >  1:37892385/1‑67 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                           >  1:21579813/1‑67 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:38170231/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:33112417/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:32095835/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:9627290/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:29536983/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:12786450/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:19308863/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTAACtctc                         >  1:11729119/1‑68 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:30996047/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:8369284/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:36731446/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:7774288/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:3606890/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:15262827/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                         >  1:39489898/1‑67 (MQ=255)
                                              gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                      <  1:17964910/68‑1 (MQ=255)
                                                gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                     <  1:14032474/67‑1 (MQ=255)
                                                 ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                   >  1:18941359/1‑68 (MQ=255)
                                                   tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                 <  1:26815094/68‑1 (MQ=255)
                                                   tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                 <  1:11171146/68‑1 (MQ=255)
                                                    aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                <  1:32492562/68‑1 (MQ=255)
                                                     cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                >  1:10014583/1‑66 (MQ=255)
                                                     cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                >  1:4152931/1‑67 (MQ=255)
                                                        cGGGAGTTGCAAATTACCGGTAAACAGGTTCCGGAAATTATGCTGCTTAACTCTCACGATGGCACCAg            <  1:15519842/68‑1 (MQ=21)
                                                         tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:20608368/67‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:8610570/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:14646281/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:34916975/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:8525858/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:34264651/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:32231322/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           <  1:7428150/68‑1 (MQ=255)
                                                          ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt           >  1:6609700/1‑67 (MQ=255)
                                                           gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc         >  1:25419187/1‑68 (MQ=255)
                                                            aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc         <  1:9053158/67‑1 (MQ=255)
                                                            aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg        >  1:1501323/1‑68 (MQ=255)
                                                            aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg        >  1:33775378/1‑68 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt       <  1:19808827/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt       <  1:29308904/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt       <  1:27224908/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt       <  1:4467520/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt       <  1:15830548/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa      <  1:30533432/67‑1 (MQ=255)
                                                                gCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAc     <  1:13024462/67‑1 (MQ=255)
                                                                   |                                                                  
CACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAG  >  CP009273/2762232‑2762365

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: