Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,314 G→A 10.3% Q97Q (CAG→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,3140GA10.3% 1666.9 / 147.9 567Q97Q (CAG→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (293/213);  new base A (48/10);  total (342/225)
Fisher's exact test for biased strand distribution p-value = 1.85e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

Only 100 of 581 total aligned reads displayed.

CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381
                                                                   |                                                                   
cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                                                                     >  1:7621365/1‑68 (MQ=255)
  cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                                                                   >  1:36980699/1‑68 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:14589935/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            >  1:25777919/1‑68 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:29357111/68‑1 (MQ=255)
         gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                            <  1:340390/68‑1 (MQ=255)
             tcATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                         <  1:12559011/66‑1 (MQ=255)
                tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:16814090/1‑66 (MQ=255)
                tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:21810085/1‑66 (MQ=255)
                tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:29687609/1‑66 (MQ=255)
                 aCTAAGCATATGCTGCGCCTACGTCGGGTGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:24138420/1‑68 (MQ=255)
                 aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                       >  1:32282536/1‑65 (MQ=255)
                  cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:15943491/1‑67 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                                    >  1:13334363/1‑66 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:35581714/68‑1 (MQ=255)
                   tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:27071715/68‑1 (MQ=255)
                    aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                                  <  1:503113/64‑1 (MQ=255)
                    aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt                                                 >  1:22257882/1‑68 (MQ=255)
                     aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                                  <  1:35174342/66‑1 (MQ=255)
                     aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                                >  1:10268429/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:26891669/1‑68 (MQ=255)
                      gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:21653418/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:1628436/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:10807458/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:33582676/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:29167153/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:13153286/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:19405137/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:362725/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:305482/1‑68 (MQ=255)
                       cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                              >  1:24846290/1‑68 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:20019608/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:23403951/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:14947450/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:15116214/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:8560515/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:6699370/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:4021867/1‑67 (MQ=255)
                       cATATGCTGCGCCTACGTCGGCAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtaa                                               >  1:3862718/1‑65 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:35955362/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                               >  1:23718973/1‑66 (MQ=255)
                        atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCGTa                                               >  1:7090567/1‑66 (MQ=255)
                         tatGCTGCGCCTACGTGGGGAGGGGCAAATTACCGGTAATCAGGTTCCGGAAATTATCCTGCTTAAct                                            <  1:20778617/68‑1 (MQ=255)
                         tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                            <  1:30890128/68‑1 (MQ=255)
                          atGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           <  1:11056360/68‑1 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActct                                          <  1:27316215/68‑1 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:24399546/1‑67 (MQ=255)
                           tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                           >  1:8346460/1‑67 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:22469968/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:555491/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:25459747/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:11547232/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:36060846/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:33240702/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:18261734/1‑68 (MQ=255)
                            gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac                                         >  1:5735920/1‑66 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:28612384/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:15596212/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:19781431/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:18852042/1‑67 (MQ=255)
                             ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                         >  1:32557082/1‑67 (MQ=255)
                               gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                      <  1:18499026/68‑1 (MQ=255)
                                 gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                                    >  1:34966877/1‑68 (MQ=255)
                                  ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   >  1:6149790/1‑68 (MQ=255)
                                   cTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                                   <  1:20496099/67‑1 (MQ=255)
                                    tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                                 <  1:12342616/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:28922981/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:1003428/68‑1 (MQ=255)
                                     aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                                <  1:26003971/68‑1 (MQ=255)
                                      cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCCGAAATAATCCTGCTTAACTCTCACGATGGCa                                >  1:252142/1‑67 (MQ=255)
                                       gTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCAcc                              >  1:17261160/1‑68 (MQ=255)
                                         cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                            <  1:3191129/68‑1 (MQ=255)
                                          gggAGGTGCAAATTACCGGTAAACAGGTTCGGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:29487265/68‑1 (MQ=255)
                                          gggAGGGGCATATTACCGGTATAAAGGTTCCGGAATTTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:13582758/68‑1 (MQ=21)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:14305421/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:34300834/68‑1 (MQ=255)
                                          gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACACTCACGATGGCACCAGt                           <  1:9757184/68‑1 (MQ=255)
                                          gggACGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           <  1:34691368/68‑1 (MQ=255)
                                           ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                           >  1:16514371/1‑67 (MQ=255)
                                            gAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTc                         >  1:26651254/1‑68 (MQ=255)
                                             aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                        >  1:4719760/1‑68 (MQ=255)
                                              ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       >  1:9313480/1‑68 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:28162932/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:17938620/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:4574447/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                       <  1:23036097/66‑1 (MQ=255)
                                                ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                      <  1:36522308/67‑1 (MQ=255)
                                                 gCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTAc                     <  1:31470148/67‑1 (MQ=255)
                                                   aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg                  <  1:26424325/68‑1 (MQ=255)
                                                   aaaTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAg                  <  1:19155021/68‑1 (MQ=255)
                                                       tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgc              <  1:33448114/68‑1 (MQ=255)
                                                        aCCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgct             <  1:24604743/68‑1 (MQ=255)
                                                         ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg            <  1:22034622/68‑1 (MQ=255)
                                                         ccGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgctg            <  1:19074332/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:27840472/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:11878761/68‑1 (MQ=255)
                                                              aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg       <  1:10482520/68‑1 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     <  1:30234445/68‑1 (MQ=255)
                                                                aCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCa     <  1:8876275/68‑1 (MQ=255)
                                                                   ggTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGt  >  1:3303204/1‑68 (MQ=255)
                                                                   |                                                                   
CCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGGCATGT  >  CP009273/2762247‑2762381

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: