Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,147,776	G→A	intergenic (+71/+8)	rpoD → / ← ygjF	RNA polymerase sigma factor/G/U mismatch‑specific DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,147,776	0	G	A	100.0%	191.5 / NA	50	intergenic (+71/+8)	rpoD/ygjF	RNA polymerase sigma factor/G/U mismatch‑specific DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/15);  total (35/15)

CGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTG  >  REL606/3147707‑3147848                                                                      |                                                                         cGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGttt                                                                    >  2:744697/1‑76 (MQ=255)   gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGttttt                                                                  >  2:1675456/1‑76 (MQ=255)    tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGtttttt                                                                 >  1:1344273/1‑76 (MQ=255)     aGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTa                                                                <  1:1298728/76‑1 (MQ=255)     aGGCCGGATCAGGCGTTACGCCGCAACCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTa                                                                >  2:494415/1‑76 (MQ=255)      ggCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTAt                                                               >  1:1194919/1‑76 (MQ=255)        ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCg                                                             <  1:142796/76‑1 (MQ=255)            aTCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCa                                                         >  2:45022/1‑76 (MQ=255)               aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGc                                                      >  2:1346908/1‑76 (MQ=255)               aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGc                                                      >  2:1628897/1‑76 (MQ=255)                ggCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCa                                                     >  2:1107321/1‑76 (MQ=255)                ggCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCa                                                     >  2:708812/1‑76 (MQ=255)                  cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACt                                                   <  2:70924/76‑1 (MQ=255)                   gTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTa                                                  <  1:1158721/76‑1 (MQ=255)                     tACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTAcc                                                >  1:956541/1‑76 (MQ=255)                      aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCa                                               <  1:1504145/76‑1 (MQ=255)                           gcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGcc                                          >  2:249956/1‑76 (MQ=255)                           gcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGcc                                          >  2:721667/1‑76 (MQ=255)                            cACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCt                                         >  2:1330893/1‑76 (MQ=255)                             aCCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                        >  2:157096/1‑76 (MQ=255)                             aCCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                        >  2:1209508/1‑76 (MQ=255)                               ccGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGt                                      >  2:1313468/1‑76 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCg                                                             >  1:1216724/1‑49 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCg                                                             <  2:1216724/49‑1 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                                  >  2:129176/1‑76 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                                  >  2:1271231/1‑76 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                                  <  1:748509/76‑1 (MQ=255)                                   cACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAg                                  <  1:1200815/76‑1 (MQ=255)                                      tAGGCCCTCTGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTc                               >  2:926497/1‑76 (MQ=255)                                       aGGCCCTCTGCACAAACGCCACNTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcg                              >  2:233856/1‑76 (MQ=255)                                            ctctGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATa                         <  1:1350428/76‑1 (MQ=255)                                              ctGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACg                       >  1:1661327/1‑76 (MQ=255)                                              ctGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACg                       >  2:297630/1‑76 (MQ=255)                                              ctGCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACg                       >  2:685417/1‑76 (MQ=255)                                                gCACAAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCt                     >  1:652629/1‑76 (MQ=255)                                                 cacaAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCtt                    >  2:614963/1‑76 (MQ=255)                                                  acaAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTc                   <  2:1122784/76‑1 (MQ=255)                                                   caAACGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCa                  <  1:926809/76‑1 (MQ=255)                                                    aaaCGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCaa                 >  1:803781/1‑76 (MQ=255)                                                     aaCGCCACCTTTTCGGTAGCGTNTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAc                >  2:228611/1‑76 (MQ=255)                                                      aCGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAAcc               <  2:1524250/76‑1 (MQ=255)                                                       cGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCa              >  2:1257797/1‑76 (MQ=255)                                                       cGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCa              >  2:634498/1‑76 (MQ=255)                                                       cGCCACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCa              >  2:118251/1‑76 (MQ=255)                                                         ccACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGt            >  1:640222/1‑76 (MQ=255)                                                          cACCTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGtt           >  2:1285623/1‑76 (MQ=255)                                                             cTTTTCGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTtct        <  2:804772/76‑1 (MQ=255)                                                                  cGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGt   >  1:281041/1‑76 (MQ=255)                                                                  cGGTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGt   <  2:145366/76‑1 (MQ=255)                                                                   ggTAGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTg  <  1:323940/76‑1 (MQ=255)                                                                   ggTAGCGTTTTGTATCGCCCACGCCCTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTg  <  1:323936/76‑1 (MQ=255)                                                                      |                                                                         CGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTG  >  REL606/3147707‑3147848

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