Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 283,483 T→C intergenic (‑526/‑1054) ykgM ← / → eaeH 50S ribosomal protein L31/attaching and effacing protein, pathogenesis factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606283,4830TC97.7% 878.5 / 4.7 213intergenic (‑526/‑1054)ykgM/eaeH50S ribosomal protein L31/attaching and effacing protein, pathogenesis factor
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base C (120/88);  total (124/89)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 217 total aligned reads displayed.

TATTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAA  >  REL606/283410‑283558
                                                                         |                                                                           
tATTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGctt                                                                           >  2:3926404/1‑76 (MQ=255)
tATTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGctt                                                                           >  2:28349/1‑76 (MQ=255)
 aTTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTg                                                                          >  1:4819606/1‑76 (MQ=255)
 aTTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTg                                                                          >  2:1231588/1‑76 (MQ=255)
 aTTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTg                                                                          >  2:5982224/1‑76 (MQ=255)
  ttCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTgg                                                                         >  1:379876/1‑76 (MQ=255)
  ttCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTgg                                                                         >  1:342120/1‑76 (MQ=255)
  ttCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTgg                                                                         <  1:3396956/76‑1 (MQ=255)
   tctcTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGt                                                                        <  2:131474/76‑1 (MQ=255)
    ctctTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTc                                                                       >  1:5387210/1‑76 (MQ=255)
     tctTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCa                                                                      <  2:317228/76‑1 (MQ=255)
     tctTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTAGGTCa                                                                      >  2:5871705/1‑76 (MQ=255)
      ctTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGTTTGGTCaa                                                                     >  2:367790/1‑76 (MQ=255)
       tttCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAAc                                                                    >  1:2554229/1‑76 (MQ=255)
       tttCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAAc                                                                    <  2:4483577/76‑1 (MQ=255)
        ttCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAAcc                                                                   <  1:3119785/76‑1 (MQ=255)
        ttCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAAcc                                                                   >  1:2670408/1‑76 (MQ=255)
         tCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCa                                                                  <  1:325963/76‑1 (MQ=255)
          cATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCaa                                                                 >  2:5361857/1‑76 (MQ=255)
          cATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCaa                                                                 >  1:3167486/1‑76 (MQ=255)
          cATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCaa                                                                 >  1:3620741/1‑76 (MQ=255)
             tGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAttt                                                              >  1:5185519/1‑76 (MQ=255)
             tGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAttt                                                              <  1:4746560/76‑1 (MQ=255)
             tGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAttt                                                              >  2:4302068/1‑76 (MQ=255)
              ggTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAtttt                                                             >  1:3868534/1‑76 (MQ=255)
               gTACTTCCTCGCTTTAGAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAttttt                                                            <  1:1257085/76‑1 (MQ=255)
               gTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAAttttt                                                            >  2:834731/1‑76 (MQ=255)
                 aCTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTg                                                          <  2:6286195/76‑1 (MQ=255)
                 aCTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTg                                                          >  1:1725238/1‑76 (MQ=255)
                  cTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGa                                                         <  1:1465989/76‑1 (MQ=255)
                  cTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGa                                                         <  2:930722/76‑1 (MQ=255)
                   ttCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGAc                                                        >  2:6254111/1‑76 (MQ=255)
                   ttCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGAc                                                        <  1:1216908/76‑1 (MQ=255)
                   ttCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGAc                                                        >  1:5822361/1‑76 (MQ=255)
                   ttCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGAc                                                        >  1:6002481/1‑76 (MQ=255)
                   ttCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGAc                                                        >  2:2239576/1‑76 (MQ=255)
                      cTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATg                                                     <  1:5740216/76‑1 (MQ=255)
                        cGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTa                                                   >  1:1111427/1‑76 (MQ=255)
                            ttAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATcac                                               >  2:1752127/1‑76 (MQ=255)
                            ttAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATcac                                               >  1:1156745/1‑76 (MQ=255)
                            ttAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATcac                                               >  1:1841292/1‑76 (MQ=255)
                            ttAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATcac                                               >  1:4515973/1‑76 (MQ=255)
                             tAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATcaca                                              >  2:4519658/1‑76 (MQ=255)
                              aaaaaaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACaa                                             >  2:6086034/1‑76 (MQ=255)
                              aaaaaaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACaa                                             <  1:2993308/76‑1 (MQ=255)
                               aaaaaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACaaa                                            >  2:4889704/1‑76 (MQ=255)
                                aaaaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAAt                                           >  2:1490419/1‑76 (MQ=255)
                                  aaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAAttt                                         >  1:2088286/1‑76 (MQ=255)
                                  aaGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAAttt                                         >  2:4736272/1‑76 (MQ=255)
                                  aaGAGTGCACATCGTAAGTGACCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAAttt                                         <  2:1439033/76‑1 (MQ=255)
                                    gagTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGa                                       <  1:4906940/76‑1 (MQ=255)
                                      gTGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAAt                                     <  1:3972992/76‑1 (MQ=255)
                                       tGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATa                                    >  1:6318700/1‑76 (MQ=255)
                                       tGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATa                                    >  2:4023571/1‑76 (MQ=255)
                                       tGCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATa                                    >  2:3210384/1‑76 (MQ=255)
                                        gCACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAg                                   <  2:1830149/76‑1 (MQ=255)
                                         cACTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGa                                  >  1:1802136/1‑76 (MQ=255)
                                          aCTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGAt                                 >  2:2662180/1‑76 (MQ=255)
                                           cTTCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATg                                <  2:277320/76‑1 (MQ=255)
                                            ttCGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGt                               <  1:974612/76‑1 (MQ=255)
                                              cGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAt                             >  1:2833676/1‑76 (MQ=255)
                                              cGTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAt                             >  2:4155117/1‑76 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATg                                <  2:5513508/72‑1 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATg                                >  1:5513508/1‑72 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAtt                            <  2:2080952/76‑1 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAtt                            >  1:4426630/1‑76 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAtt                            >  1:754524/1‑76 (MQ=255)
                                               gTAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTAtt                            <  1:3677427/76‑1 (MQ=255)
                                                tAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTa                           <  1:1276212/76‑1 (MQ=255)
                                                tAAGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTa                           >  2:5234150/1‑76 (MQ=255)
                                                 aaGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTAc                          >  1:1189388/1‑76 (MQ=255)
                                                 aaGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTAc                          >  1:4145542/1‑76 (MQ=255)
                                                  aGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACa                         <  1:5035401/76‑1 (MQ=255)
                                                  aGTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACa                         <  1:3713048/76‑1 (MQ=255)
                                                   gTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACAt                        <  2:2514762/76‑1 (MQ=255)
                                                   gTGCCCTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACAt                        >  1:3906952/1‑76 (MQ=255)
                                                    tGCCCTTATATAAATAACGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATc                       >  1:5883276/1‑76 (MQ=255)
                                                      cccTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCaa                     <  2:1376859/76‑1 (MQ=255)
                                                      cccTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCaa                     >  2:5116051/1‑76 (MQ=255)
                                                        cTTATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACt                   >  2:2162001/1‑76 (MQ=255)
                                                         ttATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTa                  <  2:3952569/76‑1 (MQ=255)
                                                         ttATATAAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTa                  >  1:4292524/1‑76 (MQ=255)
                                                           atataAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATc                >  1:4357436/1‑76 (MQ=255)
                                                            tataAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCt               >  1:5637038/1‑76 (MQ=255)
                                                            tataAATAACGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCt               <  2:3492689/76‑1 (MQ=255)
                                                             ataaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCtt              >  2:498789/1‑76 (MQ=255)
                                                             ataaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCtt              <  2:4592690/76‑1 (MQ=255)
                                                             ataaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCtt              >  1:3442528/1‑76 (MQ=255)
                                                              taaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCttt             >  1:2657196/1‑76 (MQ=255)
                                                               aaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCtttt            <  1:889412/76‑1 (MQ=255)
                                                               aaataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCtttt            <  2:3888203/76‑1 (MQ=255)
                                                                aataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTa           <  2:5937925/76‑1 (MQ=255)
                                                                 ataaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTAt          >  2:438594/1‑76 (MQ=255)
                                                                  taaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTAtt         <  2:1181409/76‑1 (MQ=255)
                                                                  taaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTAtt         >  1:3940333/1‑76 (MQ=255)
                                                                  taaCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTAtt         <  2:3546106/76‑1 (MQ=255)
                                                                    aCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGc       >  2:3380893/1‑76 (MQ=255)
                                                                    aCGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGc       <  2:2006468/76‑1 (MQ=255)
                                                                     cGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCa      <  2:1043902/76‑1 (MQ=255)
                                                                     cGAGCTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCa      >  1:6123526/1‑76 (MQ=255)
                                                                         |                                                                           
TATTCTCTTTCATTGGTACTTCCTCGCTTTAAAAAAGAGTGCACTTCGTAAGTGCCCTTATATAAATAACGAGTTTGGTCAACCAATTTTTTGACATGTATCACAAATTTGAATAGATGTATTACATCAACTATCTTTTATTGCACCAA  >  REL606/283410‑283558

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