Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,308 T→G 10.4% G95G (GGT→GGG)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,3080TG10.4% 1965.8 / 158.9 600G95G (GGT→GGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (314/223);  new base G (55/8);  total (369/231)
Fisher's exact test for biased strand distribution p-value = 3.75e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 619 total aligned reads displayed.

GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCG  >  CP009273/2762241‑2762374
                                                                   |                                                                  
gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                    >  1:17598364/1‑68 (MQ=255)
gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                    <  1:20619257/68‑1 (MQ=255)
gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                                    >  1:8894107/1‑68 (MQ=255)
  gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                  <  1:27934808/68‑1 (MQ=255)
  gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                  >  1:20262975/1‑68 (MQ=255)
   aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                  >  1:6522539/1‑67 (MQ=255)
    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                                                 <  1:3242139/67‑1 (MQ=255)
    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                >  1:19743859/1‑68 (MQ=255)
      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                              >  1:6076309/1‑68 (MQ=255)
      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                              >  1:21969064/1‑68 (MQ=255)
      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                              >  1:12250712/1‑68 (MQ=255)
          tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                           >  1:21605357/1‑67 (MQ=255)
               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                     <  1:25997970/68‑1 (MQ=255)
                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTAt                                                >  1:33333247/1‑66 (MQ=255)
                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGAAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                >  1:26302821/1‑65 (MQ=255)
                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                             >  1:534547/1‑67 (MQ=255)
                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                             >  1:31355089/1‑67 (MQ=255)
                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                            <  1:937904/68‑1 (MQ=255)
                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                                            <  1:34441504/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                           <  1:22161951/68‑1 (MQ=255)
                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                           <  1:834361/68‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                           <  1:20907736/64‑1 (MQ=255)
                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                           <  1:16922034/64‑1 (MQ=255)
                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                          >  1:5635965/1‑68 (MQ=255)
                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                          >  1:31541674/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:15258948/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:34808626/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:10034764/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:27074642/1‑68 (MQ=255)
                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:12696702/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:30905236/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:32016544/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:29343515/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:13645974/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:19468235/1‑68 (MQ=255)
                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                       >  1:17193476/1‑68 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTaa                                       >  1:33949800/1‑68 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:28941115/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:2856538/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:15553206/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:24078713/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:16690053/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:17849792/1‑67 (MQ=255)
                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:18470683/1‑67 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                        >  1:30268990/1‑66 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                      <  1:16374527/68‑1 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                      <  1:23284781/68‑1 (MQ=255)
                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                      <  1:14584906/68‑1 (MQ=255)
                               tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                      >  1:24728653/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                    >  1:21373711/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                    >  1:14940136/1‑67 (MQ=255)
                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                    >  1:31237878/1‑67 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:7251605/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:29620160/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:11793163/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:26733988/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:26649578/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:32609843/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:22695066/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:18059418/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:15996758/1‑68 (MQ=255)
                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCTTAACtctc                                  >  1:13787432/1‑68 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:9800992/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:14234876/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:10718794/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                  >  1:11295086/1‑67 (MQ=255)
                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATACTGCATAACtctc                                  >  1:13044381/1‑67 (MQ=255)
                                    tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                                 <  1:10902952/67‑1 (MQ=255)
                                    tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                >  1:2991751/1‑68 (MQ=255)
                                    tgcGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                >  1:25637414/1‑68 (MQ=255)
                                       gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                             >  1:11542573/1‑68 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                          <  1:28279957/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                          <  1:9160478/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                          <  1:7788422/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                          <  1:6909918/68‑1 (MQ=255)
                                          tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                          <  1:4451475/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTTAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                         <  1:26475924/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                         <  1:23057966/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                         <  1:10311697/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                         <  1:30599487/68‑1 (MQ=255)
                                           aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                         <  1:27763072/68‑1 (MQ=255)
                                              tCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCa                      >  1:7970673/1‑68 (MQ=255)
                                               cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                     <  1:2870581/68‑1 (MQ=255)
                                                gggAGGGGCAATTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:18671489/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:33528953/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:11094361/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:20128413/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:402710/68‑1 (MQ=255)
                                                gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                    <  1:3573912/68‑1 (MQ=255)
                                                  gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                  >  1:1911712/1‑68 (MQ=255)
                                                   aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                  <  1:4972355/67‑1 (MQ=255)
                                                   aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                 >  1:26847937/1‑68 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                <  1:17251444/66‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa               <  1:21023737/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa               <  1:6350563/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa               <  1:25451937/67‑1 (MQ=255)
                                                      ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa               <  1:8757092/67‑1 (MQ=255)
                                                             tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt         >  1:25142604/1‑66 (MQ=255)
                                                              aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct      <  1:23690361/68‑1 (MQ=255)
                                                              aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct      <  1:24903795/68‑1 (MQ=255)
                                                                   |                                                                  
GTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCG  >  CP009273/2762241‑2762374

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