Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,278 G→A 9.2% L85L (CTG→CTA) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2780GA 9.2% 1222.0 / 90.3 402L85L (CTG→CTA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (207/158);  new base A (24/13);  total (231/171)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

Only 100 of 556 total aligned reads displayed.

GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345
                                                                  |                                                                   
gCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    >  1:23027558/1‑68 (MQ=255)
 ccATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    >  1:3125641/1‑67 (MQ=255)
  cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                    <  1:29218558/66‑1 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                                >  1:13618681/1‑67 (MQ=255)
     tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                               <  1:27228461/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:31630845/68‑1 (MQ=255)
      cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                              <  1:10032809/68‑1 (MQ=255)
       tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                             >  1:11559863/1‑68 (MQ=255)
        ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                            <  1:22664346/68‑1 (MQ=255)
          gCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCg                                                          >  1:27862253/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:20403532/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:212229/1‑68 (MQ=255)
           ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                         >  1:21870817/1‑68 (MQ=255)
                cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                    <  1:297422/68‑1 (MQ=255)
                    cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                >  1:16539442/1‑68 (MQ=255)
                     acGGGTACGTGCCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               >  1:19331344/1‑68 (MQ=255)
                     acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                               <  1:8936032/68‑1 (MQ=255)
                           aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                          <  1:13117900/67‑1 (MQ=255)
                            cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                        >  1:26441977/1‑68 (MQ=255)
                            cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                        <  1:16205691/68‑1 (MQ=255)
                             gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                       >  1:3708005/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:31843501/1‑68 (MQ=255)
                              tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                      >  1:32183421/1‑68 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:25489408/1‑68 (MQ=255)
                                 cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                   >  1:20301986/1‑68 (MQ=255)
                                   cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                 >  1:26061048/1‑68 (MQ=255)
                                    ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                <  1:3565559/68‑1 (MQ=255)
                                      gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                              <  1:12028300/68‑1 (MQ=255)
                                       tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                              >  1:10364194/1‑67 (MQ=255)
                                       tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGtt                              >  1:16465419/1‑67 (MQ=255)
                                           cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                         >  1:19081080/1‑68 (MQ=255)
                                            gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                        <  1:12451692/68‑1 (MQ=255)
                                                   tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                   >  1:15490234/1‑66 (MQ=255)
                                                    aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                <  1:3253376/68‑1 (MQ=255)
                                                      tAAGCTTATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:22183774/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:11703265/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:20001750/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:21641196/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:33697976/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:2949069/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:15728144/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:5937424/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:14182409/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:24389089/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:14790800/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:25154990/68‑1 (MQ=255)
                                                      tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc              <  1:33318406/68‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:3027256/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:13231038/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:17195611/64‑1 (MQ=255)
                                                       aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc              <  1:18228964/64‑1 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:34094856/1‑68 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:8150180/1‑68 (MQ=255)
                                                       aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt             >  1:23774747/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:12625600/1‑68 (MQ=255)
                                                         gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:9493130/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:11259431/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          >  1:30754710/1‑68 (MQ=255)
                                                          cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa          <  1:7816833/68‑1 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTATCCTGCTTa           >  1:29930086/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:4403429/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:17680959/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:10982897/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:5073050/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:576194/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:34372039/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:10696727/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:15899599/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:28596817/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:18866903/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:7420498/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:27495277/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:6982336/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGCTTCCGGAAATTATCCTGCTTa           >  1:4745800/1‑67 (MQ=255)
                                                          cATATGCTGCGCCTAAGTCGGGAGGGGCAAATTACCGGGAAACAGGTTCCGGAAATTATCCTGCttt           >  1:28023321/1‑66 (MQ=25)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa           >  1:1959787/1‑66 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:15177738/68‑1 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:24761432/68‑1 (MQ=255)
                                                           atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc         <  1:20962080/68‑1 (MQ=255)
                                                             atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       <  1:32811250/68‑1 (MQ=255)
                                                             atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       <  1:3863056/68‑1 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:7097306/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:9818991/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:6627186/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:4253471/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:33420592/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:29015933/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:28357346/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:26156823/1‑67 (MQ=255)
                                                              tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc       >  1:12852906/1‑67 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATc                  >  1:30875966/1‑55 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:4191224/1‑68 (MQ=255)
                                                               gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:18406827/1‑68 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:5369383/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:6212445/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:26739706/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:24049103/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc     >  1:12204461/1‑67 (MQ=255)
                                                                ctgcGCCTACGTCGGGAGGGGCAAATTACCCGTAAACAGGTTCCGGAAAATATCCTGCTTCACtctc     >  1:26921056/1‑67 (MQ=255)
                                                                 tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa    <  1:14514245/67‑1 (MQ=255)
                                                                  |                                                                   
GCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACG  >  CP009273/2762212‑2762345

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: