Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,308,106	G→A	A6A (GCG→GCA)	yhcB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,308,106	0	G	A	100.0%	275.8 / NA	81	A6A (GCG→GCA)	yhcB	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/41);  total (40/41)

CTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCG  >  REL606/3308058‑3308153                                                 |                                                cTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACt                                               <  1:479419/51‑1 (MQ=255)  tACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTa                                              <  1:6038378/51‑1 (MQ=255)  tACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTa                                              >  1:4144844/1‑51 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:4269558/51‑1 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:1179254/51‑1 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             >  1:6911426/1‑51 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             >  1:3397451/1‑51 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:4558825/51‑1 (MQ=255)   aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTaa                                             <  1:1949863/51‑1 (MQ=255)    cGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAt                                            <  1:2931817/51‑1 (MQ=255)     ggAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAtt                                           <  1:5059443/51‑1 (MQ=255)      gAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                          >  1:5374061/1‑51 (MQ=255)      gAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                          <  1:7192504/51‑1 (MQ=255)        aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                        <  1:4423004/51‑1 (MQ=255)        aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                        >  1:6062233/1‑51 (MQ=255)         aCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGt                                       >  1:5051921/1‑51 (MQ=255)          cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:4948040/1‑51 (MQ=255)          cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:6768410/1‑51 (MQ=255)          cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGtt                                      >  1:1958060/1‑51 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     <  1:4006964/51‑1 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     <  1:5448649/51‑1 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:889268/1‑51 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:2237172/1‑51 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:4825146/1‑51 (MQ=255)           aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTa                                     >  1:5439388/1‑51 (MQ=255)            aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    <  1:3268049/51‑1 (MQ=255)            aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    >  1:5662585/1‑51 (MQ=255)            aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAg                                    >  1:1412667/1‑51 (MQ=255)             aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   >  1:7236411/1‑51 (MQ=255)             aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   <  1:3916041/51‑1 (MQ=255)             aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                   <  1:1394189/51‑1 (MQ=255)              aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:4475438/1‑51 (MQ=255)              aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:327147/1‑51 (MQ=255)              aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  <  1:5826930/51‑1 (MQ=255)              aGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                  >  1:30401/1‑51 (MQ=255)                aCAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgt                                <  1:6802174/51‑1 (MQ=255)                 cAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgtc                               <  1:6574489/51‑1 (MQ=255)                 cAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgtc                               <  1:3870763/51‑1 (MQ=255)                  aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              <  1:5184232/51‑1 (MQ=255)                  aaCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGtcgtcg                              <  1:3258823/51‑1 (MQ=255)                   aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             >  1:3102172/1‑51 (MQ=255)                   aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             >  1:1920932/1‑51 (MQ=255)                   aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             <  1:4976322/51‑1 (MQ=255)                   aCGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCgg                             >  1:6138568/1‑51 (MQ=255)                     gggAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGca                           >  1:3950778/1‑51 (MQ=255)                      ggAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcat                          <  1:5894073/51‑1 (MQ=255)                        agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        <  1:6010/51‑1 (MQ=255)                        agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        >  1:6051567/1‑51 (MQ=255)                        agaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatca                        <  1:5449326/51‑1 (MQ=255)                         gaTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatcat                       <  1:590000/51‑1 (MQ=255)                          aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCatt                      <  1:4675863/51‑1 (MQ=255)                          aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCatt                      >  1:3896898/1‑51 (MQ=255)                            gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcatc                         <  1:5880414/46‑1 (MQ=255)                            gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCattat                    >  1:4931266/1‑51 (MQ=255)                               cATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTgg                 <  1:5517785/51‑1 (MQ=255)                                 tGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg               >  1:1515610/1‑51 (MQ=255)                                 tGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg               <  1:4779777/51‑1 (MQ=255)                                 tGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg               <  1:1693295/51‑1 (MQ=255)                                   aCCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt             <  1:5697433/51‑1 (MQ=255)                                   aCCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt             >  1:2059598/1‑51 (MQ=255)                                    ccTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtg            <  1:485044/51‑1 (MQ=255)                                     cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:1320321/1‑51 (MQ=255)                                     cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           >  1:3718431/1‑51 (MQ=255)                                     cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt           <  1:6063877/51‑1 (MQ=255)                                      tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:5976856/1‑51 (MQ=255)                                      tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:1529778/1‑51 (MQ=255)                                      tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          <  1:1101398/51‑1 (MQ=255)                                      tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg          >  1:722898/1‑51 (MQ=255)                                      tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtg            >  1:995331/1‑49 (MQ=255)                                       gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTgg         <  1:6939116/51‑1 (MQ=255)                                        ggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGc        <  1:5018023/51‑1 (MQ=255)                                         gAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGcc       >  1:6952292/1‑51 (MQ=255)                                          aaTATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCa      >  1:2871948/1‑51 (MQ=255)                                          aaTATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCa      <  1:2461494/51‑1 (MQ=255)                                           atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:3952858/1‑51 (MQ=255)                                           atatGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt     >  1:7188263/1‑51 (MQ=255)                                             atGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGc   >  1:3484280/1‑51 (MQ=255)                                             atGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGc   <  1:7064885/51‑1 (MQ=255)                                              tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  <  1:2998446/51‑1 (MQ=255)                                              tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  >  1:494198/1‑51 (MQ=255)                                              tGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCg  <  1:5462152/51‑1 (MQ=255)                                                 |                                                CTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCG  >  REL606/3308058‑3308153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑