Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,762,276 C→T 6.0% 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,2760CT 6.0% 1402.1 / 53.9 421L85L (CTG→TTG) ‡yfjQCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (243/150);  new base T (17/8);  total (262/159)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 563 total aligned reads displayed.

CCAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACTCTC  >  CP009273/2762209‑2762342
                                                                   |                                                                  
 cAGCCATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATAtgct                                                                   <  1:28354747/68‑1 (MQ=255)
     cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCTTATGCTgcg                                                                >  1:19849293/1‑67 (MQ=255)
     cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcg                                                                <  1:14834996/67‑1 (MQ=255)
     cATTCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTgcg                                                                <  1:10005801/67‑1 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                            <  1:10254139/68‑1 (MQ=255)
        tCTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTa                                                            <  1:31961446/68‑1 (MQ=255)
         cTTTGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                           <  1:429534/68‑1 (MQ=255)
          tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                          <  1:18578062/68‑1 (MQ=255)
          tttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                          <  1:23531530/68‑1 (MQ=255)
          tttGCCAGTCAGACACGGGTACGTGACCCCCGGCGTCGTGAGCATACTAAGCATATGCTGCGCCTACg                                                          <  1:20178447/68‑1 (MQ=255)
           ttGCCTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGt                                                         <  1:25118929/68‑1 (MQ=255)
              ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      >  1:26236776/1‑68 (MQ=255)
              ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      >  1:10534563/1‑68 (MQ=255)
              ccTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      >  1:18021883/1‑68 (MQ=255)
               cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCggg                                                     >  1:27638953/1‑68 (MQ=255)
               cTGTCAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCgg                                                      <  1:20884837/67‑1 (MQ=255)
                   cAGACACGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAggg                                                 <  1:1922489/68‑1 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                               <  1:23784887/66‑1 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGc                                               <  1:12411139/66‑1 (MQ=255)
                       cacGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaa                                             >  1:22379234/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:24514287/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:2959820/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:30353269/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:225546/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:31705856/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:6825810/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            >  1:10733509/1‑68 (MQ=255)
                        acGGGTACGTGACCCCCGTCGTCGTGACCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCaaa                                            <  1:14561942/68‑1 (MQ=255)
                           ggTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAAtt                                          >  1:202586/1‑67 (MQ=255)
                            gTACGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTa                                         >  1:9830526/1‑67 (MQ=255)
                              aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTAc                                                           <  1:13682562/47‑1 (MQ=255)
                              aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                       <  1:17007509/67‑1 (MQ=255)
                              aCGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTAcc                                       <  1:3786955/67‑1 (MQ=255)
                               cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                     <  1:22124892/68‑1 (MQ=255)
                               cGTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCgg                                     <  1:16028364/68‑1 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCggg                                    >  1:27856056/1‑67 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    <  1:26412415/68‑1 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    >  1:9304331/1‑68 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    >  1:1571467/1‑68 (MQ=255)
                                gTGACCCCCGTCGTCGTGAGCATACTAAGCAAATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGt                                    <  1:28672687/68‑1 (MQ=255)
                                 tGACCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTa                                   >  1:10891220/1‑68 (MQ=255)
                                   aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                  >  1:24049204/1‑67 (MQ=255)
                                   aCCCCCGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTaa                                  >  1:5487331/1‑67 (MQ=255)
                                    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                >  1:27958568/1‑68 (MQ=255)
                                    cccccGTCGTCGTGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAAc                                >  1:31269068/1‑68 (MQ=255)
                                        cgtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGt                            >  1:21334682/1‑68 (MQ=255)
                                         gtcgtcgtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGtt                           <  1:26719593/68‑1 (MQ=255)
                                               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                     <  1:25862134/68‑1 (MQ=255)
                                               gtGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaa                     <  1:574117/68‑1 (MQ=255)
                                                tGAGCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGaaa                    >  1:21498323/1‑68 (MQ=255)
                                                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAAtt                  <  1:5336729/67‑1 (MQ=255)
                                                   gCATACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTa                 >  1:3426250/1‑68 (MQ=255)
                                                     aTACTAATCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                >  1:4665783/1‑67 (MQ=255)
                                                      tACTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTAt                >  1:7204489/1‑66 (MQ=255)
                                                       aCTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             <  1:15487433/68‑1 (MQ=255)
                                                        cTAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:25670969/1‑67 (MQ=255)
                                                         tAATCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:27171983/68‑1 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:9143131/1‑66 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:18921601/1‑66 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCt             >  1:12685527/1‑66 (MQ=255)
                                                         tAAGCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGc           <  1:25493962/68‑1 (MQ=255)
                                                          aaacATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGc           <  1:1180728/64‑1 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:17436054/1‑68 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:8155729/1‑68 (MQ=255)
                                                            gCATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:17748006/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:13271773/1‑68 (MQ=255)
                                                             cATATGTTACGTCTTCGCAGGGAAGGGCATATTAACGGGCAGGAAGTTCCTGAAATTATCCTGCTTaa       >  1:29833905/1‑68 (MQ=255)
                                                             cATATGGTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:23233669/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:6244944/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:13825825/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:3609000/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:3346192/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:11209283/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:11572183/1‑67 (MQ=255)
                                                             cATATGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTa        >  1:22925856/1‑67 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc      <  1:7388428/68‑1 (MQ=255)
                                                              atatGCTGCGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAAc      <  1:20645720/68‑1 (MQ=255)
                                                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc    >  1:7602667/1‑67 (MQ=255)
                                                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc    >  1:16149932/1‑67 (MQ=255)
                                                                 tgctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAActc    >  1:25366935/1‑67 (MQ=255)
                                                                 tgctgcGCCTACGTCGGGAGGGGCAAATAACCCGAAAACAGGTACCGGAAAATATCCACCTTAActa    >  1:19483885/1‑66 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:2923695/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:18235166/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:8568365/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:4805348/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:3908229/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:30719007/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:29368616/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:26916616/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:14199512/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:15033412/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:164305/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:20500494/1‑68 (MQ=255)
                                                                  gctgcGCCTACGTCGGGAGGGGCAAATTACCGGAAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:24890791/1‑68 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTACCtctc  >  1:6026439/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:28126751/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:6532716/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:16639342/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTACCGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:14134667/1‑67 (MQ=255)
                                                                   ctgcGCCTACGTCGGGAGGGGCAAATTAACGGTAAACAGGTTCCGGAAATTATCCTGCTTAACtctc  >  1:19679265/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: