Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,554,355 G→T 8.6% intergenic (+98/‑372) yffL → / → yffM CPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,554,3550GT 8.6% 1043.8 / 8.5 401intergenic (+98/‑372)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (160/194);  new base T (10/36);  total (170/231)
Fisher's exact test for biased strand distribution p-value = 2.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 418 total aligned reads displayed.

TCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAA  >  CP009273/2554288‑2554421
                                                                   |                                                                  
tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCg                                                                    <  1:42878516/68‑1 (MQ=255)
tCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGGGGGTCAAAAGTTGGCg                                                                    >  1:27519152/1‑68 (MQ=255)
 cTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCga                                                                   <  1:13051245/68‑2 (MQ=255)
  tttGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCttt                                                                  >  1:2134416/1‑68 (MQ=255)
    tGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTaa                                                                <  1:7048211/68‑1 (MQ=255)
     gCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAg                                                               <  1:10386131/68‑1 (MQ=255)
      ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTAAAAAGTTGCAGTTAAgg                                                              >  1:38807271/1‑68 (MQ=25)
      ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAgg                                                              >  1:46698843/1‑68 (MQ=255)
      ccTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAATTTTGCCTTTAAgg                                                              >  1:36967218/1‑68 (MQ=25)
       cTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCTTTAAgga                                                             >  1:14898465/1‑67 (MQ=21)
          aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                          <  1:19968201/68‑1 (MQ=255)
          aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                          <  1:43549415/68‑1 (MQ=255)
          aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                          <  1:47115272/68‑1 (MQ=255)
          aaTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                                          <  1:50253949/68‑1 (MQ=255)
           aTTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                                         <  1:41912191/68‑1 (MQ=255)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                        >  1:26892735/1‑68 (MQ=255)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                        >  1:12435739/1‑68 (MQ=255)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                                        >  1:22276604/1‑68 (MQ=255)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAATTTGCCGTTAAGGTGGCTc                                                        >  1:35298795/1‑68 (MQ=21)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCCGTTAAGGTGGGTc                                                        >  1:25199296/1‑68 (MQ=255)
            ttCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTAAAAAGTTGCAGTTAAGGTGGGTc                                                        >  1:30135399/1‑68 (MQ=255)
             tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTATTGGTGGGTCAAAAGTTGCCGTTAAGGTGGCTCa                                                       >  1:32886029/1‑68 (MQ=255)
             tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTTGGAAAAAAGTTGCAGTAAAGGAGTGTCa                                                       >  1:5436255/1‑68 (MQ=21)
             tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                                       >  1:44793191/1‑68 (MQ=255)
             tCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCa                                                       >  1:22760329/1‑68 (MQ=255)
                 tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAg                                                   >  1:39209672/1‑68 (MQ=255)
                     ttCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTg                                                >  1:49776238/1‑67 (MQ=255)
                      tCAGGTGGGTCAAAATTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:30491998/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:3978743/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:37658991/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:10765069/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:20324310/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:7755068/68‑1 (MQ=255)
                      tCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGAGGGTCAAAAGTTGcc                                              >  1:31396427/1‑68 (MQ=255)
                      tCAGGTGGGTCAAAAGGTGCCGTTAGTGGGGGGTCAAAATATTCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:46975135/68‑1 (MQ=17)
                      tCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:4028833/68‑1 (MQ=255)
                      tCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:45834674/68‑1 (MQ=255)
                      tCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:34228425/68‑1 (MQ=255)
                      tCAGGTGGGGCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:9548506/68‑1 (MQ=255)
                      tCAGGGGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAATAGTTGcc                                              <  1:18208580/68‑1 (MQ=255)
                      tCAGGGGGGTCAAAAGTAGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGcc                                              <  1:8758664/68‑1 (MQ=255)
                         ggTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTGAAGGAGGGTCAAAAGTTGCAGtt                                           >  1:24619346/1‑68 (MQ=17)
                           tGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAATTTGCCGTTaa                                         >  1:50744923/1‑68 (MQ=255)
                           tGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGttta                                         >  1:8076740/1‑66 (MQ=255)
                              gTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAgg                                       <  1:35956334/67‑1 (MQ=255)
                                 aaaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                   <  1:19056111/68‑1 (MQ=255)
                                 aaaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTggg                                   <  1:23795967/68‑1 (MQ=255)
                                  aaaGTTGCCGTTAGTGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                  <  1:25776261/68‑1 (MQ=255)
                                  aaaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGt                                  <  1:19434257/68‑1 (MQ=255)
                                   aaGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTc                                 >  1:13513556/1‑68 (MQ=255)
                                    aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                >  1:29671082/1‑68 (MQ=255)
                                    aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                >  1:48199973/1‑68 (MQ=255)
                                    aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                >  1:28950335/1‑68 (MQ=255)
                                    aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCa                                >  1:6576247/1‑68 (MQ=255)
                                    aGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGACGATAAGGTGGGGCa                                >  1:24037979/1‑68 (MQ=18)
                                       tGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCaaaa                             <  1:21110140/68‑1 (MQ=255)
                                             tAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGc                        <  1:4945712/67‑1 (MQ=255)
                                               gtggtgGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAATGTGGGTCAAAAGTTGcata                     >  1:28354063/1‑65 (MQ=255)
                                                tggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAg                    >  1:38102645/1‑68 (MQ=255)
                                                 tgtggGTCATTAGTTGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                   <  1:26218677/67‑1 (MQ=17)
                                                 ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                   >  1:20703522/1‑68 (MQ=255)
                                                 ggtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCAGTTAAGGTGGGTCAAAAGTTGCCGAg                    <  1:23293805/67‑1 (MQ=255)
                                                 ggtgcGTCAAACGTTGCCGTTACGGTGGGGCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGAg                    <  1:37315217/67‑1 (MQ=14)
                                                  gtggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGt                   <  1:47821408/67‑1 (MQ=255)
                                                   tggGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTaa                 >  1:32171498/1‑68 (MQ=255)
                                                   tggGTCAAAAGGTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                  <  1:235346/67‑1 (MQ=255)
                                                   tggGTAAATTGTTGCCGTTAAGGTGTGTCAAAAGTTGCCGTTAAGGGGGGTCAAAAGTTGCGGAGTa                  <  1:51107953/67‑1 (MQ=17)
                                                   tggGTAAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTa                  <  1:6212264/67‑1 (MQ=255)
                                                     ggTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa               >  1:31680634/1‑68 (MQ=255)
                                                      gTCAAAAGTTGGCGTTTAGGCGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa               <  1:34869721/67‑1 (MQ=18)
                                                      gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa               <  1:43816427/67‑1 (MQ=255)
                                                      gTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGa               <  1:47437086/67‑1 (MQ=255)
                                                        cAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAt             >  1:44178978/1‑67 (MQ=255)
                                                        cAAAAGTAGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAAt             <  1:37202506/67‑1 (MQ=255)
                                                          aatgTTGCCTTTAAGGTGGGTCAATAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          <  1:29446448/65‑1 (MQ=21)
                                                          aaaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          <  1:1008916/68‑1 (MQ=255)
                                                          aaaGTAGCCGTTAAGGTGGGTCAACAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          <  1:33780127/68‑1 (MQ=255)
                                                          aaaGCTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          <  1:15448123/68‑1 (MQ=255)
                                                           aaGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          >  1:5842176/1‑67 (MQ=255)
                                                            acgtGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctac        <  1:17656609/65‑1 (MQ=255)
                                                            aGTTGCCGTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCATAAGTTGCCGAGTAAGAATTctac        <  1:3101031/68‑1 (MQ=21)
                                                            aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctac        <  1:21928267/68‑1 (MQ=255)
                                                            aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          >  1:2434078/1‑66 (MQ=255)
                                                            aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTct          >  1:40814308/1‑66 (MQ=255)
                                                            aGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGCGCAAAAGTTGACGAGTAAGAATTct          >  1:32486027/1‑66 (MQ=21)
                                                             gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       >  1:16445048/1‑68 (MQ=255)
                                                             gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       >  1:46155099/1‑68 (MQ=255)
                                                             gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       >  1:5135597/1‑68 (MQ=255)
                                                             gTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       >  1:26557091/1‑68 (MQ=255)
                                                             cttGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       <  1:43274882/67‑1 (MQ=255)
                                                               tGCCTTTTAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:42064619/68‑1 (MQ=255)
                                                               tGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:15906635/68‑1 (MQ=255)
                                                               tGCCTTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:42521011/68‑1 (MQ=255)
                                                               tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTctact       <  1:28730755/66‑1 (MQ=255)
                                                               tGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:9086085/68‑1 (MQ=255)
                                                               cgCCTTTAAGGTGGGTCAAAAGTTGCCCTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:16894477/67‑1 (MQ=21)
                                                               cgCCGTTAAGTTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTAt     <  1:39479242/67‑1 (MQ=255)
                                                                  cGTTACGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTa   <  1:49237957/67‑1 (MQ=255)
                                                                  cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTaa  >  1:48582702/1‑68 (MQ=255)
                                                                  cGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTa   <  1:7531281/67‑1 (MQ=255)
                                                                   |                                                                  
TCTTTGCCTTAATTCCATGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAA  >  CP009273/2554288‑2554421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: