Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,284 A→T 8.4% L87L (CTA→CTT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2840AT 8.4% 1993.9 / 128.1 633L87L (CTA→CTTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (423/156);  new base T (38/15);  total (462/171)
Fisher's exact test for biased strand distribution p-value = 8.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

Only 100 of 647 total aligned reads displayed.

TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351
                                                                   |                                                                   
tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                                     <  1:24769964/68‑1 (MQ=255)
 cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                    <  1:25480666/68‑1 (MQ=255)
 cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                                    <  1:34898260/68‑1 (MQ=255)
   ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                  <  1:22921454/68‑1 (MQ=255)
   ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                                  <  1:17962240/68‑1 (MQ=255)
     gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                                <  1:14111529/68‑1 (MQ=255)
      ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               >  1:27302998/1‑68 (MQ=255)
       cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               <  1:12321643/67‑1 (MQ=255)
           cAGACCCGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                           <  1:24287793/67‑1 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:5600058/68‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                        <  1:18556785/66‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                      >  1:38174597/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGTAGGGGCaaa                                                     >  1:3903720/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     <  1:10039887/68‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:4738890/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:15762944/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:8052469/1‑68 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:5903515/1‑67 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:400852/1‑67 (MQ=255)
                      aCGTGACCCCCGTCGGCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                <  1:22179674/67‑1 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:6274670/1‑68 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             <  1:30255016/68‑1 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCTGTa                                            >  1:5071032/1‑68 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                            >  1:29532381/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:13820386/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:17826409/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:35400497/1‑68 (MQ=255)
                           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:4243609/1‑67 (MQ=255)
                           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:15987214/1‑67 (MQ=255)
                           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:23634805/1‑67 (MQ=255)
                           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:26642169/1‑67 (MQ=255)
                           aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:38454370/1‑67 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                         >  1:23915625/1‑68 (MQ=255)
                            cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                         >  1:19756477/1‑68 (MQ=255)
                                  tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                    >  1:13170270/1‑67 (MQ=255)
                                   cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                                  >  1:15492339/1‑68 (MQ=255)
                                      cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                               >  1:32141972/1‑68 (MQ=255)
                                      cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                               >  1:32967020/1‑68 (MQ=255)
                                           gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTa                          >  1:34242273/1‑68 (MQ=255)
                                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:36466560/1‑65 (MQ=255)
                                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:32634548/1‑65 (MQ=255)
                                               aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:11611314/1‑65 (MQ=255)
                                                cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                     <  1:17212806/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:2926060/1‑66 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:31866954/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:18309472/68‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:7676932/64‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:37939870/64‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:9204068/64‑1 (MQ=255)
                                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:16904944/1‑68 (MQ=255)
                                                    gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:2795211/1‑68 (MQ=255)
                                                    acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                 >  1:35035228/2‑68 (MQ=255)
                                                    acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                 >  1:26324433/2‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:27032757/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:11459902/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:28912256/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:21753957/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:9624550/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:37052456/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:8481953/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:35798874/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:19393385/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:28409686/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:27741480/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:18100352/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:21215341/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:11138135/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:14524379/1‑67 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAATTTATCCTGCTTa                 >  1:10812979/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:8279610/1‑66 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:25271490/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:31321880/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:21579813/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:11986069/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:33856206/1‑67 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:33393002/1‑67 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:20033362/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:14973123/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:15212673/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:8744137/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:10161998/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:3744360/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:20890270/1‑68 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:38664988/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:34515959/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:24542710/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:6819477/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:12959965/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:3606890/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:16256160/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:17467332/1‑67 (MQ=255)
                                                            tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc         >  1:30775219/1‑68 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:2255698/68‑1 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:19101686/68‑1 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:32377299/68‑1 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa       <  1:12653297/67‑1 (MQ=255)
                                                                ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     >  1:2603640/1‑68 (MQ=255)
                                                                ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     >  1:29830575/1‑68 (MQ=255)
                                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:7419544/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:20401150/68‑1 (MQ=255)
                                                                   |                                                                   
TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351

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: