Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,276 C→T 12.9% L85L (CTG→TTG) ‡ yfjQ → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,762,2760CT12.9% 1232.7 / 152.5 445L85L (CTG→TTG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (246/140);  new base T (47/10);  total (295/150)
Fisher's exact test for biased strand distribution p-value = 4.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 539 total aligned reads displayed.

CCAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTC  >  CP009273/2762209‑2762342
                                                                   |                                                                  
     cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                 <  1:35033666/66‑1 (MQ=255)
     cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgctgc                                                                 <  1:13125236/66‑1 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                             >  1:1025966/1‑67 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                             >  1:22336950/1‑67 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCt                                                             >  1:28919315/1‑67 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                            <  1:7007191/68‑1 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                            <  1:36624722/68‑1 (MQ=255)
          tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                          <  1:29197408/68‑1 (MQ=255)
          tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCCCCTACg                                                          >  1:20940317/1‑68 (MQ=255)
          tttGCCTGTCAGACACGGGTACGTGACCCCCGTCATCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                          <  1:17019363/68‑1 (MQ=255)
              ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      >  1:18808922/1‑68 (MQ=255)
              ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      >  1:19953562/1‑68 (MQ=255)
               cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      <  1:27901665/67‑1 (MQ=255)
                    aGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                 <  1:6558230/67‑1 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                             <  1:27533089/68‑1 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                             >  1:21497031/1‑68 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                             >  1:6076826/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTGGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:35253383/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            <  1:10078144/68‑1 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            <  1:161930/68‑1 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:8967265/1‑68 (MQ=255)
                           ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGTGAGGGGCAAAtt                                          >  1:5388421/1‑67 (MQ=255)
                           ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                          >  1:965517/1‑67 (MQ=255)
                              aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                       <  1:29982088/67‑1 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    <  1:16820201/68‑1 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCAAATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    >  1:28419807/1‑68 (MQ=255)
                                 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                   >  1:36864570/1‑68 (MQ=255)
                                 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                   >  1:2390368/1‑68 (MQ=255)
                                      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                              >  1:23248266/1‑68 (MQ=255)
                                      cccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAg                              >  1:33431269/1‑68 (MQ=255)
                                       ccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAgg                             <  1:20593609/68‑1 (MQ=255)
                                          tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                           >  1:31167317/1‑67 (MQ=255)
                                          tcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                           >  1:17419008/1‑67 (MQ=255)
                                           cgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTcc                         >  1:7927543/1‑68 (MQ=255)
                                               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                     <  1:20849463/68‑1 (MQ=255)
                                               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                     <  1:13489389/68‑1 (MQ=255)
                                               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                     <  1:14589935/68‑1 (MQ=255)
                                                   tcATACTAAGCATATGCTGCGCCTACGTCGGGAGTGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                  <  1:36388259/66‑1 (MQ=255)
                                                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                 >  1:13986759/1‑68 (MQ=255)
                                                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                 >  1:18273587/1‑68 (MQ=255)
                                                     aTACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACACGTTCCGGAAATTa                 <  1:36143109/66‑1 (MQ=255)
                                                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                >  1:21810085/1‑66 (MQ=255)
                                                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                >  1:4417212/1‑66 (MQ=255)
                                                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                >  1:20346590/1‑65 (MQ=255)
                                                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:15943491/1‑67 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:35617877/1‑66 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:13227282/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:24665125/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:34167215/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:32255920/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:10996108/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:8587115/68‑1 (MQ=255)
                                                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc           <  1:3939147/64‑1 (MQ=255)
                                                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc           <  1:10612359/64‑1 (MQ=255)
                                                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc           <  1:24586760/64‑1 (MQ=255)
                                                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt          >  1:30853483/1‑68 (MQ=255)
                                                          aaGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCt          >  1:25439470/1‑68 (MQ=255)
                                                           aGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCtt         >  1:9480107/1‑68 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:27072766/1‑68 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:3386587/1‑68 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:32782061/1‑68 (MQ=255)
                                                            acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa        >  1:14813110/2‑68 (MQ=255)
                                                            acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa        >  1:31678122/2‑68 (MQ=255)
                                                            acATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTa        >  1:5042458/2‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:32363425/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:11723600/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:24846290/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:25038051/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:33582676/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:22511800/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:10807458/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:15594930/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:4716601/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTCaa       >  1:33178475/1‑68 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGATATTATCCTGCTTa        >  1:7473879/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:16578234/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:37121378/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:23403951/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:23568613/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:19400578/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:25865058/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:15116214/1‑67 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:25617469/1‑66 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:1908363/1‑66 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc      <  1:18093018/68‑1 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCAGCTTa        >  1:1963288/1‑66 (MQ=255)
                                                               tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct     <  1:12936914/68‑1 (MQ=255)
                                                               tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAct     <  1:8404659/68‑1 (MQ=255)
                                                               tatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc      >  1:24252180/1‑67 (MQ=255)
                                                               aatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTCCTTAAc      >  1:29423744/2‑67 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:3464985/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:26430725/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCGGCTTAACtctc  >  1:36017455/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAACCTGCTTAACactc  >  1:2687985/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAAATACCGGAAAACCGGTTCCGGAAATTAACCTGCTTAActcac  >  1:5735920/1‑66 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:22009054/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:229569/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:12553667/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:6850811/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:29572030/1‑67 (MQ=255)
                                                                   |                                                                  
CCAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTC  >  CP009273/2762209‑2762342

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: