Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,202,603 A→G intergenic (‑44/‑156) iclR ← / → metH DNA‑binding transcriptional repressor/B12‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,202,6030AG98.3% 244.8 / ‑4.3 59intergenic (‑44/‑156)iclR/metHDNA‑binding transcriptional repressor/B12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (30/28);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGT  >  REL606/4202530‑4202674
                                                                         |                                                                       
gCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTCTTAtt                                                                       <  2:1508095/76‑1 (MQ=255)
 ccGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTttgttg                                                                      >  2:1592914/1‑76 (MQ=255)
    cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTa                                                                   >  1:487921/1‑76 (MQ=255)
     gTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAg                                                                  <  1:1187609/76‑1 (MQ=255)
      tttCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGc                                                                 <  2:325467/76‑1 (MQ=255)
       ttCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCt                                                                >  2:736466/1‑76 (MQ=255)
         cgcgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTaa                                                              >  1:590334/1‑76 (MQ=255)
          gcgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATCTTCATTTTTGTTGTTAGCTAAt                                                             <  2:924997/76‑1 (MQ=255)
           cgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATg                                                            <  2:748451/76‑1 (MQ=255)
           cgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATg                                                            <  2:23815/76‑1 (MQ=255)
             ggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCa                                                          >  1:606066/1‑76 (MQ=255)
             ggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCa                                                          >  1:1616288/1‑76 (MQ=255)
               aaTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAAt                                                        >  2:1362470/1‑76 (MQ=255)
                aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATa                                                       >  1:1027616/1‑76 (MQ=255)
                aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATa                                                       <  2:594336/76‑1 (MQ=255)
                 tGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAg                                                      >  2:1426966/1‑76 (MQ=255)
                      gCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCt                                                 <  2:683341/76‑1 (MQ=255)
                        gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGa                                               <  2:1325757/76‑1 (MQ=255)
                           cATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACt                                            >  2:452733/1‑76 (MQ=255)
                            aTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTg                                           >  1:1738491/1‑76 (MQ=255)
                              gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGAt                                         >  1:174718/1‑76 (MQ=255)
                              gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGAt                                         <  1:574158/76‑1 (MQ=255)
                              gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGAt                                         <  2:355226/76‑1 (MQ=255)
                               aCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATc                                        >  1:95826/1‑76 (MQ=255)
                                 aGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCg                                      >  1:742000/1‑76 (MQ=255)
                                   tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGAt                                    >  1:515827/1‑76 (MQ=255)
                                   tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGAt                                    >  1:1392891/1‑76 (MQ=255)
                                   tctcCTTTTTCTGTATCGTGGAAATCAATTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGAt                                    <  1:1479060/76‑1 (MQ=255)
                                    ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATg                                   <  1:1474212/76‑1 (MQ=255)
                                    ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATg                                   >  1:666924/1‑76 (MQ=255)
                                    ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATg                                   <  1:330446/76‑1 (MQ=255)
                                    ctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATg                                   <  1:80492/76‑1 (MQ=255)
                                     tcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGa                                  >  1:483414/1‑76 (MQ=255)
                                      ccTTTTTCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAg                                 <  1:1002987/76‑1 (MQ=255)
                                        tttttCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGtt                               <  2:1358426/76‑1 (MQ=255)
                                         ttttCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTa                              <  1:1737565/76‑1 (MQ=255)
                                           ttCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAAt                            <  1:1642346/76‑1 (MQ=255)
                                            tCTGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATg                           >  1:1496445/1‑76 (MQ=255)
                                              tGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGtt                         >  1:987363/1‑76 (MQ=255)
                                              tGTATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGtt                         >  2:1486468/1‑76 (MQ=255)
                                              tGNATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGtt                         >  2:1246634/1‑76 (MQ=255)
                                                tATCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGa                       <  1:1275609/76‑1 (MQ=255)
                                                 aTCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGaa                      >  1:771810/1‑76 (MQ=255)
                                                 aTCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGaa                      >  1:364081/1‑76 (MQ=255)
                                                  tCGTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAAc                     <  2:817843/76‑1 (MQ=255)
                                                    gTGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACaa                   <  1:126492/76‑1 (MQ=255)
                                                     tGGAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACaaa                  >  2:1370914/1‑76 (MQ=255)
                                                      ggAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAAt                 >  2:1270996/1‑76 (MQ=255)
                                                      ggAAATCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAAt                 >  1:1093332/1‑76 (MQ=255)
                                                        aaaTCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAAtct               >  1:779404/1‑76 (MQ=255)
                                                         aaTCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAAtctc              <  1:1592914/76‑1 (MQ=255)
                                                          aTCATTTTCATTTTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCa             <  1:1696184/76‑1 (MQ=255)
                                                           tcattttcatttTTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAt            <  1:1222291/76‑1 (MQ=255)
                                                             attttcattttTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGt          <  1:484304/76‑1 (MQ=255)
                                                             attttcattttTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGt          <  1:755714/76‑1 (MQ=255)
                                                              ttttcattttTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGtt         <  2:1616288/76‑1 (MQ=255)
                                                                ttcattttTGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGc       >  1:873153/1‑76 (MQ=255)
                                                                    tttttGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg   >  1:1196517/1‑76 (MQ=255)
                                                                     ttttGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggt  <  2:245026/76‑1 (MQ=255)
                                                                     ttttGTTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggt  >  1:958653/1‑76 (MQ=255)
                                                                         |                                                                       
GCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGT  >  REL606/4202530‑4202674

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