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

TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCC  >  REL606/3928187‑3928331                                                                          |                                                                        tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa                                                                       >  2:1205918/1‑76 (MQ=255)  ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                      >  2:934928/1‑76 (MQ=255)   tACCCCAAGTTTAAGAACTCACAACACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                     <  2:171102/76‑1 (MQ=255)    aCCCCAAGTTTAAGAACTCCCACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc                                                                    >  2:1402869/1‑76 (MQ=255)      cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActctn                                                                  <  1:1401526/76‑2 (MQ=255)      cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                  <  1:579178/76‑1 (MQ=255)       ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg                                                                 <  2:1195743/76‑1 (MQ=255)            tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                            <  2:1310086/76‑1 (MQ=255)            tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa                                                            >  1:522441/1‑76 (MQ=255)             ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                           <  1:300570/76‑1 (MQ=255)              tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                          >  1:1020778/1‑76 (MQ=255)                aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                        >  2:911749/1‑76 (MQ=255)                    cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                    >  1:595354/1‑76 (MQ=255)                     tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg                                                   >  2:583350/1‑76 (MQ=255)                       acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt                                                 <  2:943768/76‑1 (MQ=255)                        caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                >  2:786333/1‑76 (MQ=255)                        caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                <  2:52948/76‑1 (MQ=255)                        caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                <  2:1451823/76‑1 (MQ=255)                         accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                               <  2:790051/76‑1 (MQ=255)                         accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                               <  1:829327/76‑1 (MQ=255)                             cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa                                           >  2:116380/1‑76 (MQ=255)                               aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                         <  1:1491986/76‑1 (MQ=255)                                  aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                      >  1:1036323/1‑76 (MQ=255)                                  aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                      >  1:1036333/1‑76 (MQ=255)                                  aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                      >  2:555133/1‑76 (MQ=255)                                     cTTACCGAATTAAAGAAGACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                   <  1:824/76‑1 (MQ=255)                                      ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                  >  1:1273034/1‑76 (MQ=255)                                         ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                               <  2:14306/76‑1 (MQ=255)                                         ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                               <  1:1292882/76‑1 (MQ=255)                                          cGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg                              <  1:789193/76‑1 (MQ=255)                                           gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATgn                             <  1:1327640/76‑2 (MQ=255)                                           gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                             <  1:1486049/76‑1 (MQ=255)                                              ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                          <  2:1171532/76‑1 (MQ=255)                                              ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                          <  2:69544/76‑1 (MQ=255)                                              ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                          <  2:583506/76‑1 (MQ=255)                                              ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                          <  2:1425948/76‑1 (MQ=255)                                               tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                         >  1:293041/1‑76 (MQ=255)                                               tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                         <  2:146275/76‑1 (MQ=255)                                               tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCAGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGGaa                         >  1:293049/1‑76 (MQ=255)                                               tAAAGAAAACGCCGGTTTCTTAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                         <  2:146259/76‑1 (MQ=255)                                                aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                        >  2:681367/1‑76 (MQ=255)                                                aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                        <  2:1486462/76‑1 (MQ=255)                                                aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                        >  2:67155/1‑76 (MQ=255)                                                 aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                       <  2:781864/76‑1 (MQ=255)                                                 aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                       <  1:861028/76‑1 (MQ=255)                                                  aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                      >  1:924715/1‑76 (MQ=255)                                                  aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                      <  2:333121/76‑1 (MQ=255)                                                   gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                     <  2:1123074/76‑1 (MQ=255)                                                    aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg                    <  2:1363004/76‑1 (MQ=255)                                                     aTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa                   <  2:1207718/76‑1 (MQ=255)                                                     aTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa                   >  1:617607/1‑76 (MQ=255)                                                        cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                <  1:849938/76‑1 (MQ=255)                                                         gCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatt               <  1:17404/76‑1 (MQ=255)                                                           cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat             <  1:87128/76‑1 (MQ=255)                                                           cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat             >  2:886786/1‑76 (MQ=255)                                                              tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt          <  1:791258/76‑1 (MQ=255)                                                              tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACACTCtttt          >  2:1403770/1‑76 (MQ=255)                                                                 cTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc       <  2:213507/76‑1 (MQ=255)                                                                    aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACATGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt    <  2:162699/76‑1 (MQ=255)                                                                     gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc   >  2:536991/1‑76 (MQ=255)                                                                     gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATc   <  1:1091586/76‑1 (MQ=255)                                                                      cTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc  >  1:695202/1‑76 (MQ=255)                                                                      cTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc  >  1:199464/1‑76 (MQ=255)                                                                          |                                                                        TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCC  >  REL606/3928187‑3928331

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