Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,299 A→T 10.3% Q92H (CAA→CAT yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2990AT10.3% 1790.3 / 142.8 578Q92H (CAA→CATyfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (328/190);  new base T (48/12);  total (376/202)
Fisher's exact test for biased strand distribution p-value = 9.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 608 total aligned reads displayed.

CACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGA  >  CP009273/2762232‑2762366
                                                                   |                                                                   
 acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                    >  1:2051090/1‑68 (MQ=255)
 acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                    >  1:22380958/1‑68 (MQ=255)
 acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                    >  1:12991322/1‑68 (MQ=255)
 acGGGTACGTGACCCCCGTCGTAGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                                                    >  1:24001203/1‑68 (MQ=255)
        cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                                             <  1:7762291/68‑1 (MQ=255)
         gagACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                            >  1:7301453/3‑68 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                            <  1:4450918/68‑1 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                            >  1:14388433/1‑68 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                            >  1:1496588/1‑68 (MQ=255)
         gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                                            <  1:10883701/68‑1 (MQ=255)
          tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                           >  1:13190941/1‑68 (MQ=255)
          tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                           >  1:27430538/1‑68 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                          >  1:6096671/1‑68 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                          >  1:31900032/1‑68 (MQ=255)
           gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                          >  1:6554061/1‑68 (MQ=255)
            aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                          >  1:25150127/1‑67 (MQ=255)
               cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                      >  1:29984880/1‑68 (MQ=255)
               cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                      >  1:23248266/1‑68 (MQ=255)
                 cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                                    >  1:33353133/1‑68 (MQ=255)
                 cgtcgtcgtGAGCATACTAAGCATATGCTGCCCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                                    <  1:32973487/68‑1 (MQ=255)
                   tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                   >  1:17419008/1‑67 (MQ=255)
                    cgtcgtGAGCATACTAAGCATATGCGGAGCCTACGTCGGGAGGGGGAAATTACCGGTAAACAGGTTcc                                                 >  1:34260853/1‑68 (MQ=25)
                        gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                             >  1:25777919/1‑68 (MQ=255)
                         tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCATATTACCGGTAAACAGGTTCCGGaa                                             <  1:13601350/67‑1 (MQ=255)
                         tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                                            >  1:19390086/1‑68 (MQ=255)
                            tcATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                          <  1:12559011/66‑1 (MQ=255)
                            gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                          <  1:19525569/67‑1 (MQ=255)
                            gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                         >  1:33867895/1‑68 (MQ=255)
                            gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGATCCGGAAATTa                                         >  1:19106866/1‑68 (MQ=255)
                              aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                        <  1:20162177/67‑1 (MQ=255)
                               tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                        >  1:21810085/1‑66 (MQ=255)
                               tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                        >  1:16814090/1‑66 (MQ=255)
                                aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                        >  1:31113095/1‑65 (MQ=255)
                                  tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                     >  1:1580181/1‑66 (MQ=255)
                                  tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:4076359/68‑1 (MQ=255)
                                  tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:18105882/68‑1 (MQ=255)
                                  tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:8695773/68‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                   <  1:503113/64‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                   <  1:10612359/64‑1 (MQ=255)
                                   aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                  >  1:25478932/4‑68 (MQ=255)
                                    aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                   <  1:35174342/66‑1 (MQ=255)
                                     gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:21653418/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:24846290/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:11723600/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:14660618/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:1628436/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:35655160/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:36052163/1‑68 (MQ=255)
                                      cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                               >  1:18397998/1‑68 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:6974892/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:17029860/1‑67 (MQ=255)
                                      cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:30470906/1‑67 (MQ=255)
                                       atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:23718973/1‑66 (MQ=255)
                                       atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:32257828/1‑66 (MQ=255)
                                       atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                >  1:9701255/1‑66 (MQ=255)
                                        tatGCTGCGCCTACGTGGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTATCCTGCTTAAct                             <  1:20778617/68‑1 (MQ=255)
                                        tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                             <  1:8404659/68‑1 (MQ=255)
                                        aatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTCCTTAAc                              >  1:29423744/2‑67 (MQ=255)
                                         atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            <  1:33553168/68‑1 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:10230942/1‑67 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:35539236/1‑67 (MQ=255)
                                          tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                            >  1:24356355/1‑67 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:938151/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:31585970/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:555491/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:11084528/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:3464985/1‑68 (MQ=255)
                                           gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:36441355/1‑68 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACCGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:15379430/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:29571947/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:22851683/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:2461882/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:20905274/1‑67 (MQ=255)
                                            ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                          >  1:28612384/1‑67 (MQ=255)
                                             tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCa                         <  1:13356642/67‑1 (MQ=255)
                                             tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                        >  1:9028214/1‑68 (MQ=255)
                                              gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                       <  1:34888839/68‑1 (MQ=255)
                                              gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                       <  1:22203588/68‑1 (MQ=255)
                                                gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                     <  1:3904895/68‑1 (MQ=255)
                                                gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                     >  1:27087110/1‑68 (MQ=255)
                                                   tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                  <  1:23482886/68‑1 (MQ=255)
                                                   tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                  <  1:26174265/68‑1 (MQ=255)
                                                   tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                  <  1:12342616/68‑1 (MQ=255)
                                                    aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 <  1:1003428/68‑1 (MQ=255)
                                                    aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 <  1:28922981/68‑1 (MQ=255)
                                                     cgacGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACACTCACGATGGCa                 >  1:26996696/4‑67 (MQ=255)
                                                     cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 >  1:36800643/1‑67 (MQ=255)
                                                     cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                 >  1:16020129/1‑67 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:37071597/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:19814948/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:5764617/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:26506878/68‑1 (MQ=255)
                                                         gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            <  1:32519709/68‑1 (MQ=255)
                                                          ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt            >  1:16514371/1‑67 (MQ=255)
                                                            aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCCCCAGTTCg         >  1:2923817/1‑68 (MQ=255)
                                                            aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg         >  1:4719760/1‑68 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:36262659/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:26771139/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt        <  1:28162932/66‑1 (MQ=255)
                                                               ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa       <  1:18873064/67‑1 (MQ=255)
                                                                   |                                                                   
CACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGA  >  CP009273/2762232‑2762366

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: