Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,520,689 (T)8→7 intergenic (‑101/‑45) dapA ← / → gcvR dihydrodipicolinate synthase/DNA‑binding transcriptional repressor, regulatory protein accessory to GcvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,520,6820T.98.0% 240.4 / ‑1.7 51intergenic (‑94/‑52)dapA/gcvRdihydrodipicolinate synthase/DNA‑binding transcriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (21/29);  total (21/30)
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.97e-01

ACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACA  >  REL606/2520607‑2520753
                                                                           |                                                                       
aCAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGt                                                                         >  1:1618903/1‑76 (MQ=255)
 cAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtt                                                                        <  1:652143/76‑1 (MQ=255)
  aaGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGttt                                                                       >  1:1403157/1‑76 (MQ=255)
   aGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                                                      >  1:1150658/1‑76 (MQ=255)
   aGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                                                      >  1:883053/1‑76 (MQ=255)
   aGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                                                      >  1:588475/1‑76 (MQ=255)
   aGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttt                                                                      >  1:1501650/1‑76 (MQ=255)
        ctcAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttat                                                                 >  1:434704/1‑74 (MQ=255)
         tcAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGtttttttatg                                                                <  2:1636978/76‑4 (MQ=255)
            aTGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCtt                                                            >  2:1197578/1‑76 (MQ=255)
                 aCGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCtt                                                       >  1:404706/1‑76 (MQ=255)
                  cGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTa                                                      <  2:691812/76‑1 (MQ=255)
                   gTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCATTa                                                      <  1:1252076/76‑1 (MQ=255)
                   gTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTaa                                                     >  1:258108/1‑76 (MQ=255)
                   gTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTaa                                                     <  1:926945/76‑1 (MQ=255)
                   gTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTaa                                                     >  2:118132/1‑76 (MQ=255)
                   gTTTGGAATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTaa                                                     >  2:1603586/1‑76 (MQ=255)
                    tttGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAg                                                    >  1:1319863/1‑76 (MQ=255)
                     ttGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGa                                                   <  2:215225/76‑1 (MQ=255)
                      tGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaa                                                  >  1:297913/1‑76 (MQ=255)
                      tGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaa                                                  <  1:729853/76‑1 (MQ=255)
                       ggTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaac                                                 >  2:588162/1‑76 (MQ=255)
                       ggTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaac                                                 >  1:847371/1‑76 (MQ=255)
                         tATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaacaa                                               <  1:1373003/76‑1 (MQ=255)
                         tATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaacaa                                               <  2:1633651/76‑1 (MQ=255)
                          aTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaacaac                                              <  1:301489/76‑1 (MQ=255)
                          aTGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGaacaac                                              <  1:1437901/76‑1 (MQ=255)
                           tGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACt                                             >  1:398347/1‑76 (MQ=255)
                           tGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACt                                             <  2:1436009/76‑1 (MQ=255)
                            ggCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTc                                            <  1:622003/76‑1 (MQ=255)
                            ggCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTc                                            >  1:1152479/1‑76 (MQ=255)
                            ggCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTc                                            <  1:622007/76‑1 (MQ=255)
                                  aaaaGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCt                                      <  2:1507779/76‑1 (MQ=255)
                                   aaaGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt                                     <  2:1162620/76‑1 (MQ=255)
                                   aaaGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt                                     <  1:1630206/76‑1 (MQ=255)
                                   aaaGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCtt                                     >  2:847152/1‑76 (MQ=255)
                                    aagcaagcAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTa                                    <  2:972414/76‑1 (MQ=255)
                                    aagcaagcAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTa                                    >  2:1010348/1‑76 (MQ=255)
                                      gcaagcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTaaa                                  <  1:240603/76‑1 (MQ=255)
                                      gcaagcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTaaa                                  <  2:772705/76‑1 (MQ=255)
                                       caagcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAg                                 <  2:1591487/76‑1 (MQ=255)
                                          gcaGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGaa                              <  1:1019269/76‑1 (MQ=255)
                                                agaACCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCa                        <  1:883528/76‑1 (MQ=255)
                                                  aaCCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGt                      <  1:603064/76‑1 (MQ=255)
                                                  aaCCGTTCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGt                      <  1:1166048/76‑1 (MQ=255)
                                                        tCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGaca                <  1:180548/76‑1 (MQ=255)
                                                        tCTGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGaca                <  1:180530/76‑1 (MQ=255)
                                                          tGATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACt              >  2:1451709/1‑76 (MQ=255)
                                                           gATTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTg             >  1:869406/1‑76 (MQ=255)
                                                            aTTGTTGTATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGt            >  2:1463980/1‑76 (MQ=255)
                                                              tgttgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCa          >  1:1305800/1‑76 (MQ=255)
                                                               gttgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCAt         >  2:170721/1‑76 (MQ=255)
                                                               gttgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCAt         >  2:879079/1‑76 (MQ=255)
                                                                 tgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCg       >  2:415967/1‑76 (MQ=255)
                                                                 tgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCg       <  2:713660/76‑1 (MQ=255)
                                                                 tgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCg       <  2:713663/76‑1 (MQ=255)
                                                                 tgtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCg       <  1:1147156/76‑1 (MQ=255)
                                                                  gtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGc      <  1:795723/76‑1 (MQ=255)
                                                                  gtATGCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGc      >  1:994571/1‑76 (MQ=255)
                                                                      gCATGTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACa  <  2:279412/76‑1 (MQ=255)
                                                                           |                                                                       
ACAAGTGTCTCAATGGTACGTTTGGTATGGCATTAAAAGCAAGCAGACAGAACCGTTCTGATTGTTGTATGCATGTTTTTTTTATGCTTTCCTTAAGAACAACTCACCCCTTAAAGGAATAACCAGTTTGACACTGTCATCGCAACA  >  REL606/2520607‑2520753

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