Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,309 A→C 10.3% 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,3090AC10.3% 1770.1 / 158.9 550K96Q (AAA→CAA) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (296/196);  new base C (47/10);  total (344/206)
Fisher's exact test for biased strand distribution p-value = 7.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 588 total aligned reads displayed.

TGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCGGG  >  CP009273/2762242‑2762376
                                                                   |                                                                   
tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                                                     >  1:27430538/1‑68 (MQ=255)
 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTGGGGAGGGGCAAATAACCGGTaa                                                                    >  1:22636264/1‑68 (MQ=255)
 gACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                    >  1:6828903/1‑68 (MQ=255)
  aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                                                    >  1:4145281/1‑67 (MQ=255)
      ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                                                               <  1:20593609/68‑1 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:31167317/1‑67 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:17419008/1‑67 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:28027358/1‑67 (MQ=255)
         tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                                                             >  1:4670145/1‑67 (MQ=255)
          cgtcgtGAGCATACTAAGCAAATGCTGCGCCTACCTCGGGAGGGGCAAATTACCGGTAAACAGGTTc                                                            >  1:24811767/1‑67 (MQ=255)
              gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                                                       <  1:18418446/68‑1 (MQ=255)
               tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCATATTACCGGTAAACAGGTTCCGGaa                                                       <  1:13601350/67‑1 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                                                    <  1:25854650/67‑1 (MQ=255)
                  gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                                                   >  1:29090370/1‑68 (MQ=255)
                    aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACACGTTCCGGAAATTa                                                   <  1:36143109/66‑1 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                                                  >  1:32282536/1‑65 (MQ=255)
                      aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               <  1:29268407/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:13334363/1‑66 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt                                               >  1:35617877/1‑66 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:26802737/68‑1 (MQ=255)
                        tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc                                             <  1:34167215/68‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:10612359/64‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc                                             <  1:24586760/64‑1 (MQ=255)
                         aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCt                                            >  1:18146492/4‑68 (MQ=255)
                          aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt                                           >  1:10268429/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:19142116/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:3386587/1‑68 (MQ=255)
                           gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:28770668/1‑68 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          >  1:35069360/2‑68 (MQ=255)
                           acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa                                          >  1:28461694/2‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:37088402/1‑68 (MQ=255)
                            cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa                                         >  1:25038051/1‑68 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTATCCTGCTTa                                          >  1:7473879/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:26310310/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:8560515/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:29424587/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:13176436/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:5965967/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:16578234/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:19400578/1‑67 (MQ=255)
                            cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:16989675/1‑67 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCTAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:30666162/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:23718973/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:33394033/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:25617469/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa                                          >  1:3913387/1‑66 (MQ=255)
                             atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCTTa                                          >  1:1963288/1‑66 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:12936914/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct                                       <  1:3639688/68‑1 (MQ=255)
                              tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc                                        >  1:24252180/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:23438089/1‑67 (MQ=255)
                                tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc                                      >  1:8346460/1‑67 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:5676544/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:11084528/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:31585970/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:15458161/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:15816364/1‑68 (MQ=255)
                                 gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCGGCTTAACtctc                                    >  1:36017455/1‑68 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:32557082/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:11595461/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:5085057/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:29572030/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:36717156/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc                                    >  1:20905274/1‑67 (MQ=255)
                                  ctgcGCCTACGTCGGGAGGGGCAAAATACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACactc                                    >  1:21835511/1‑67 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                  >  1:2018962/1‑68 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                  >  1:6315253/1‑68 (MQ=255)
                                   tgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCAc                                  >  1:24008604/1‑68 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:14465818/68‑1 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:18935074/68‑1 (MQ=255)
                                    gcgcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACg                                 <  1:34880528/68‑1 (MQ=255)
                                      gcCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGAt                               >  1:34573297/1‑68 (MQ=255)
                                       ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                              >  1:20001341/1‑68 (MQ=255)
                                       ccTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATg                              >  1:8950224/1‑68 (MQ=255)
                                         tACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGc                            <  1:21734208/68‑1 (MQ=255)
                                          aCGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           <  1:1003428/68‑1 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCCCGATGGCa                           >  1:4932492/1‑67 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           >  1:26531393/1‑67 (MQ=255)
                                           cGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCa                           >  1:16020129/1‑67 (MQ=255)
                                              cGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAg                       <  1:3191129/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:11970689/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:30308297/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      <  1:27155470/68‑1 (MQ=255)
                                               gggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACACTCACGATGGCACCAGt                      <  1:9757184/68‑1 (MQ=255)
                                                ggAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGt                      >  1:27036776/1‑67 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTATCTCTCACGATGGCACCAGTTCg                   >  1:16380786/1‑68 (MQ=255)
                                                  aGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCg                   >  1:22182597/1‑68 (MQ=255)
                                                   ggggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  >  1:12468666/1‑68 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:25319879/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGt                  <  1:14974772/66‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                 <  1:35490350/67‑1 (MQ=255)
                                                     ggCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTa                 <  1:934027/67‑1 (MQ=255)
                                                          aTTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt           >  1:7060932/1‑68 (MQ=255)
                                                           ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt           >  1:22258817/1‑67 (MQ=255)
                                                           ttACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAt           >  1:33174182/1‑67 (MQ=255)
                                                            tACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGAtgc         <  1:33491095/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:14620273/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:10880517/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:18675501/68‑1 (MQ=255)
                                                                   aaaCAGGTTCCGGAAATTATCCTGCTTAACTCTCACGATGGCACCAGTTCGTACCAGATGCTGCCggg  <  1:23079330/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: