Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,320 G→T 7.8% P99P (CCG→CCT)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3200GT 7.8% 2272.2 / 130.7 672P99P (CCG→CCT) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (376/242);  new base T (38/15);  total (414/258)
Fisher's exact test for biased strand distribution p-value = 1.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 672 total aligned reads displayed.

CGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAG  >  CP009273/2762255‑2762385
                                                                 |                                                                 
cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                                 >  1:34836713/1‑68 (MQ=255)
cgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGa                                                                 >  1:26559188/1‑68 (MQ=255)
 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                                <  1:8120639/68‑1 (MQ=255)
 gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                                <  1:16665553/68‑1 (MQ=255)
     gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                            >  1:9716805/1‑68 (MQ=255)
          cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                        >  1:10123082/1‑67 (MQ=255)
          cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                        >  1:20192164/1‑67 (MQ=255)
          cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTAGCCt                                                        >  1:27040984/1‑67 (MQ=255)
           tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                      <  1:11592911/68‑1 (MQ=255)
            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                      <  1:9204068/64‑1 (MQ=255)
            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                      <  1:35842517/64‑1 (MQ=255)
            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                                     >  1:25296695/4‑68 (MQ=255)
            aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                                     >  1:17929572/4‑68 (MQ=255)
             aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                                <  1:27295295/56‑1 (MQ=255)
             aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                    >  1:38647340/1‑68 (MQ=255)
              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                   >  1:18116929/1‑68 (MQ=255)
              gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                                                   >  1:18432790/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:3541254/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:24309621/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:23946976/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:21478570/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:37400364/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:29734959/1‑68 (MQ=255)
               cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                                  >  1:16008889/1‑68 (MQ=255)
               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                   >  1:16996040/1‑67 (MQ=255)
               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                   >  1:28434766/1‑67 (MQ=255)
               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                   >  1:21215341/1‑67 (MQ=255)
               cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                                   >  1:21947613/1‑67 (MQ=255)
                atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                                 <  1:15615714/68‑1 (MQ=255)
                  atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               <  1:2366073/68‑1 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:18696107/1‑67 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:39061877/1‑67 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:36767365/1‑67 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:17304808/1‑67 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:31321880/1‑67 (MQ=255)
                   tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                               >  1:24799153/1‑67 (MQ=255)
                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:20890270/1‑68 (MQ=255)
                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:29148810/1‑68 (MQ=255)
                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:6224096/1‑68 (MQ=255)
                    gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGATCCGGAAATTATCCTGCTTAACtctc                                             >  1:32165073/1‑68 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:32604596/1‑67 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:2995183/1‑67 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:30494063/1‑67 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:3607579/1‑67 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:38310979/1‑67 (MQ=255)
                     ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                             >  1:39489898/1‑67 (MQ=255)
                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                          <  1:2255698/68‑1 (MQ=255)
                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                          <  1:19227813/68‑1 (MQ=255)
                       gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                          <  1:16307567/68‑1 (MQ=255)
                         gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                         <  1:17632702/67‑1 (MQ=255)
                         gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                         <  1:14032474/67‑1 (MQ=255)
                         gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGa                                         <  1:12653297/67‑1 (MQ=255)
                         gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                        >  1:26295836/1‑68 (MQ=255)
                          ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                       >  1:33060804/1‑68 (MQ=255)
                            tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                     <  1:7419544/68‑1 (MQ=255)
                            tACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                     <  1:6493028/68‑1 (MQ=255)
                              cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                                    >  1:10014583/1‑66 (MQ=255)
                              cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                    >  1:12371385/1‑67 (MQ=255)
                              cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                    >  1:4152931/1‑67 (MQ=255)
                                gcGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAGCCTGCTTAACTCTCACGAGGGCACCa                                 >  1:32434524/2‑68 (MQ=21)
                                 cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                                >  1:835633/1‑68 (MQ=255)
                                  tggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:20608368/67‑1 (MQ=255)
                                  gggAGTGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:36493671/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:8525858/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:13061926/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:34082761/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:34264651/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:14646281/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:8610570/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:38817270/68‑1 (MQ=255)
                                  gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               <  1:19443584/68‑1 (MQ=255)
                                   ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                               >  1:27959606/1‑67 (MQ=255)
                                    gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                             >  1:11951661/1‑68 (MQ=255)
                                      ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           >  1:3354214/1‑68 (MQ=255)
                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           <  1:4467520/66‑1 (MQ=255)
                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           <  1:29308904/66‑1 (MQ=255)
                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           <  1:19808827/66‑1 (MQ=255)
                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           <  1:11307036/66‑1 (MQ=255)
                                        ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                           <  1:24588039/66‑1 (MQ=255)
                                              ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt                    >  1:5871974/1‑67 (MQ=255)
                                                aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct                 <  1:31996007/68‑1 (MQ=255)
                                                  cGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctgc               <  1:37879755/68‑1 (MQ=255)
                                                      aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg           <  1:25763108/68‑1 (MQ=255)
                                                      aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg           <  1:10873994/68‑1 (MQ=255)
                                                      aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg           <  1:23049513/68‑1 (MQ=255)
                                                       aaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGc          >  1:38043915/1‑68 (MQ=255)
                                                        aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa         >  1:30825691/1‑68 (MQ=255)
                                                        aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa         >  1:3515285/1‑68 (MQ=255)
                                                        aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa         <  1:13467405/68‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:7018257/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:7552714/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:3383432/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:15258181/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:5197807/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:11116356/66‑1 (MQ=255)
                                                             ttCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:15036908/66‑1 (MQ=255)
                                                             ttCAGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt      <  1:1287854/66‑1 (MQ=255)
                                                              tCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTc    <  1:3460758/67‑1 (MQ=255)
                                                              tCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTc    <  1:4812062/67‑1 (MQ=255)
                                                              tCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTc    <  1:21603999/67‑1 (MQ=255)
                                                                 |                                                                 
CGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGTTCAG  >  CP009273/2762255‑2762385

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: ‑