Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 284,986 A→C 11.3% intergenic (‑113/+26) yagI ← / ← argF CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273284,9860AC11.3% 509.7 / 4.7 167intergenic (‑113/+26)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base A (135/6);  new base C (3/15);  total (138/29)
Fisher's exact test for biased strand distribution p-value = 9.67e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTA  >  CP009273/284923‑285051
                                                                |                                                                  
ggATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                               >  1:9184095/1‑68 (MQ=255)
 gATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                               >  1:26844825/1‑67 (MQ=255)
       aaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAc                                                         >  1:17729473/1‑67 (MQ=255)
       aaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAAc                                                         >  1:26090692/1‑67 (MQ=255)
         aTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAAAAAcgc                                                       >  1:5653256/1‑67 (MQ=255)
                 aaTGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                             <  1:3429026/68‑1 (MQ=255)
                  aTGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAc                                              >  1:12374173/1‑67 (MQ=255)
                    ggAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCa                                           >  1:9802201/1‑68 (MQ=255)
                    ggAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCa                                           >  1:4819560/1‑68 (MQ=255)
                      aGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCACCCACCCGCAGCGGACCGaa                                          >  1:22207909/1‑67 (MQ=18)
                       gCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATc                                        >  1:26361478/1‑68 (MQ=255)
                       gCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATc                                        >  1:21314982/1‑68 (MQ=255)
                        ccGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATc                                        >  1:1367375/1‑67 (MQ=255)
                        ccGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATc                                        >  1:21093139/1‑67 (MQ=255)
                         cGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATc                                        >  1:15101558/1‑66 (MQ=255)
                           gtgtTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGTACCAACGCGCAGCGGACCCAATCACt                                     >  1:18039040/1‑67 (MQ=255)
                           gtgtTCACCCGCCCCGTGGGGTCGTGGGCGGGGGTGCGGCCACCAACGCGCAGCGGACCCAATCACt                                     <  1:18693097/67‑1 (MQ=255)
                           gtgtCCCCCCTCCCCGGGGGGGAGGGGGCGGGGGGGCGGCCACCAACGCGCAGCGGACCCAATCACt                                     <  1:28782200/67‑1 (MQ=255)
                            tgtTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTc                                    >  1:7528002/1‑67 (MQ=255)
                              ttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcccc                                 >  1:27503070/1‑68 (MQ=255)
                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCCCTCCCCAAGc                             >  1:28047584/1‑68 (MQ=255)
                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGc                             >  1:2267049/1‑68 (MQ=255)
                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGc                             >  1:15807458/1‑68 (MQ=255)
                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGc                             >  1:1385127/1‑68 (MQ=255)
                                  ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGc                             >  1:11335219/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCACTCACTCCCCAAGCGt                           >  1:24864452/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:14009834/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:15517803/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:9666471/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:10382083/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:17706092/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:6435185/1‑68 (MQ=255)
                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGt                           >  1:1280725/1‑68 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCCGCGGCCCCAATCACTCCCCACGCGGTg                         >  1:6915096/1‑67 (MQ=21)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCCGCGGACCCAATCACTCCCCAAGCGTTg                         >  1:18097050/1‑67 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                         >  1:8049289/1‑67 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                         >  1:10120097/1‑67 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                         >  1:14077662/1‑67 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                         >  1:2082720/1‑67 (MQ=255)
                                       tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCCTTg                         >  1:19925065/1‑67 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAcc                                             <  1:4237563/45‑1 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCCCTCCCCAAGCGTTGCca                      >  1:13603508/1‑68 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                      >  1:23174184/1‑68 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                      >  1:10736490/1‑68 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                      >  1:19837475/1‑68 (MQ=255)
                                         cctGTGGGAGAGGGCCGGGGTGACGGCACCAACGCGCAGCGGGCCCATTCACTCCCCAAGCGATGCca                      >  1:24071939/1‑68 (MQ=25)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:21125832/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:28394966/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:14782845/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:356218/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:27282898/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:11155490/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:4874968/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:8293737/1‑67 (MQ=255)
                                             tgGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCcatca                   >  1:15942164/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:872657/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCCCTCCCCAAGCATTGCcctcata                 >  1:6378713/1‑66 (MQ=37)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:7266686/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:4594541/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:23148914/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:24732918/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:9050423/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:18202770/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:18269116/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:8942406/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:14211531/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:22303489/1‑67 (MQ=255)
                                               ggAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                 >  1:26176619/1‑67 (MQ=255)
                                                gagaGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                >  1:4980811/1‑67 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:19510860/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:18466264/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:740456/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:11333498/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:10640570/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:1817961/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:8968452/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:17775756/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:161708/1‑68 (MQ=255)
                                                 agagGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCAcc              >  1:15115342/1‑68 (MQ=255)
                                                   agGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCg             >  1:10230550/1‑67 (MQ=255)
                                                   agGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCg             >  1:25956579/1‑67 (MQ=255)
                                                    gggCGGGGGGGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:5383686/68‑1 (MQ=255)
                                                     ggCGGGGGTGGGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:9535034/67‑1 (MQ=255)
                                                     ggCGGGGGTGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:16297547/67‑1 (MQ=255)
                                                     ggCGGGGGTGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:26644292/67‑1 (MQ=255)
                                                     ggCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                      <  1:163334/56‑1 (MQ=255)
                                                      gggggggTGGGGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:9546520/64‑1 (MQ=255)
                                                      gCCGGGCGTCAGGCCCCCATCGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGcc           <  1:8934826/67‑1 (MQ=255)
                                                      gCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCtt         >  1:17271524/1‑68 (MQ=255)
                                                         ggggTGTGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTaa       <  1:18697993/67‑1 (MQ=255)
                                                         ggggTGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTaa       <  1:9419062/67‑1 (MQ=255)
                                                          gggTGTGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAAt      <  1:15067433/67‑1 (MQ=255)
                                                          gggTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAAt      <  1:23805159/67‑1 (MQ=255)
                                                          gggTGCGGCCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAAt      <  1:28663249/67‑1 (MQ=255)
                                                           ggTGAGGGCACCAACGCGCAGCGGACCGAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCg    >  1:5077361/1‑68 (MQ=255)
                                                           ggTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCCCCGCCTTAATCg    >  1:2638463/1‑68 (MQ=255)
                                                           ggTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCg    >  1:14868008/1‑68 (MQ=255)
                                                           ggTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCg    >  1:28799887/1‑68 (MQ=255)
                                                             gTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGGCTTAATCGt   >  1:11587324/1‑68 (MQ=255)
                                                              tGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCATTAATCGTa  >  1:13128398/1‑68 (MQ=255)
                                                                |                                                                  
GGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTA  >  CP009273/284923‑285051

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: