Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,411,844 G→T F130L (TTC→TTA ECB_04076 ← putative oxoacyl reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,411,8440GT93.5% 150.6 / 1.8 47F130L (TTC→TTAECB_04076putative oxoacyl reductase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (19/24);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

CTGACGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTGAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTTCGAG  >  REL606/4411776‑4411915
                                                                    |                                                                       
cTGACGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCCGTAAGCGCTTACTTAATGCCt                                                                  >  1:1532438/1‑76 (MQ=255)
 tGACGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTg                                                                 >  2:288375/1‑76 (MQ=255)
 tGACGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCCGCCGTGCCAGCAAACGCCTACTTAATGcccg                                                                 >  1:822321/1‑74 (MQ=255)
 tGACGAAAATAATCGAGCGGTGCGGCAGCGCCGCTTCCGGGTTTGGCTGGGGCAGATCGCGTTTACTTCATGccgg                                                                 >  2:229277/1‑74 (MQ=255)
   aCGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAg                                                               <  2:409858/76‑1 (MQ=255)
      aaaaTAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGGGCCAGTAAGCGCTTACTTAATGCCTGAGCaa                                                            >  2:1383304/1‑76 (MQ=255)
       aaaTAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAg                                                           <  1:716437/76‑1 (MQ=255)
       aaaTAAGCGAGCGGTGCGGCAGGTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAg                                                           <  1:700858/76‑1 (MQ=255)
         ataatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGga                                                         >  2:782509/1‑76 (MQ=255)
         ataatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGga                                                         <  2:147022/76‑1 (MQ=255)
          taatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaa                                                        <  2:1607825/76‑1 (MQ=255)
           aatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaag                                                       <  1:1745043/76‑1 (MQ=255)
           aatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaag                                                       <  2:963071/76‑1 (MQ=255)
           aatCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaag                                                       <  1:744417/76‑1 (MQ=255)
           aatCGAGAGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaag                                                       >  1:1170987/1‑76 (MQ=255)
            atCGAGCGGTGCGGCAGCTCCGCTTCCGGGTTCGGCTGTGCCAGTAAGCGATTACTTAATGCCTGAGCAAGgaaga                                                      >  2:155790/1‑76 (MQ=255)
             tCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaagaa                                                     <  1:1755968/76‑1 (MQ=255)
             tCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGgaagaa                                                     <  1:1184245/76‑1 (MQ=255)
              cGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAAc                                                    >  2:1092174/1‑76 (MQ=255)
              cGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTGAATGCCTGAGCAAGGAAGAAc                                                    <  1:273386/76‑1 (MQ=255)
               gAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGGCAGGAAGCGCCTACTTAATGCCTGAGGAAGggagaaag                                                   >  2:386010/1‑74 (MQ=255)
                 gCGGTGCGGCAGCTCCGCTTCCGGTGTCGGCTGTGCTAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGgc                                                 <  2:99633/76‑1 (MQ=255)
                  cGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAATAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGgcg                                                >  2:52825/1‑76 (MQ=255)
                   ggTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTTCTTAATGCCTGAGCCAGGAAGAACGgcgc                                               >  2:722411/1‑76 (MQ=255)
                   ggTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGgcgc                                               <  1:506873/76‑1 (MQ=255)
                    gTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGgcgcg                                              >  1:197515/1‑76 (MQ=255)
                       cGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGATGAACGGCGCGCgg                                           <  1:1201612/76‑1 (MQ=255)
                        ggCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCggg                                          >  1:749012/1‑76 (MQ=255)
                         gCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGgt                                         <  2:1388768/76‑1 (MQ=255)
                         gCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGgt                                         <  2:669982/76‑1 (MQ=255)
                           aGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTGAATGCCTGAGCAAGGAAGAACGGCGCGCGGgtgt                                       >  1:126987/1‑76 (MQ=255)
                             cTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTa                                     <  2:1051733/76‑1 (MQ=255)
                                cGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATg                                  <  2:1744272/76‑1 (MQ=255)
                                cGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATg                                  >  2:577409/1‑76 (MQ=255)
                                  cTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGc                                <  2:755966/76‑1 (MQ=255)
                                    tCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCaa                              <  1:585402/76‑1 (MQ=255)
                                         tttCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGttt                         >  1:480124/1‑76 (MQ=255)
                                         tttCGGCTGTGCCACTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGttt                         >  1:983756/1‑76 (MQ=255)
                                          ttCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGtttt                        >  1:425421/1‑76 (MQ=255)
                                           tCGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTg                       <  2:1626076/76‑1 (MQ=255)
                                            cGGCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGa                      <  1:940055/76‑1 (MQ=255)
                                              gCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATc                    <  1:408576/76‑1 (MQ=255)
                                              gCTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATc                    <  1:1424458/76‑1 (MQ=255)
                                               cTGTGCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCg                   >  1:467288/1‑76 (MQ=255)
                                                tgtgCCAGTAAGCGCTTACTGAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATTTTTTGATCGa                  <  2:180660/76‑1 (MQ=255)
                                                 gtgCCAGTAAGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGaa                 >  1:789600/1‑76 (MQ=255)
                                                         aaGCGCTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTg         >  2:219946/1‑76 (MQ=255)
                                                            cgcTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGtt      <  2:1065530/76‑1 (MQ=255)
                                                            cgcTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGtt      >  1:87350/1‑76 (MQ=255)
                                                            cgcTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGtt      >  1:1083217/1‑76 (MQ=255)
                                                             gcTTACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTtc     >  2:2146/1‑76 (MQ=255)
                                                              cttacttaATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTtcg    <  1:309573/76‑1 (MQ=255)
                                                               ttACTCAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTtcga   <  2:1812099/76‑1 (MQ=255)
                                                                tACTTAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTtcgag  >  1:905370/1‑76 (MQ=255)
                                                                    |                                                                       
CTGACGAAAATAATCGAGCGGTGCGGCAGCTCCGCTTCCGGTTTCGGCTGTGCCAGTAAGCGCTTACTGAATGCCTGAGCAAGGAAGAACGGCGCGCGGGTGTTAATGGCAATGTTTTGATCGAAGCTGTCTGGTTCGAG  >  REL606/4411776‑4411915

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