Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,108,409 T→C K32E (AAG→GAG)  metJ ← transcriptional repressor protein MetJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,108,4090TC95.7% 274.8 / 0.6 69K32E (AAG→GAG) metJtranscriptional repressor protein MetJ
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (31/35);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATA  >  REL606/4108336‑4108479
                                                                         |                                                                      
tCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCaa                                                                      <  1:1695581/76‑1 (MQ=255)
tCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCaa                                                                      <  1:1390642/76‑1 (MQ=255)
 cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag                                                                     >  2:1458643/1‑76 (MQ=255)
 cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag                                                                     <  2:1819120/76‑1 (MQ=255)
 cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag                                                                     >  2:1740359/1‑76 (MQ=255)
 cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTATCGGTGAGGATTTTTAACACCTCAaga                                                                    <  1:724212/76‑1 (MQ=255)
    tGTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAgg                                                                  >  1:495482/1‑76 (MQ=255)
    ngTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAgg                                                                  >  1:110525/2‑76 (MQ=255)
     gTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGa                                                                 <  1:1067199/76‑1 (MQ=255)
      ttGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGaa                                                                >  2:1572340/1‑76 (MQ=255)
       tGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAAt                                                               <  1:505196/76‑1 (MQ=255)
         gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg                                                             <  1:842911/76‑1 (MQ=255)
         gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg                                                             <  1:1212055/76‑1 (MQ=255)
         gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg                                                             <  2:213376/76‑1 (MQ=255)
         gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg                                                             <  1:1484735/76‑1 (MQ=255)
          tAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGa                                                            >  1:1019985/1‑76 (MQ=255)
           aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGaa                                                           >  1:1215601/1‑76 (MQ=255)
           aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGaa                                                           >  2:308063/1‑76 (MQ=255)
             cGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAAc                                                         >  2:1455336/1‑76 (MQ=255)
              gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAAcc                                                        <  2:77844/76‑1 (MQ=255)
                 aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTa                                                     <  2:1095718/76‑1 (MQ=255)
                  cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTaa                                                    <  1:659084/76‑1 (MQ=255)
                  cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTaa                                                    <  1:18462/76‑1 (MQ=255)
                   ncAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAAt                                                   >  1:896390/2‑76 (MQ=255)
                     aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCt                                                 <  2:359982/76‑1 (MQ=255)
                     aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCt                                                 >  2:323533/1‑76 (MQ=255)
                       gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt                                               >  1:1049242/1‑76 (MQ=255)
                       gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt                                               >  1:544903/1‑76 (MQ=255)
                       gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt                                               <  2:747678/76‑1 (MQ=255)
                       gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt                                               >  1:799798/1‑76 (MQ=255)
                         tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt                                             <  1:1134293/76‑1 (MQ=255)
                         tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt                                             <  1:898696/76‑1 (MQ=255)
                         tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt                                             >  1:1407709/1‑76 (MQ=255)
                          gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTg                                            <  2:687399/76‑1 (MQ=255)
                          gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTg                                            >  1:291112/1‑76 (MQ=255)
                           ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGa                                           <  1:1454262/76‑1 (MQ=255)
                           ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGa                                           <  1:516797/76‑1 (MQ=255)
                               ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTg                                       <  2:1407709/76‑1 (MQ=255)
                                 tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGtt                                     >  2:564790/1‑76 (MQ=255)
                                 tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGtt                                     >  1:1245366/1‑76 (MQ=255)
                                  gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc                                    >  2:578843/1‑76 (MQ=255)
                                    cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc                                  <  1:1521117/76‑1 (MQ=255)
                                      acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc                                >  1:842197/1‑76 (MQ=255)
                                        gcgtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtctt                              <  2:486953/76‑1 (MQ=255)
                                          gtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                            <  2:873698/76‑1 (MQ=255)
                                          gtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct                            >  1:55595/1‑76 (MQ=255)
                                            gcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTg                          <  2:1115910/76‑1 (MQ=255)
                                             cgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                         >  2:350500/1‑76 (MQ=255)
                                             cgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc                         <  1:242853/76‑1 (MQ=255)
                                              gtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc                        <  2:1111795/76‑1 (MQ=255)
                                              gtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc                        >  1:1471748/1‑76 (MQ=255)
                                               ttCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg                       <  1:485690/76‑1 (MQ=255)
                                                    cGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc                  <  1:1045355/76‑1 (MQ=255)
                                                       tGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGc               <  2:48807/76‑1 (MQ=255)
                                                        gAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg              >  1:1775946/1‑76 (MQ=255)
                                                        gAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg              >  1:1775940/1‑76 (MQ=255)
                                                         aGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGt             <  2:1562840/76‑1 (MQ=255)
                                                             tttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATgg         <  1:6670/76‑1 (MQ=255)
                                                             tttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATgg         <  1:224196/76‑1 (MQ=255)
                                                              ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg        >  1:1540483/1‑76 (MQ=255)
                                                              ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg        >  1:883999/1‑76 (MQ=255)
                                                              ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg        >  1:1434591/1‑76 (MQ=255)
                                                               tttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGc       >  1:824396/1‑76 (MQ=255)
                                                                 tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg     <  2:470171/76‑1 (MQ=255)
                                                                 tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg     <  2:671992/76‑1 (MQ=255)
                                                                 tAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg     >  2:488127/1‑76 (MQ=255)
                                                                  aaCACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa    >  2:1881324/1‑76 (MQ=255)
                                                                   acacCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGat   <  1:1261106/76‑1 (MQ=255)
                                                                    cacCTCAAGAGGAATGGAAACCGTAATCGTTGTGACTTGTTCACTCTGCTTGCCGTGCTCAGCGTATGGGCTGata  <  2:329375/76‑1 (MQ=255)
                                                                         |                                                                      
TCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATA  >  REL606/4108336‑4108479

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: