Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,928,259	A→T	I15F (ATC→TTC)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,259	0	A	T	100.0%	132.0 / NA	35	I15F (ATC→TTC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/18);  total (17/18)

ATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTT  >  REL606/3928186‑3928324                                                                          |                                                                  aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTTAGCTGTTc                                                                 <  2:639941/76‑1 (MQ=255) aTTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTc                                                                 <  2:541553/76‑1 (MQ=255)    tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct                                                              >  2:1021077/1‑76 (MQ=255)    tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct                                                              >  2:1003258/1‑76 (MQ=255)       cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACtctc                                                           <  2:797675/76‑1 (MQ=255)                aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAat                                                  <  1:529593/76‑1 (MQ=255)                   aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGATAATATg                                               >  1:1182041/1‑76 (MQ=255)                   aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATg                                               >  2:332011/1‑76 (MQ=255)                     cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATggg                                             >  2:1790/1‑76 (MQ=255)                     cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATggg                                             >  1:719931/1‑76 (MQ=255)                     cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATggg                                             >  1:719945/1‑76 (MQ=255)                       cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc                                           >  2:379644/1‑76 (MQ=255)                         caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTg                                         <  2:834383/76‑1 (MQ=255)                          accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTgg                                        >  1:1040924/1‑76 (MQ=255)                          accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTgg                                        >  1:892215/1‑76 (MQ=255)                          accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTgg                                        >  1:121524/1‑76 (MQ=255)                          accacTATGAATCTTACAGAATGAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGCAAATATGGGGCTgg                                        >  2:188330/1‑76 (MQ=255)                           ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGa                                       <  1:500196/76‑1 (MQ=255)                              cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGaaaa                                    >  2:292907/1‑76 (MQ=255)                                   aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                               >  2:316655/1‑76 (MQ=255)                                       ttACCGAATTAANGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                           >  1:16865/1‑76 (MQ=255)                                         aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                         >  1:770398/1‑76 (MQ=255)                                          ccGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                        <  1:599111/76‑1 (MQ=255)                                              aTTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                    <  2:740088/76‑1 (MQ=255)                                                 aaaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                 <  2:365922/76‑1 (MQ=255)                                                  aaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                <  2:935710/76‑1 (MQ=255)                                                     aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg             <  1:803139/76‑1 (MQ=255)                                                       tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc           <  1:1044901/76‑1 (MQ=255)                                                       tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc           <  1:405451/76‑1 (MQ=255)                                                        aCGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACa          <  1:946353/76‑1 (MQ=255)                                                        aCGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACa          <  1:946342/76‑1 (MQ=255)                                                           ccGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta       >  1:55909/1‑76 (MQ=255)                                                            cGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat      <  1:379644/76‑1 (MQ=255)                                                             ggTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt     <  1:396925/76‑1 (MQ=255)                                                               tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt   <  2:86852/76‑1 (MQ=255)                                                                ttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttttt  >  1:620334/1‑76 (MQ=255)                                                                          |                                                                  ATTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTT  >  REL606/3928186‑3928324

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