Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,140,926 T→G E107A (GAA→GCA)  ECB_02035 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,140,9260TG100.0% 159.2 / NA 39E107A (GAA→GCA) ECB_02035hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (12/27);  total (12/27)

ATCACAGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTTCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTAACAACTTGCACCCA  >  REL606/2140857‑2140996
                                                                     |                                                                      
aTCACAGCTTCTTCTATAACATTGGTCATTGTCGGCGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGgtg                                                                  <  2:326737/76‑1 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              >  2:1127977/1‑76 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              >  2:609207/1‑76 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              >  2:658936/1‑76 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              >  2:814792/1‑76 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              >  2:824146/1‑76 (MQ=255)
    caGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGtata                                                              <  2:505493/76‑1 (MQ=255)
      gCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACg                                                            <  1:617869/76‑1 (MQ=255)
      gCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACg                                                            <  1:617884/76‑1 (MQ=255)
        ttcttcTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGt                                                          <  2:1052700/76‑1 (MQ=255)
           ttctATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAgg                                                       <  1:986034/76‑1 (MQ=255)
              tataACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTGTACGGTAGGACa                                                    <  1:1201876/76‑1 (MQ=255)
               ataACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGGCaa                                                   >  1:588160/1‑76 (MQ=255)
               ataACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACaa                                                   >  1:423255/1‑76 (MQ=255)
                 aaCATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAata                                                 >  2:751266/1‑76 (MQ=255)
                 aaCATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAata                                                 <  1:40660/76‑1 (MQ=255)
                   cATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAAtata                                               <  1:977239/76‑1 (MQ=255)
                    aTTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGGCAATATaa                                              >  1:940427/1‑76 (MQ=255)
                           aTTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCg                                       <  2:280438/76‑1 (MQ=255)
                            ttGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGt                                      <  1:841500/76‑1 (MQ=255)
                            ttGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGt                                      <  1:1126799/76‑1 (MQ=255)
                             tGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTc                                     <  1:136588/76‑1 (MQ=255)
                             tGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTc                                     >  2:488492/1‑76 (MQ=255)
                              gTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCa                                    >  2:1340517/1‑76 (MQ=255)
                               tCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCAt                                   <  2:1289999/76‑1 (MQ=255)
                               tCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCAt                                   <  1:1461624/76‑1 (MQ=255)
                                 gACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATcc                                 >  2:951599/1‑76 (MQ=255)
                                 gACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATcc                                 <  2:629259/76‑1 (MQ=255)
                                   cGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCt                               <  2:817228/76‑1 (MQ=255)
                                     cTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCttt                             <  2:959943/76‑1 (MQ=255)
                                     cTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCttt                             <  2:36563/76‑1 (MQ=255)
                                     cTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCttt                             <  2:959924/76‑1 (MQ=255)
                                     cTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCttt                             <  2:270671/76‑1 (MQ=255)
                                      tGGCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTa                            <  2:915514/76‑1 (MQ=255)
                                       ggCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTAc                           >  1:1114796/1‑76 (MQ=255)
                                       ggCGCTTCTTCCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTAc                           <  1:217499/76‑1 (MQ=255)
                                            ttcttcCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGttct                      <  2:1470955/76‑1 (MQ=255)
                                              cttcCGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGttcttc                    <  1:617449/76‑1 (MQ=255)
                                                  cGGGGCTTCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTa                >  1:209255/1‑76 (MQ=255)
                                                         tCTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTAACAACtt         >  2:606312/1‑76 (MQ=255)
                                                          cTGACTTCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTAACAACTTg        <  2:12808/76‑1 (MQ=255)
                                                                tCGCTGCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTAACAACTTGCACcca  <  2:1497962/76‑1 (MQ=255)
                                                                     |                                                                      
ATCACAGCTTCTTCTATAACATTGGTCATTGTCGACGCTGGCGCTTCTTCCGGGGCTTCTGACTTCGCTTCTGGTGTATACGGTAGGACAATATAAATGTTCGTCATCCCTTTACGTTCTTCACTAACAACTTGCACCCA  >  REL606/2140857‑2140996

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