Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,288 C→A 8.4% R89R (CGG→AGG)  yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2880CA 8.4% 2197.3 / 136.0 628R89R (CGG→AGG) yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (431/143);  new base A (38/15);  total (469/159)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 640 total aligned reads displayed.

GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355
                                                                  |                                                                   
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:31634339/1‑68 (MQ=255)
 ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                                   >  1:33716458/1‑68 (MQ=255)
     tCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAgg                                                               >  1:5365458/1‑68 (MQ=255)
      cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                              <  1:36111433/68‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                                            <  1:18556785/66‑1 (MQ=255)
          cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                                          >  1:38174597/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:22541279/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:15207833/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         >  1:13134951/1‑68 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         <  1:6708865/68‑1 (MQ=255)
           acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                         <  1:7637894/68‑1 (MQ=255)
            cGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAt                                                        <  1:4523827/68‑1 (MQ=255)
                 aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                                    <  1:9134562/67‑1 (MQ=255)
                   gTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                 <  1:2779193/68‑1 (MQ=255)
                    tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                >  1:15468172/1‑68 (MQ=255)
                     gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:17826409/1‑68 (MQ=255)
                      aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:23634805/1‑67 (MQ=255)
                      aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                               >  1:26642169/1‑67 (MQ=255)
                       cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                             >  1:19756477/1‑68 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:34714029/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:32656537/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:31917439/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:20368943/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:21548920/68‑1 (MQ=255)
                                  gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                  <  1:30149167/68‑1 (MQ=255)
                                   tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                 >  1:5803357/1‑68 (MQ=255)
                                    gAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAt                                <  1:21197241/68‑1 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                               <  1:7941152/67‑1 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                              >  1:35192212/1‑68 (MQ=255)
                                      gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                              >  1:4846397/1‑68 (MQ=255)
                                        aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             <  1:217453/67‑1 (MQ=255)
                                         tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:24732509/1‑66 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:11611314/1‑65 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                             >  1:23891082/1‑65 (MQ=255)
                                          aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          <  1:4210693/68‑1 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:10123082/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                          >  1:24162564/1‑67 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTg                         <  1:17212806/68‑1 (MQ=255)
                                           cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTAGCCt                          >  1:27040984/1‑67 (MQ=255)
                                            tAAGCATATGCTGCGCCTACGGCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                        <  1:18309472/68‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:38657030/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                        <  1:1346438/64‑1 (MQ=255)
                                             aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                       >  1:17929572/4‑68 (MQ=255)
                                             aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                       >  1:13932463/1‑68 (MQ=255)
                                             aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCt                       >  1:35414272/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:18116929/1‑68 (MQ=255)
                                               gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:34112183/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:11459902/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:22011970/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:3997369/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:8282571/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:12311248/1‑68 (MQ=255)
                                                cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                    >  1:712474/1‑68 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:12880185/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:873338/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:3892408/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:37543677/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:6107103/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:11138135/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:9624550/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:28198606/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:28434766/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:26344559/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:24519543/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:19393385/1‑67 (MQ=255)
                                                cATATGCTGCGCCTACGTCGGGAGGCGCAAATTACCGGTAAACAGGATCCGGAAAATATCCAGCttt                     >  1:20847657/1‑66 (MQ=25)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAATTTATCCTGCTTa                     >  1:10812979/1‑66 (MQ=255)
                                                 atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                     >  1:36974872/1‑66 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:32420721/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:11986069/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:34244580/1‑67 (MQ=255)
                                                    tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                 >  1:25271490/1‑67 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:20033362/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:12584660/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:29148810/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:14801249/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:34965355/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:29536983/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:1746369/1‑68 (MQ=255)
                                                     gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGATCCGGAAATTATCCTGCTTAACtctc               >  1:32165073/1‑68 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTaaa                                              >  1:15902398/1‑36 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTATCCTGCTGAACtctc               >  1:16227905/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:14309690/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:22907423/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:27376120/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:38453663/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:30996047/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:29298502/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:16818229/1‑67 (MQ=255)
                                                      ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc               >  1:30494063/1‑67 (MQ=255)
                                                       tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                  >  1:8472195/1‑63 (MQ=255)
                                                        gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg            <  1:19101686/68‑1 (MQ=255)
                                                        gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg            <  1:37937187/68‑1 (MQ=255)
                                                           ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:2603640/1‑68 (MQ=255)
                                                           ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:33060804/1‑68 (MQ=255)
                                                            cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg         >  1:15695252/1‑67 (MQ=255)
                                                             tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc       <  1:35885234/68‑1 (MQ=255)
                                                               cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGcc      >  1:10014583/1‑66 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  <  1:2559142/68‑1 (MQ=255)
                                                                  cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg  >  1:36498219/1‑68 (MQ=255)
                                                                  |                                                                   
GCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAG  >  CP009273/2762222‑2762355

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: