Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,928,515	A→C	K100T (AAG→ACG)	rho →	transcription termination factor Rho

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,928,515	0	A	C	100.0%	91.0 / NA	25	K100T (AAG→ACG)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/9);  total (16/9)

CTGATGACATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAC  >  REL606/3928443‑3928586                                                                         |                                                                        cTGATGACATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGAt                                                                      >  2:1029067/1‑76 (MQ=255)  tgatgaCATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGAtt                                                                     <  1:1189095/76‑1 (MQ=255)    atgaCATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTcg                                                                   >  1:529560/1‑76 (MQ=255)      gaCATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTcgcc                                                                 <  2:782011/76‑1 (MQ=255)        cATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTcgcccg                                                               <  1:126962/76‑1 (MQ=255)          tCTACGTTTCCCTTAGCCANANCCGCCGTNTCAACCTCCGCACTGGTGATANCANCTCTGGTACGATTcgcccgcc                                                             <  2:257017/76‑1 (MQ=255)           cTACGTTCCCCGTACCCAACCCCGCCTTTTCACCCTCCGCCCTGGTGATACCATCTCTGGTACGATTCGCccgccg                                                            <  2:720782/76‑1 (MQ=255)            tACTTTCCCCTTACCCAAACCCGCCGTTTCACCCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGa                                                           <  2:1145361/76‑1 (MQ=255)              cGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGaaa                                                         >  1:525166/1‑76 (MQ=255)                tttCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAaaga                                                       >  2:819257/1‑76 (MQ=255)                ttccccTTACCCAAACCCGCCGTTTCACCCCCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAaaga                                                       <  2:1182059/73‑1 (MQ=255)                  tCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAaagaag                                                     >  1:1362799/1‑76 (MQ=255)                  ccccTTCCCCAACCCCGCCTTTTCACCCCCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAaagaag                                                     <  1:1258767/75‑1 (MQ=255)                     cTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACNATCTCTGGTACGATTCGCNCGNCGAAAGAAGGTg                                                  >  1:1440037/1‑76 (MQ=255)                     cTTACCCAACCCCGCCGTTCCACCCCCCGCGCTGGTAATCCCATCTCGGGTAGGATTCGCCCGCCGAAAGAAGGTg                                                  <  2:877756/76‑1 (MQ=255)                        aGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAAc                                               >  1:1679040/1‑76 (MQ=255)                              aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTAtt                                         >  2:1087273/1‑76 (MQ=255)                                        tCACCCCCCGCACGGGTAATCCCATCTCTGGTACGTTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                               <  2:84647/76‑1 (MQ=255)                                         cacccTCCGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                              <  2:1494404/73‑1 (MQ=255)                                         cAACCTCCGCACCGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                              >  2:462676/1‑76 (MQ=255)                                            cccccGCACTGGTAATACCATCCCGGTTACGATCCCCCCGCCGAAAGAAGGTGAACGCTATTGTGCGCTGCTGaaa                           <  1:955007/73‑1 (MQ=255)                                               ccGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                        >  2:195132/1‑76 (MQ=255)                                                cGCACTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                       >  1:1381264/1‑76 (MQ=255)                                                cGCACGGGTGATCCCCTCTCTGGTACGATCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTa                       <  1:1029080/76‑1 (MQ=255)                                                    cTGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                   >  2:807037/1‑76 (MQ=255)                                                     tGGTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGaa                  <  2:1647086/76‑1 (MQ=255)                                                      ggTGATACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                 >  2:969880/1‑76 (MQ=255)                                                        tGATCCCATCTCTGGTACGTTCCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt               <  1:884534/76‑1 (MQ=255)                                                        tGATCCCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt               <  2:148767/76‑1 (MQ=255)                                                           tACCATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc            >  1:629938/1‑76 (MQ=255)                                                             ccctctcGGGTACGTTCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt          <  1:942042/73‑1 (MQ=255)                                                             ccATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt          <  1:1253065/76‑1 (MQ=255)                                                             ccATCTCGGGTGCGAGTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCAGCTGAAAGTNAACGAAGTTAACtt          <  1:613/76‑1 (MQ=255)                                                              cATCTCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc         >  2:314178/1‑76 (MQ=255)                                                              cATCTCTGGTACGATCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc         <  2:834073/76‑1 (MQ=255)                                                               aTCGCTGGTACGATTCGCCCGCCGAAAGAAGGTGAACGCTAGGTTGCGCTGCTGNAAGTTAACTANGTTTACTTCg        <  2:1734430/76‑1 (MQ=255)                                                                tctcGGGTACGATCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGa       <  2:757731/76‑1 (MQ=255)                                                                  tctGGTGCGATCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACa     <  2:1226040/76‑1 (MQ=255)                                                                   tgggTAAGATCCCCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa    <  1:1373427/74‑1 (MQ=255)                                                                     ggTACGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc  >  2:1053695/1‑76 (MQ=255)                                                                         |                                                                        CTGATGACATCTACGTTTCCCTTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAC  >  REL606/3928443‑3928586

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