Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 548,535 T→A intergenic (‑257/‑332) nmpC ← / → ybcR outer membrane porin, at locus of qsr prophage; b0553_1/predicted phage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606548,5350TA98.1% 672.6 / 2.0 161intergenic (‑257/‑332)nmpC/ybcRouter membrane porin, at locus of qsr prophage; b0553_1/predicted phage lysis protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base A (73/85);  total (74/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 184 total aligned reads displayed.

ATAAACCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAATCAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTGTGGTGGTGTT  >  REL606/548462‑548607
                                                                         |                                                                        
aTAAACCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACa                                                                        <  2:5568020/76‑1 (MQ=255)
aTAAACCTATTAGCTAATTTATATTAATGGCTGGTATTTATGGGGGGTCTATAATTTAGTGGTTTAATTTTAAACa                                                                        >  2:6214912/1‑76 (MQ=255)
 tAAACCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACaa                                                                       <  2:2441116/76‑1 (MQ=255)
  aaaCCTATTAGCTAATTTATATTAATGGCTGTTATTTTTGGGGGGTCTATAATTTAGTGGGTTAATTTTAAACaaa                                                                      >  2:3758212/1‑75 (MQ=255)
  aaaCCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAaacaac                                                                      <  2:2386324/76‑1 (MQ=255)
  aaaCCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAaacaac                                                                      <  1:2310148/76‑1 (MQ=255)
   aaCCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACt                                                                     <  1:775460/76‑1 (MQ=255)
   aaCCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACt                                                                     <  1:1324231/76‑1 (MQ=255)
      cTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTaaa                                                                  >  1:4977119/1‑76 (MQ=255)
      cTATTAGCTAATTTATATTAATGGCTGGTTTTTTTGGGGGGTTTATAATTTAGTGGGTTTATTTTAAAAAACTaaa                                                                  >  1:4661145/1‑76 (MQ=255)
      cTATTAACTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTaaa                                                                  <  1:7289641/76‑1 (MQ=255)
       tATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTaaaa                                                                 >  1:1403351/1‑76 (MQ=255)
       tATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTaaaa                                                                 <  2:4064216/76‑1 (MQ=255)
        aTTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTaaaaa                                                                <  1:7241167/76‑1 (MQ=255)
        aTTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTaaaaa                                                                >  1:5804049/1‑76 (MQ=255)
         ttAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAAt                                                               >  1:1807459/1‑76 (MQ=255)
           aGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATaa                                                             <  1:265636/76‑1 (MQ=255)
           aGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATaa                                                             <  1:4618151/76‑1 (MQ=255)
           aGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATaa                                                             <  2:2556072/76‑1 (MQ=255)
           aGCTAATTTATATTAATGGCTGTTATTTATGGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATaa                                                             >  1:6535308/1‑76 (MQ=255)
            gCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAAc                                                            <  2:153787/76‑1 (MQ=255)
             cTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACg                                                           >  1:6528630/1‑76 (MQ=255)
             cTAATTTATATTAATGGCTGTTATTTATGGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACg                                                           >  1:4076666/1‑76 (MQ=255)
              tAATTTATATTAATGGCTGTTATTTATGGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGc                                                          >  1:2413112/1‑76 (MQ=255)
              tAATTTATATTAATGGCTGGTTTTTTTTGGGGGTCTTTAATTTCGTGGTTTTATTTTAAACAAACAAAAATAAAGc                                                          >  2:6697959/1‑76 (MQ=255)
               aaTTTATATTAATGGCTGTTTTTTTTGGGGGGTTTATAATTCCGTGGGTTTATTTTAAACAACTAAAAATAAAGcc                                                         >  2:1293369/1‑76 (MQ=255)
               aaTTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGcc                                                         >  1:1034020/1‑76 (MQ=255)
               aaTTTATATTAATGGCTGGTTTTTTTGGGGGGTTTATAATTCAGTGGGTTAATTTAAAACAAATAAAAATAAAGcc                                                         >  2:790457/1‑76 (MQ=255)
                 tttATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCgg                                                       <  1:746159/76‑1 (MQ=255)
                 tttATATTAATGGCTGGTTTTTTTGGGGGGTCTTTAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCgg                                                       >  2:3554238/1‑76 (MQ=255)
                  ttATATTAATGGCTGGTATTTTTTGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACACCGGa                                                      >  2:5523487/1‑76 (MQ=255)
                   tataTTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGaa                                                     <  1:898656/76‑1 (MQ=255)
                   tataTTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGaa                                                     <  2:1210287/76‑1 (MQ=255)
                   tataTTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGaa                                                     <  2:6467306/76‑1 (MQ=255)
                   tataTTAATGGCTGGTTTTTTTGGGGGGTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAAAGCCGGaa                                                     >  2:6918596/1‑76 (MQ=255)
                    atatTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGaaa                                                    >  1:2795777/1‑76 (MQ=255)
                    atatTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGaaa                                                    >  1:119346/1‑76 (MQ=255)
                      atTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAatt                                                  <  2:4807639/76‑1 (MQ=255)
                      atTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAatt                                                  <  2:7130800/76‑1 (MQ=255)
                      atTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAatt                                                  <  2:2242308/76‑1 (MQ=255)
                       ttAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAatta                                                 <  1:4449320/76‑1 (MQ=255)
                          aTGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattt                                              <  2:5093665/76‑1 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             <  1:6266891/76‑1 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             <  1:1361355/76‑1 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             >  1:5394791/1‑76 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             >  2:6501193/1‑76 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             >  1:1076885/1‑76 (MQ=255)
                           tGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttattta                                             >  2:3055224/1‑76 (MQ=255)
                            ggCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttatttat                                            <  1:6650472/76‑1 (MQ=255)
                            ggCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAAttatttat                                            <  1:3951672/76‑1 (MQ=255)
                             gCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATtatttatt                                           >  2:3952329/1‑76 (MQ=255)
                             gCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATtatttatt                                           <  1:1413878/76‑1 (MQ=255)
                              cTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTg                                          <  1:7441254/76‑1 (MQ=255)
                               tGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTgg                                         >  1:7212459/1‑76 (MQ=255)
                               tGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTgg                                         >  1:5234840/1‑76 (MQ=255)
                                gTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGt                                        <  1:6595892/76‑1 (MQ=255)
                                 ttatttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGtt                                       <  1:6890512/76‑1 (MQ=255)
                                 ttatttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGtt                                       <  1:1613848/76‑1 (MQ=255)
                                   atttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAt                                     >  1:2464739/1‑76 (MQ=255)
                                   atttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAt                                     >  2:429935/1‑76 (MQ=255)
                                   atttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAt                                     >  1:1980535/1‑76 (MQ=255)
                                    tttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAtt                                    <  2:4992497/76‑1 (MQ=255)
                                     ttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAttt                                   >  1:3036910/1‑76 (MQ=255)
                                     ttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAttt                                   <  2:233297/76‑1 (MQ=255)
                                     ttatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTAttt                                   <  2:4417075/76‑1 (MQ=255)
                                      tatGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttg                                  >  1:2630796/1‑76 (MQ=255)
                                        tGGGGGTTCTATAATTCAGTGGTTTAATTTAAATCAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttgtt                                <  1:4791574/76‑1 (MQ=255)
                                        tGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttgtt                                >  1:1032198/1‑76 (MQ=255)
                                        tGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttgtt                                >  2:677458/1‑76 (MQ=255)
                                        tGGGGGTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttgtt                                <  2:6108112/76‑1 (MQ=255)
                                         gggggTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATttgttg                               <  1:4712731/76‑1 (MQ=255)
                                           gggTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAg                             >  2:5728713/1‑76 (MQ=255)
                                             gTTCTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGt                           >  1:230734/1‑76 (MQ=255)
                                                cTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCGGGAAATTATTTATTGGTTATTTGTTGAGGtttt                        <  2:2667591/76‑1 (MQ=255)
                                                cTATAATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGtttt                        >  2:6814321/1‑76 (MQ=255)
                                                 tataATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTc                       <  2:6384654/76‑1 (MQ=255)
                                                 tataATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTc                       >  1:3250635/1‑76 (MQ=255)
                                                   taATTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCtt                     >  2:5275783/1‑76 (MQ=255)
                                                    aaTTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTta                    >  2:4076793/1‑76 (MQ=255)
                                                     aTTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtat                   <  1:901052/76‑1 (MQ=255)
                                                     aTTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtat                   <  2:1689452/76‑1 (MQ=255)
                                                     aTTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtat                   >  1:2621946/1‑76 (MQ=255)
                                                     aTTCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtat                   <  1:3792148/76‑1 (MQ=255)
                                                      ttCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtata                  <  1:6322227/76‑1 (MQ=255)
                                                       tCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtatat                 <  1:3820404/76‑1 (MQ=255)
                                                       tCAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtatat                 >  1:4801443/1‑76 (MQ=255)
                                                        cAGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTtatata                >  2:4272820/1‑76 (MQ=255)
                                                         aGTGGTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTatatat               >  1:236486/1‑76 (MQ=255)
                                                            ggTTTAATTTAAAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTtg            >  1:7194803/1‑76 (MQ=255)
                                                              tttaatttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTtgtg          >  1:349209/1‑76 (MQ=255)
                                                               ttaatttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtgg         >  1:4354142/1‑76 (MQ=255)
                                                               ttaatttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtgg         >  1:154914/1‑76 (MQ=255)
                                                               ttaatttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtgg         >  2:3262291/1‑76 (MQ=255)
                                                                  atttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtggtgg      >  1:3403425/1‑76 (MQ=255)
                                                                   tttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtggtggt     >  2:7382526/1‑76 (MQ=255)
                                                                    ttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtggtggtg    >  1:5252898/1‑76 (MQ=255)
                                                                    ttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtggtggtg    <  1:6707113/76‑1 (MQ=255)
                                                                    ttaaAACAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTgtggtggtg    >  2:5259110/1‑76 (MQ=255)
                                                                      aaaaCAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTGTGGTGGTGtt  >  2:6278275/1‑76 (MQ=255)
                                                                      aaaaCAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTGTGGTGGTGtt  >  2:6153703/1‑76 (MQ=255)
                                                                         |                                                                        
ATAAACCTATTAGCTAATTTATATTAATGGCTGTTATTTATGGGGGTTCTATAATTCAGTGGTTTAATTTAAATCAACTAAAAATAACGCCGGAAATTATTTATTGGTTATTTGTTGAGGTTTTCTTATATATTTGTGGTGGTGTT  >  REL606/548462‑548607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: