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%	266.4 / NA	68	I15N (ATC→AAC)	rho	transcription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (26/42);  total (26/42)

TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928188‑3928332                                                                         |                                                                         ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       <  1:612183/76‑1 (MQ=255) ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       <  1:612167/76‑1 (MQ=255)  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      <  2:304602/76‑1 (MQ=255)  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      <  1:1286604/76‑1 (MQ=255)   aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc                                                                     >  1:792067/1‑76 (MQ=255)    ccccAAGTTTACGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct                                                                    >  1:34714/1‑76 (MQ=255)    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct                                                                    >  1:34698/1‑76 (MQ=255)        aggTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  2:802755/74‑1 (MQ=255)        aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  1:962425/76‑1 (MQ=255)         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               >  2:721435/1‑76 (MQ=255)           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                             >  1:326423/1‑76 (MQ=255)            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                            >  1:1194585/1‑76 (MQ=255)             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           <  1:1654877/76‑1 (MQ=255)              aaGAACTCACACCACTATTGATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          <  1:748217/76‑1 (MQ=255)              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  2:221065/1‑76 (MQ=255)               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         >  2:1723012/1‑76 (MQ=255)               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         >  2:175319/1‑76 (MQ=255)                gAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAatat                                                        >  2:383445/1‑76 (MQ=255)                 aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATg                                                       <  2:1488376/76‑1 (MQ=255)                  aCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgg                                                      <  1:1517778/76‑1 (MQ=255)                  aCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgg                                                      <  1:803556/76‑1 (MQ=255)                     cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGc                                                   <  1:325240/76‑1 (MQ=255)                      acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                  <  2:511222/76‑1 (MQ=255)                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 <  2:1082234/76‑1 (MQ=255)                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 <  1:1651893/76‑1 (MQ=255)                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                >  2:513435/1‑76 (MQ=255)                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                <  1:194660/76‑1 (MQ=255)                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                >  2:1008524/1‑76 (MQ=255)                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              <  1:789080/76‑1 (MQ=255)                           acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  1:687972/76‑1 (MQ=255)                            cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa                                            <  1:796619/76‑1 (MQ=255)                             tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc                                           >  1:972795/1‑76 (MQ=255)                                gAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTg                                        >  2:1293973/1‑76 (MQ=255)                                 aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                       <  2:515313/76‑1 (MQ=255)                                 aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                       <  2:469864/76‑1 (MQ=255)                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    <  2:751985/76‑1 (MQ=255)                                      tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                                  <  2:970091/76‑1 (MQ=255)                                      tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                                  >  1:985987/1‑76 (MQ=255)                                       aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACGCTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                                 <  2:97697/76‑1 (MQ=255)                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  1:192825/76‑1 (MQ=255)                                         cGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg                               <  1:838930/76‑1 (MQ=255)                                           aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg                             >  1:1663949/1‑76 (MQ=255)                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                            >  2:606130/1‑76 (MQ=255)                                            aTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGt                            <  1:541257/76‑1 (MQ=255)                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           <  1:1421138/76‑1 (MQ=255)                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           <  1:571063/76‑1 (MQ=255)                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  1:1618667/76‑1 (MQ=255)                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         >  1:1487715/1‑76 (MQ=255)                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        >  1:1154309/1‑76 (MQ=255)                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  2:1061672/76‑1 (MQ=255)                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       <  1:705017/76‑1 (MQ=255)                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       <  1:400680/76‑1 (MQ=255)                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  1:1005887/76‑1 (MQ=255)                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  2:1264328/76‑1 (MQ=255)                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  1:97653/76‑1 (MQ=255)                                                   aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                     <  1:721435/76‑1 (MQ=255)                                                     tACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc                   >  2:955583/1‑76 (MQ=255)                                                         ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               <  2:1114571/76‑1 (MQ=255)                                                           ggTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             <  1:691924/76‑1 (MQ=255)                                                           ggTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             >  1:115279/1‑76 (MQ=255)                                                           ggTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             <  1:923675/76‑1 (MQ=255)                                                            gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGNATGCGTAAGCAGGACATTAttt            <  2:235005/76‑1 (MQ=255)                                                             tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt           >  1:1132951/1‑76 (MQ=255)                                                             tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt           >  1:965362/1‑76 (MQ=255)                                                                cTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc        <  2:1132951/76‑1 (MQ=255)                                                                 tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc       <  1:577405/76‑1 (MQ=255)                                                                   aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt     <  1:949498/76‑1 (MQ=255)                                                                    gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc    >  2:1324790/1‑76 (MQ=255)                                                                     cTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc   >  1:444783/1‑76 (MQ=255)                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  1:1256887/1‑76 (MQ=255)                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  >  1:855621/1‑76 (MQ=255)                                                                         |                                                                         TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928188‑3928332

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