Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,309 A→C 8.7% 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 8.7% 2006.6 / 143.9 603K96Q (AAA→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (375/175);  new base C (38/15);  total (413/190)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 649 total aligned reads displayed.

TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376
                                                                   |                                                                   
tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                     >  1:15468172/1‑68 (MQ=255)
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:17903767/1‑68 (MQ=255)
   cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                                                  >  1:23915625/1‑68 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:34179812/1‑67 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:13170270/1‑67 (MQ=255)
              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:3499921/68‑1 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                   >  1:29462491/1‑68 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                   >  1:21694550/1‑68 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                   >  1:28716282/1‑68 (MQ=255)
                    aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCCCTAAACAGGTTCCGGAAATTAt                                                  <  1:18806805/67‑1 (MQ=255)
                     tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:38144021/1‑66 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:35972012/1‑67 (MQ=255)
                       cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:20192164/1‑67 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:12766302/68‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:9204068/64‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                            >  1:25296695/4‑68 (MQ=255)
                          aGCATATGCTTCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                           >  1:12386781/1‑68 (MQ=255)
                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                           >  1:35461940/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:31670131/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:7061655/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:18116929/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                                          >  1:18432790/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAGCGGGCAGGAAGTTCCTGAAATTATCCTGCGTaa                                         >  1:27071063/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:22011970/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:26048485/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:21478570/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:24309621/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:30967784/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:11459902/1‑68 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:27844137/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:11138135/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:25614001/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:15092909/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:3949536/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:21215341/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:32672185/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCAGGAAATTATCATGCTTa                                          >  1:38893700/1‑67 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:5690446/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        <  1:1597119/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTTTCCTGCTTAAc                                        >  1:35352536/1‑67 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:26729317/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        >  1:20691830/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                                     <  1:8536071/68‑1 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:11986069/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:29202900/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:24799153/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:33889850/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:16280174/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:17000313/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:36767365/1‑67 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:10161998/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:9627290/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:33064941/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:8232208/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:4502904/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:34505055/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:3744360/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:14801249/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:20890270/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:13606898/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:23660334/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTACGGAAATTATCCTGCTTAACtctc                                    >  1:11729119/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGATCCGGAAATTATCCTGCTTAACtctc                                    >  1:32165073/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAAAAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:837377/1‑68 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:30494063/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:2995183/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:38453663/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:38664988/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:14309690/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:27376120/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:7774288/1‑67 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:2255698/68‑1 (MQ=255)
                                        cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                              >  1:15695252/1‑67 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGCTGGc                            <  1:28330133/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:5937421/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:34827809/68‑1 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:26455181/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:17703128/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:32492562/68‑1 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc                           >  1:10014583/1‑66 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           >  1:4152931/1‑67 (MQ=255)
                                              cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                       <  1:13000159/68‑1 (MQ=255)
                                              cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                       >  1:36316041/1‑68 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:19443584/68‑1 (MQ=255)
                                                ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      >  1:6609700/1‑67 (MQ=255)
                                                 gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                    >  1:4958167/1‑68 (MQ=255)
                                                   ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  >  1:3354214/1‑68 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:24557926/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:19808827/66‑1 (MQ=255)
                                                        aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg             <  1:8691163/68‑1 (MQ=255)
                                                             aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct        <  1:31996007/68‑1 (MQ=255)
                                                             aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct        <  1:29787211/68‑1 (MQ=255)
                                                               cGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctgc      <  1:37879755/68‑1 (MQ=255)
                                                                ggTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGcc     <  1:36510943/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTATCTCTCTCGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:17332460/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:23049513/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:7466173/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:10873994/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAACAGATGCTGCCggg  <  1:19241994/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACCCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:16787571/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: