Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,303,925 C→A 15.9% intergenic (+98/‑86) metJ → / → yiiX DNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,303,9250CA15.9% 205.6 / 20.7 83intergenic (+98/‑86)metJ/yiiXDNA‑binding transcriptional repressor, S‑adenosylmethionine‑binding/predicted peptidoglycan peptidase
Reads supporting (aligned to +/- strand):  ref base C (37/30);  new base A (7/6);  total (47/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

GCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTTATA  >  Exported/2303857‑2303995
                                                                      |                                                                      
gCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                                                        >  1:4704022/1‑69 (MQ=255)
     aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGAt                                                                   <  1:5246844/69‑1 (MQ=255)
     aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGAt                                                                   >  1:3972239/1‑69 (MQ=255)
     aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGAt                                                                   <  1:157414/69‑1 (MQ=255)
     aTGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGAt                                                                   <  1:2440714/69‑1 (MQ=255)
     aTGCGCTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGa                                                                    >  1:3907292/1‑70 (MQ=2)
      tGCGCTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGAt                                                                   >  1:3803956/1‑70 (MQ=2)
      tGCGCTACGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGAt                                                                   >  1:3459869/1‑70 (MQ=2)
         gcTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCg                                                             >  1:3102765/1‑71 (MQ=255)
          cTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGc                                                                        <  1:3405734/59‑1 (MQ=255)
           tACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           >  1:5184782/1‑71 (MQ=255)
           tACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           >  1:4047786/1‑71 (MQ=255)
            aCGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           >  1:3891404/1‑70 (MQ=255)
             cGCTTATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAGGCGtt                                                           <  1:3029182/71‑1 (MQ=2)
               cTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATa                                                                  >  1:5691889/1‑60 (MQ=255)
                ttATCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAGGCGtt                                                           <  1:303607/68‑1 (MQ=2)
                  aTCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCAcgccg                                                    >  1:3580256/1‑71 (MQ=255)
                  aTCAGGCCTACAGGATATCTGGCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCAcgc                                                      >  1:3866421/1‑71 (MQ=2)
                     aGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGAt                                                                   <  1:3365622/53‑1 (MQ=255)
                     aGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                 <  1:5221789/71‑1 (MQ=255)
                      ggCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                <  1:3792219/71‑1 (MQ=255)
                       gCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                               >  1:3514639/1‑71 (MQ=255)
                       gCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                <  1:2469343/70‑1 (MQ=255)
                         cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGt                                                            >  1:426077/1‑56 (MQ=255)
                         cTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                             <  1:5416851/71‑1 (MQ=255)
                           aCGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGa                                                                    <  1:2987731/46‑1 (MQ=255)
                            cGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                                 >  1:3843765/1‑48 (MQ=255)
                            cGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                                 >  1:2546639/1‑48 (MQ=255)
                            cGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                           >  1:3740576/1‑70 (MQ=255)
                            cGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                           >  1:3871214/1‑70 (MQ=255)
                            cGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                           >  1:4573581/1‑70 (MQ=255)
                             gTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                                                         >  1:3675332/1‑55 (MQ=255)
                                atatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaac                                      <  1:4141506/71‑1 (MQ=255)
                                 tatTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaaca                                     <  1:856903/71‑1 (MQ=255)
                                  atTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaa                                    >  1:2532159/1‑71 (MQ=255)
                                  atTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaa                                    >  1:2453031/1‑71 (MQ=255)
                                     tggCAACTTATTAAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCAaaaaa                                      >  1:4308557/3‑67 (MQ=14)
                                     tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAaacaac                                   >  1:948933/1‑70 (MQ=14)
                                     tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAaacaac                                   >  1:3159430/1‑70 (MQ=14)
                                     tGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                                 >  1:1396753/1‑40 (MQ=255)
                                     tGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATaa                                                                 >  1:1536045/1‑40 (MQ=255)
                                     tGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAgg                                                               >  1:5499850/1‑42 (MQ=255)
                                     tGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAgg                                                               >  1:4843121/1‑42 (MQ=255)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa                                                                 <  1:1560332/39‑1 (MQ=9)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa                                                                 <  1:4310830/39‑1 (MQ=9)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa                                                                 >  1:252804/2‑40 (MQ=9)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa                                                                 <  1:2203729/39‑1 (MQ=9)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATaa                                                                 <  1:2282605/39‑1 (MQ=9)
                                      cgCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATa                                                                  >  1:4462287/2‑39 (MQ=9)
                                       gCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaac                                   >  1:127821/1‑69 (MQ=255)
                                       gCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGAt                                                                   >  1:4527265/1‑37 (MQ=9)
                                          aTTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAaacaa                                    >  1:4600911/1‑65 (MQ=255)
                                           tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGcc                              >  1:3177189/1‑69 (MQ=11)
                                           tttATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGt                                                            <  1:3047079/40‑1 (MQ=255)
                                           tttATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGt                                                            <  1:3885705/40‑1 (MQ=255)
                                           tttATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGt                                                            <  1:2223893/40‑1 (MQ=255)
                                           tttATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGt                                                            <  1:5144503/40‑1 (MQ=255)
                                            ttattaAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCa                                                  <  1:1336786/49‑1 (MQ=14)
                                            ttattaAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCa                                                  <  1:3061514/49‑1 (MQ=14)
                                            ttattaAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCa                                                  <  1:3138717/49‑1 (MQ=14)
                                            ttattaAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCACGCCGCa                                                  <  1:1191464/49‑1 (MQ=14)
                                            ttATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACg                            >  1:5221284/1‑71 (MQ=255)
                                             tATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                                                         >  1:1490483/1‑41 (MQ=38)
                                             tATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                                                         >  1:2501847/1‑41 (MQ=38)
                                             tATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCa                                                         >  1:2428901/1‑41 (MQ=38)
                                             tATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATaaa                                       >  1:4563262/1‑59 (MQ=255)
                                             tATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGt                           >  1:2382162/1‑71 (MQ=255)
                                              aTTAAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCa                                                         >  1:1024020/1‑40 (MQ=14)
                                              aTTAAAATTGCATGAACTTGTAGGACGGATAAGGCGTTCa                                                         >  1:5552732/1‑40 (MQ=14)
                                               ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt                          <  1:3852601/70‑1 (MQ=255)
                                               ttGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGtt                          <  1:1826962/70‑1 (MQ=255)
                                                 gAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           <  1:85072/35‑1 (MQ=40)
                                                 gAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           <  1:4438711/35‑1 (MQ=40)
                                                 gAATTTGCACGAACTTGTAGGCCGGATAAGGCGtt                                                           <  1:1782587/35‑1 (MQ=40)
                                                   aTTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTc                       >  1:4089437/1‑69 (MQ=255)
                                                          cGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCa                              >  1:2079034/1‑55 (MQ=255)
                                                          cGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGa                                 >  1:5269627/1‑52 (MQ=16)
                                                              cTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCAc                             <  1:3035817/52‑1 (MQ=38)
                                                              cTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCAc                             <  1:4693721/52‑1 (MQ=38)
                                                                 gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCttt       >  1:4754269/1‑71 (MQ=255)
                                                                 gTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCt         <  1:2558807/69‑1 (MQ=255)
                                                                  tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCtttt      <  1:4882940/71‑1 (MQ=255)
                                                                  tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCtttt      <  1:1175348/71‑1 (MQ=255)
                                                                  tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCttt       <  1:2698612/70‑1 (MQ=255)
                                                                   aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCTCtttt      >  1:1746273/1‑70 (MQ=255)
                                                                      ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTtata  <  1:1662872/71‑1 (MQ=255)
                                                                      ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTtata  <  1:5623803/71‑1 (MQ=255)
                                                                      ccGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTtata  <  1:4756588/71‑1 (MQ=255)
                                                                      |                                                                      
GCCTGATGCGCTACGCTTATCAGGCCTACGTCATATTGCAATTTATTGAATTTGCACGAACTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGAGCACGTTGTCTGCGACCCACCGCTTTTTATA  >  Exported/2303857‑2303995

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