Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,309 A→C 6.4% K96Q (AAA→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,3090AC 6.4% 2023.7 / 90.7 574K96Q (AAA→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (312/225);  new base C (24/13);  total (336/238)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 626 total aligned reads displayed.

TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376
                                                                   |                                                                   
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:5635480/1‑68 (MQ=255)
     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAActg                                                                <  1:6066273/68‑3 (MQ=255)
     cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                >  1:14222165/1‑68 (MQ=255)
        gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             <  1:12028300/68‑1 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGtt                                                             >  1:16465419/1‑67 (MQ=255)
             cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                        >  1:26104866/1‑68 (MQ=255)
              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:23328282/68‑1 (MQ=255)
              gggAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:30695846/66‑1 (MQ=255)
                  gCATACTAAGCATATTCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                    <  1:19576734/67‑1 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                    <  1:9605558/67‑1 (MQ=255)
                    aTACTAAGCATATGCTGCGACTACGTCGGGAGGGGCAAATTACCTGTAAAAAGGTTCCGGAAATTAt                                                  <  1:7031410/67‑1 (MQ=21)
                     tACTAAGCATATGCTGCGCCTACGTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:4359927/1‑66 (MQ=255)
                      aCTACGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               <  1:3984851/68‑1 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               <  1:27005721/68‑1 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:10054902/1‑67 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:31701875/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:14790800/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:27467311/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:15728144/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:11697434/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:22578419/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:18388551/68‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:20373379/64‑1 (MQ=255)
                         aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                            >  1:23774747/1‑68 (MQ=255)
                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                           >  1:3718783/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:12625600/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:6754655/1‑68 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          >  1:9834976/2‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:1513669/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:9152410/1‑68 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTATCCTGCTTa                                          >  1:29930086/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:22053599/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:25765953/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:34372039/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:26809866/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:27734103/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:576194/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:11214640/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:7605476/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:7859267/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:13766170/1‑67 (MQ=255)
                            cATATGCTGCGCCTAAGTCGGGAGGGGCAAATTACCGGGAAACAGGTTCCGGAAATTATCCTGCttt                                          >  1:28023321/1‑66 (MQ=25)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGTTAAACAGGTTCCGGAAATTATCCAGCTTa                                          >  1:10852022/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:24655138/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:20962080/68‑1 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:15558156/68‑1 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:19866412/68‑1 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:30271242/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:3553306/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        >  1:25191035/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       >  1:4234062/1‑50 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:2014231/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:33064135/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:19209931/1‑67 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:18708903/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:32521499/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:12801894/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:16540241/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:29242106/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:29611470/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:17761472/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGATTAACtctc                                    >  1:32772016/1‑68 (MQ=255)
                                 gctgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTGCCGGAAATTATCCTGCTTAACtctc                                    >  1:33377422/1‑68 (MQ=18)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:12308943/1‑67 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                   <  1:28461007/67‑1 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                   <  1:14514245/67‑1 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                   <  1:13015530/67‑1 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:26476535/68‑1 (MQ=255)
                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      <  1:19078129/61‑1 (MQ=255)
                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                <  1:33511931/67‑1 (MQ=255)
                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                               >  1:30980546/1‑68 (MQ=255)
                                        cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                              <  1:5200599/67‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:21574119/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:33744171/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:15978563/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:4988719/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:26370762/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:13320837/68‑1 (MQ=255)
                                            gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc                         >  1:6519765/1‑68 (MQ=255)
                                              cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                       <  1:17127602/68‑1 (MQ=255)
                                               tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:23963657/67‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:2496031/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:11511646/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:4550388/68‑1 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                   >  1:34063208/1‑68 (MQ=255)
                                                   ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  >  1:84283/1‑68 (MQ=255)
                                                   ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  >  1:32086254/1‑68 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:16157604/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:22921781/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:2735594/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:21272714/66‑1 (MQ=255)
                                                        aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg             <  1:30850869/68‑1 (MQ=255)
                                                          aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt           <  1:17509330/68‑1 (MQ=255)
                                                                ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc     <  1:10614884/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTTCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:28593936/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:24321795/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:21797624/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:7317445/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:29064974/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:10423735/68‑1 (MQ=255)
                                                                   |                                                                   
TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376

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: