Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,658,199 T→G L10* (TTA→TGA)  ydgB → short chain dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,658,1990TG100.0% 204.3 / NA 50L10* (TTA→TGA) ydgBshort chain dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/26);  total (24/26)

GTGGGATAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATTAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGAT  >  REL606/1658126‑1658271
                                                                         |                                                                        
gtggGATAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATatgaa                                                                        >  1:797781/1‑76 (MQ=255)
 tggGATAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATatgaat                                                                       >  2:928088/1‑76 (MQ=255)
  gggATAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATatgaatt                                                                      >  2:1156235/1‑76 (MQ=255)
     aTAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACt                                                                   >  1:326313/1‑76 (MQ=255)
     aTAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACt                                                                   <  2:612229/76‑1 (MQ=255)
      tAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTg                                                                  >  1:776777/1‑76 (MQ=255)
       aaCGTTCATCCCATCGTTGGATTGCNAAGGAGATAACGATGGGTAACNCCCCGCCCTTGCCAATATGAATTACTgg                                                                 <  1:183228/76‑1 (MQ=255)
          gTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCgg                                                              >  2:801586/1‑76 (MQ=255)
               tCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGgtc                                                         <  2:749781/76‑1 (MQ=255)
                 ccATCGTTGGATAGCAATCGAGATGAGGACAGGTGGAACACGGACCTTACTAATATGCATTAGTGGCGCAGctcgt                                                       >  1:1264929/1‑76 (MQ=255)
                  cATCGTTGNATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGgtggtc                                                      >  1:1220205/1‑76 (MQ=255)
                  cATCGTTGGATAGCAAAGGCGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGggcgtc                                                      >  1:1080149/1‑76 (MQ=255)
                    tCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGc                                                    <  2:419714/76‑1 (MQ=255)
                     cGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCa                                                   <  2:228060/76‑1 (MQ=255)
                     cGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCa                                                   <  2:1201971/76‑1 (MQ=255)
                      gTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCAt                                                  <  2:702427/76‑1 (MQ=255)
                       ttGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATCTGAATTACTGGCGGAGGTCGTCGCATc                                                 <  1:1125023/76‑1 (MQ=255)
                        tGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCg                                                >  1:88317/1‑76 (MQ=255)
                        tGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTAGTCGCATCg                                                >  1:627468/1‑76 (MQ=255)
                         ggATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCgg                                               <  1:857907/76‑1 (MQ=255)
                         ggATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCgg                                               >  1:982704/1‑76 (MQ=255)
                           aTAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcc                                             >  2:120775/1‑76 (MQ=255)
                              gCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGGCGCATCGGcctcc                                          >  2:913790/1‑75 (MQ=255)
                               cAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgc                                         >  2:601235/1‑76 (MQ=255)
                                aaaGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGTGGTCGTCGCATCGGcctcgcc                                        <  2:899027/76‑1 (MQ=255)
                                aaaGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgcc                                        >  1:477232/1‑76 (MQ=255)
                                   ggAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgccctc                                     <  2:692602/76‑1 (MQ=255)
                                   ggAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgccctc                                     <  1:587372/76‑1 (MQ=255)
                                   ggAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgccctc                                     <  2:269967/76‑1 (MQ=255)
                                    gagaTAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGcctcgccctcg                                    <  2:239067/76‑1 (MQ=255)
                                      gaTAACGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCa                                  <  1:638426/76‑1 (MQ=255)
                                       aTAACGATGGGTAANACCCAGCCCTTGCCAATATGAATTACTGGCGGGGGTCGTCGCATCGGCCTCGCCCNCCCAt                                 <  2:1266010/76‑1 (MQ=255)
                                          aCGATGGGTAANACNCAGCCCTTGCCTATATGAATTACTGGCGGAGGTCGTCGCATCCTCCTCGCCCTCGCATGGc                              <  2:195628/76‑1 (MQ=255)
                                          aCGAGGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGc                              <  1:682967/76‑1 (MQ=255)
                                           cGATGGGTCAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCa                             >  1:1172854/1‑76 (MQ=255)
                                           cGATGGGTAAAACCCAGCCCTTGCCAATATGAATTACCGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCa                             >  1:1062567/1‑76 (MQ=255)
                                             aTGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCTTCGCATCGGCCTCGCCCTCGCATGGCAtt                           <  2:540861/76‑1 (MQ=255)
                                              tGGGTAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCAttt                          >  2:824062/1‑76 (MQ=255)
                                                  tAAAACCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCCTGGCATTTCAtt                      >  1:1005419/1‑76 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTa                     <  2:913735/76‑1 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTa                     <  2:913750/76‑1 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTa                     >  2:1076787/1‑76 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTa                     <  1:530720/76‑1 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTa                     <  1:786142/76‑1 (MQ=255)
                                                   aaaaCCCAGCCCTTGCCAATATGAATTACTGGCGCAGGTCGTCGCCGCGGCCACGCCCTCGCATTGCGTGTCAtca                     >  2:1076781/1‑74 (MQ=255)
                                                       cccAGCCCTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCa                 <  1:210815/76‑1 (MQ=255)
                                                        ccAGCCCTTGCCAATATGAATTACTGGCGGAGGTAGTCGCANCGGCCTCGCCCTCGCCCGNGANTTCATTAATCaa                >  1:1032253/1‑76 (MQ=255)
                                                             ccTTGGCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCa           >  1:81336/1‑76 (MQ=255)
                                                             ccTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCa           >  1:213816/1‑76 (MQ=255)
                                                              cTTGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCaa          >  2:1177890/1‑76 (MQ=255)
                                                                tGCCAATATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAAcc        <  2:213968/76‑1 (MQ=255)
                                                                    aaTATGAATTACTGGCGGAGTTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTTATCAAAAGCAACCGGTg    <  1:1253188/76‑1 (MQ=255)
                                                                    aaTATGAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTg    >  1:1025936/1‑76 (MQ=255)
                                                                      tatGAATTACTGGCGGAGGTCGTCGCATCGGCCTAGCNCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGAt  >  1:1232336/1‑76 (MQ=255)
                                                                         |                                                                        
GTGGGATAACGTTCATCCCATCGTTGGATAGCAAAGGAGATAACGATGGGTAAAACCCAGCCCTTGCCAATATTAATTACTGGCGGAGGTCGTCGCATCGGCCTCGCCCTCGCATGGCATTTCATTAATCAAAAGCAACCGGTGAT  >  REL606/1658126‑1658271

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