Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,284 A→T 10.3% 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,2840AT10.3% 1664.2 / 144.5 561L87L (CTA→CTTyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (348/155);  new base T (48/10);  total (396/165)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 576 total aligned reads displayed.

TCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCA  >  CP009273/2762217‑2762351
                                                                   |                                                                   
  tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                                   <  1:29197408/68‑1 (MQ=255)
  tttGCCTGTCAGACACGGGTACGTGACCCCCGTCATCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                                   <  1:17019363/68‑1 (MQ=255)
      ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               >  1:7096778/1‑68 (MQ=255)
      ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                               >  1:27797008/1‑68 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:13874118/68‑1 (MQ=255)
           cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                          <  1:33864798/68‑1 (MQ=255)
               cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                        <  1:36525876/66‑1 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:32109393/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:23352254/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:8967265/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:22380958/1‑68 (MQ=255)
                acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                     >  1:20699785/1‑68 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:965517/1‑67 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:2959935/1‑67 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:32760632/1‑67 (MQ=255)
                   ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                                   >  1:12700037/1‑67 (MQ=255)
                    gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                                  >  1:28919407/1‑67 (MQ=255)
                       cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                              <  1:7762291/68‑1 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             <  1:5225803/68‑1 (MQ=255)
                        gTGACCCCCGTCGTCGTGAGCATACTAAGCAAATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                             >  1:28419807/1‑68 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                            >  1:27430538/1‑68 (MQ=255)
                         tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                            >  1:13190941/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTGGGGAGGGGCAAATAACCGGTaa                                           >  1:22636264/1‑68 (MQ=255)
                          gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                           >  1:12997664/1‑68 (MQ=255)
                              cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                       >  1:33431269/1‑68 (MQ=255)
                                cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                     <  1:35951686/68‑1 (MQ=255)
                                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              <  1:14559307/68‑1 (MQ=255)
                                       gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                              >  1:25777919/1‑68 (MQ=255)
                                             aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         <  1:19811091/67‑1 (MQ=255)
                                             aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         <  1:20162177/67‑1 (MQ=255)
                                              tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:16814090/1‑66 (MQ=255)
                                              tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:5675739/1‑66 (MQ=255)
                                              tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:16042058/1‑66 (MQ=255)
                                              tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                         >  1:36380048/1‑66 (MQ=255)
                                               aCTAAGCATATGCTGCGCCTACGTCGGGTGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:24138420/1‑68 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                      >  1:33565807/1‑66 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:26352666/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:26802737/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:27071715/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:35581714/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:24665125/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:8587115/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:12321554/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:10996108/68‑1 (MQ=255)
                                                 tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:19005533/68‑1 (MQ=255)
                                                  aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                    <  1:23548057/64‑1 (MQ=255)
                                                   aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                    <  1:35174342/66‑1 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:25038051/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:37088402/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:32363425/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:244936/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:10358584/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:24846290/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                >  1:10807458/1‑68 (MQ=255)
                                                     cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTCaa                >  1:33178475/1‑68 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:6974892/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:34133907/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:16437649/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:23834694/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:31718595/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:2181347/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:13399908/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:23086310/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:19249605/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:5965967/1‑67 (MQ=255)
                                                     cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:20887974/1‑67 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:4651883/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:3913387/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                 >  1:21701616/1‑66 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:15003050/68‑1 (MQ=255)
                                                      atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc               <  1:18093018/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct              <  1:30890128/68‑1 (MQ=255)
                                                       tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct              <  1:8404659/68‑1 (MQ=255)
                                                       aatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTCCTTAAc               >  1:29423744/2‑67 (MQ=255)
                                                        atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             <  1:18535626/68‑1 (MQ=255)
                                                         tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc             >  1:21150940/1‑67 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:11547232/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:17371699/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:1467280/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:15458161/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:34571967/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:15816364/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:18293567/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:36060846/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:25459747/1‑68 (MQ=255)
                                                          gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:938151/1‑68 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:36717156/1‑67 (MQ=255)
                                                           ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc           >  1:6576684/1‑67 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:34880528/68‑1 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:6235986/68‑1 (MQ=255)
                                                             gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg        <  1:22203588/68‑1 (MQ=255)
                                                               gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt      >  1:4082943/1‑68 (MQ=255)
                                                                ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     >  1:31348907/1‑68 (MQ=255)
                                                                ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTACCGGAAATTATCCTGCTTAACTCTCACGATg     >  1:1947982/1‑68 (MQ=255)
                                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     <  1:30227307/67‑1 (MQ=255)
                                                                 cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg     <  1:20496099/67‑1 (MQ=255)
                                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:10176499/68‑1 (MQ=255)
                                                                  tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc   <  1:4936709/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:26003971/68‑1 (MQ=255)
                                                                   aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa  <  1:1003428/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: