Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,928,383	A→C	E56A (GAG→GCG)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,383	0	A	C	100.0%	180.7 / NA	43	E56A (GAG→GCG)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/24);  total (19/24)

CAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTA  >  REL606/3928310‑3928455                                                                          |                                                                         cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAGGATATCTTTGGTGATGGCGTACTGgcga                                                                        >  2:1294210/1‑76 (MQ=255)  aGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGat                                                                       >  1:1259508/1‑76 (MQ=255)   ggACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGata                                                                      <  1:802537/76‑1 (MQ=255)    gACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGatat                                                                     <  1:188240/76‑1 (MQ=255)     aCATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATAtt                                                                    <  2:1264964/76‑1 (MQ=255)               tGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAggatgga                                                          >  1:157751/1‑76 (MQ=255)               tGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAggatgga                                                          >  2:513725/1‑76 (MQ=255)                   aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTg                                                      <  2:140055/76‑1 (MQ=255)                   aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTg                                                      <  1:1260880/76‑1 (MQ=255)                   aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTg                                                      <  1:936275/76‑1 (MQ=255)                      cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGtt                                                   <  2:205765/76‑1 (MQ=255)                      cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGtt                                                   <  2:401116/76‑1 (MQ=255)                      cTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGtt                                                   >  1:461614/1‑76 (MQ=255)                        gAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtc                                                 <  2:872499/76‑1 (MQ=255)                         aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcc                                                >  2:399639/1‑76 (MQ=255)                         aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcc                                                >  2:11997/1‑76 (MQ=255)                         aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcc                                                >  2:437025/1‑76 (MQ=255)                          agcagcACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcct                                               >  2:1287973/1‑76 (MQ=255)                            cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcctcc                                             <  2:464325/76‑1 (MQ=255)                            cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTtcctcc                                             <  2:575336/76‑1 (MQ=255)                                   cAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGc                                      <  2:539697/76‑1 (MQ=255)                                     aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAg                                    <  1:976285/76‑1 (MQ=255)                                     aaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAg                                    >  2:928812/1‑76 (MQ=255)                                       gagTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGAc                                  <  2:713400/76‑1 (MQ=255)                                        agTGGCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACa                                 >  1:944512/1‑76 (MQ=255)                                           ggCGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCt                              >  1:487096/1‑76 (MQ=255)                                             cGAAGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTcc                            <  1:640034/76‑1 (MQ=255)                                               aaGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTa                          <  2:1532663/76‑1 (MQ=255)                                                aGATATCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAc                         >  1:1091925/1‑76 (MQ=255)                                                  atatCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCt                       >  2:133455/1‑76 (MQ=255)                                                   tatCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                      <  2:1233931/76‑1 (MQ=255)                                                   tatCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                      <  1:1552985/76‑1 (MQ=255)                                                    atCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCg                     >  1:1526110/1‑76 (MQ=255)                                                     tCTTTGGTGATGGCGTACTGGCGATATTGCAGGATGGATTTGTTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGc                    <  2:483779/76‑1 (MQ=255)                                                           gTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc              <  2:1051350/76‑1 (MQ=255)                                                           gTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc              >  1:574722/1‑76 (MQ=255)                                                           gTGATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc              >  1:1285649/1‑76 (MQ=255)                                                             gATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtg            <  2:325816/76‑1 (MQ=255)                                                             gATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtg            <  2:187667/76‑1 (MQ=255)                                                             gATGGCGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtg            <  2:1166548/76‑1 (MQ=255)                                                                  cGTACTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGAc       >  1:140169/1‑76 (MQ=255)                                                                     aCTGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATc    <  1:159645/76‑1 (MQ=255)                                                                       tGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTa  >  2:443383/1‑76 (MQ=255)                                                                       tGGCGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTa  >  1:1130571/1‑76 (MQ=255)                                                                          |                                                                         CAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTA  >  REL606/3928310‑3928455

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