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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,260	0	T	A	100.0%	163.8 / NA	43	I15N (ATC→AAC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (16/27);  total (16/27)

CCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCA  >  REL606/3928192‑3928328                                                                     |                                                                     cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                               <  2:1150076/76‑1 (MQ=255) cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                               >  2:1123104/1‑76 (MQ=255)  ccAAGTTTAAGAACTCACACCACTGTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg                                                              <  2:1275750/76‑1 (MQ=255)    aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGTTTTCTGAGCTGAACACTCTCGGc                                                            <  2:89964/76‑1 (MQ=255)     aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                           >  1:713259/1‑76 (MQ=255)      gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGa                                                          <  1:1162083/76‑1 (MQ=255)      gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGa                                                          <  1:1162069/76‑1 (MQ=255)        ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTAGGCGaaa                                                        >  1:958791/1‑76 (MQ=255)          aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                      <  1:1236259/76‑1 (MQ=255)           aGAACTCACGCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                     <  1:963185/76‑1 (MQ=255)               cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                 <  1:423249/76‑1 (MQ=255)                tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                >  2:555693/1‑76 (MQ=255)                tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                >  2:886303/1‑76 (MQ=255)                    accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                            <  1:405214/76‑1 (MQ=255)                       acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                         <  1:650497/76‑1 (MQ=255)                          aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                      >  2:671671/1‑76 (MQ=255)                           tGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCt                                     <  1:244773/76‑1 (MQ=255)                            gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                                    >  2:1163976/1‑76 (MQ=255)                             aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                   <  2:531500/76‑1 (MQ=255)                                cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGCAAAACTGGCTc                                >  2:398840/1‑76 (MQ=255)                                cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                <  1:952455/76‑1 (MQ=255)                                cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                >  1:1232021/1‑76 (MQ=255)                                  tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                              >  2:165134/1‑76 (MQ=255)                                  tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                              <  1:971711/76‑1 (MQ=255)                                         ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                       <  2:1266720/76‑1 (MQ=255)                                         ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                       <  2:925200/76‑1 (MQ=255)                                         ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                       >  1:1180339/1‑76 (MQ=255)                                         ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                       <  1:180019/76‑1 (MQ=255)                                         ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                       >  1:801267/1‑76 (MQ=255)                                          tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                      <  1:474471/76‑1 (MQ=255)                                           aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                     <  2:1082693/76‑1 (MQ=255)                                           aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                     <  2:1189004/76‑1 (MQ=255)                                           aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                     >  2:1239646/1‑76 (MQ=255)                                            aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                    <  1:224924/76‑1 (MQ=255)                                              gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                  <  2:424373/76‑1 (MQ=255)                                               aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCggg                 <  2:1314604/76‑1 (MQ=255)                                               aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                 >  2:190364/1‑76 (MQ=255)                                               aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                 <  2:219372/76‑1 (MQ=255)                                                   cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat             <  2:1149331/76‑1 (MQ=255)                                                    gCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatt            >  2:891145/1‑76 (MQ=255)                                                      cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat          <  1:173573/76‑1 (MQ=255)                                                        gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt        <  1:282236/76‑1 (MQ=255)                                                              gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa  >  1:443015/1‑76 (MQ=255)                                                                     |                                                                     CCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCA  >  REL606/3928192‑3928328

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