Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,928,901	G→T	A229S (GCT→TCT)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,901	0	G	T	100.0%	136.0 / NA	38	A229S (GCT→TCT)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (23/15);  total (23/15)

GATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCC  >  REL606/3928828‑3928973                                                                          |                                                                         gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCt                                                                        <  1:665988/76‑1 (MQ=255) gATGGTTCTGCTGATCAACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCt                                                                        >  2:923982/1‑76 (MQ=255)     gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTa                                                                    <  2:1000253/76‑1 (MQ=255)       tCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTAcc                                                                  >  1:994575/1‑76 (MQ=255)        ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCt                                                                 <  2:800661/76‑1 (MQ=255)        ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCt                                                                 >  1:333560/1‑76 (MQ=255)         tgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCtt                                                                >  2:987333/1‑76 (MQ=255)            tgATCGACNAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGTGGAAGTTGGTTCTTCTAACTTTGa                                                             >  1:1428537/1‑76 (MQ=255)            tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGCGAAGTTGTTTCTTCTACCTTTGa                                                             >  1:1096191/1‑76 (MQ=255)             gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGAc                                                            <  1:885819/76‑1 (MQ=255)             gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGAc                                                            <  2:321108/76‑1 (MQ=255)              aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACg                                                           >  1:1307692/1‑76 (MQ=255)                   cgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAAccc                                                      >  2:280379/1‑76 (MQ=255)                   cgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAAccc                                                      >  1:331830/1‑76 (MQ=255)                      aCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCa                                                   >  2:1242742/1‑76 (MQ=255)                         tCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCAtct                                                >  1:434368/1‑76 (MQ=255)                           cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCg                                              <  2:1100165/76‑1 (MQ=255)                           cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCg                                              >  2:1138188/1‑76 (MQ=255)                             gaagaaGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGcc                                            >  2:790968/1‑76 (MQ=255)                              aagaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCa                                           <  1:46154/76‑1 (MQ=255)                               agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCAc                                          <  1:555199/76‑1 (MQ=255)                                 aagTAACCGAGATGCAGNCTCTGGTAAAAGGGGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATATCGCCACGt                                        >  1:190721/1‑76 (MQ=255)                                     aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAg                                    >  1:1079627/1‑76 (MQ=255)                                        cGAGATGCAGCGTCTGGTAAAAGGAGAANGTGTTTCTTCTACCTCTGAAGNACCCGCATCTCGCCACGTTCAGGtt                                 >  1:211850/1‑76 (MQ=255)                                        cGAGANGGAGCGTCTGGTAAAAGGTGAACATGTTTCTTCTACCTTTGACGAACCCGCATCTNGCCACGTTCAGGtt                                 >  2:1484169/1‑76 (MQ=255)                                         gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTg                                <  2:1311071/76‑1 (MQ=255)                                            aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCgg                             >  1:739403/1‑76 (MQ=255)                                            aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCgg                             <  1:1020783/76‑1 (MQ=255)                                                   gTCTGGTAAAAGGTGAAGTTGTTTCTTNTANCTNTGGCGAANCCTCATGTCGCCACGTTCAGGTTGCGGAAATGGt                      <  1:1287634/76‑1 (MQ=255)                                                     cTGGTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGa                    >  1:999459/1‑76 (MQ=255)                                                       ggTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGATCCCGCATCTCGNCACGTTCAGGTTGCGGAAATGGTGATc                  <  1:1292244/76‑1 (MQ=255)                                                        gTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCg                 >  1:1050945/1‑76 (MQ=255)                                                        gTAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCg                 >  2:1316424/1‑76 (MQ=255)                                                         tAAAAGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCga                >  1:1412505/1‑76 (MQ=255)                                                         tAAAAAGTGAAGTTGTTTCTTCTACCTTTGACGAGCCCGCATCTCGCCACGTTCAGGTTGCGGTAATGGTGATCga                <  2:55911/76‑1 (MQ=255)                                                          aataGGTGAAGTTGTTTCTTCTACCTTTGACGAGCCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgag               >  1:1354361/4‑76 (MQ=255)                                                          aaaaGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgag               <  1:642325/76‑1 (MQ=255)                                                             aGGTGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAg            >  1:665479/1‑76 (MQ=255)                                                                tGAAGTTGTTTCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGcc         <  2:1394855/76‑1 (MQ=255)                                                                     ttgtttCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg    >  1:703302/1‑76 (MQ=255)                                                                     ttgtttCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg    <  1:1378164/76‑1 (MQ=255)                                                                       gtttCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGcc  <  2:177024/76‑1 (MQ=255)                                                                          |                                                                         GATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCC  >  REL606/3928828‑3928973

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: ‑