Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,207,861 T→G N47T (AAC→ACC)  ycfC ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,207,8610TG100.0% 193.8 / NA 51N47T (AAC→ACC) ycfChypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (16/35);  total (16/35)

CGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCA  >  REL606/1207813‑1207908
                                                |                                               
cgcAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTc                                               <  1:1763645/51‑1 (MQ=255)
cgcAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTc                                               <  1:2057545/51‑1 (MQ=255)
  cAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCat                                             <  1:1282291/51‑1 (MQ=255)
  cAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCat                                             >  1:257440/1‑51 (MQ=255)
  cAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCat                                             <  1:1985690/51‑1 (MQ=255)
    ggTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCatat                                           <  1:3897400/51‑1 (MQ=255)
       tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaa                                        <  1:4266263/51‑1 (MQ=255)
        ggCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaat                                       >  1:814528/1‑51 (MQ=255)
          cTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     <  1:5449065/51‑1 (MQ=255)
          cTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     <  1:6651616/51‑1 (MQ=255)
          cTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     <  1:3780928/51‑1 (MQ=255)
          cTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     <  1:4544115/51‑1 (MQ=255)
          cTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     >  1:2268918/1‑51 (MQ=255)
              gCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACt                                 <  1:5813794/51‑1 (MQ=255)
               cTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTg                                <  1:1320862/51‑1 (MQ=255)
                tACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGt                               >  1:4368485/1‑51 (MQ=255)
                 aCCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGtt                              <  1:3335823/51‑1 (MQ=255)
                  ccGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGtt                              <  1:5250709/50‑1 (MQ=255)
                   cGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataa                                     <  1:5364581/42‑1 (MQ=255)
                    gCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgag                           >  1:21041/1‑51 (MQ=255)
                      cAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgagtg                         <  1:5224164/51‑1 (MQ=255)
                       aaaaaCCGCCAGCGTCGAGCTGGGGGTCatat                                           <  1:335092/32‑1 (MQ=255)
                       aaaaaCCGCCAGCGTCGAGCTGGGGGTCATATCaa                                        <  1:1008522/35‑1 (MQ=255)
                       aaaaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgagtga                        <  1:3443975/51‑1 (MQ=255)
                       aaaaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgagtga                        >  1:3104040/1‑51 (MQ=255)
                       aaaaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtg                             <  1:123817/46‑1 (MQ=255)
                         aaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTgaga                      >  1:7006797/1‑51 (MQ=255)
                         aaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTgaga                      <  1:6735412/51‑1 (MQ=255)
                         aaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTgaga                      <  1:5744597/51‑1 (MQ=255)
                          aaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGAc                     >  1:529642/1‑51 (MQ=255)
                           aCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACg                    <  1:3849179/51‑1 (MQ=255)
                           aCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACg                    <  1:4117452/51‑1 (MQ=255)
                            ccGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACgt                   >  1:551025/1‑51 (MQ=255)
                            ccGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACgt                   >  1:2871200/1‑51 (MQ=255)
                             cGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACgtg                  <  1:3104247/51‑1 (MQ=255)
                             cGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACgtg                  <  1:2557184/51‑1 (MQ=255)
                              gCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACgtgt                 <  1:4697481/51‑1 (MQ=255)
                               ccAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTa                <  1:1952719/51‑1 (MQ=255)
                               ccAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTa                <  1:1639693/51‑1 (MQ=255)
                                cAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAg               <  1:689174/51‑1 (MQ=255)
                                 aGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAgc              <  1:1811976/51‑1 (MQ=255)
                                 aGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAgc              >  1:4665917/1‑51 (MQ=255)
                                   cGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAgcgc            <  1:7007324/51‑1 (MQ=255)
                                   cGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAgcgc            >  1:1632732/1‑51 (MQ=255)
                                     tCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCAt          <  1:6946603/51‑1 (MQ=255)
                                       gAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCg        >  1:7179082/1‑51 (MQ=255)
                                        aGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCgg       <  1:5691805/51‑1 (MQ=255)
                                        aGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCgg       >  1:2597292/1‑51 (MQ=255)
                                        aGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCgg       >  1:2756202/1‑51 (MQ=255)
                                            ggggGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATc   <  1:5073642/51‑1 (MQ=255)
                                            ggggGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATc   >  1:391125/1‑51 (MQ=255)
                                             ggggTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATca  <  1:342759/51‑1 (MQ=255)
                                                |                                               
CGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCA  >  REL606/1207813‑1207908

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