Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,268,951	T→G	I129S (ATC→AGC)	atoC →	fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,268,951	0	T	G	100.0%	229.5 / NA	60	I129S (ATC→AGC)	atoC	fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (35/25);  total (35/25)

TAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACA  >  REL606/2268906‑2268998                                              |                                                tAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTc                                            <  1:6035263/51‑1 (MQ=255) tAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTc                                            <  1:175047/51‑1 (MQ=255) tAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTc                                            <  1:1294/51‑1 (MQ=255) tAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTc                                            <  1:2595714/51‑1 (MQ=255) tAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTc                                            <  1:5667755/51‑1 (MQ=255)  aaTCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCa                                           <  1:6581955/51‑1 (MQ=255)  aaTCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCa                                           <  1:478226/51‑1 (MQ=255)  aaTCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCa                                           >  1:5220150/1‑51 (MQ=255)  aaTCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCa                                           <  1:3289062/51‑1 (MQ=255)    tCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATc                                         <  1:1089639/51‑1 (MQ=255)    tCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATc                                         <  1:6988799/51‑1 (MQ=255)     cGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCt                                        <  1:5945129/51‑1 (MQ=255)     cGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCt                                        <  1:1521603/51‑1 (MQ=255)      gTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTg                                       >  1:2099986/1‑51 (MQ=255)      gTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTg                                       >  1:5785031/1‑51 (MQ=255)       ttCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGc                                      >  1:6087213/1‑51 (MQ=255)       ttCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGc                                      >  1:4008178/1‑51 (MQ=255)        tCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCa                                     >  1:4461191/1‑51 (MQ=255)        tCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCa                                     <  1:4073453/51‑1 (MQ=255)        tCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCa                                     <  1:6399602/51‑1 (MQ=255)         cAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCAc                                    >  1:3979151/1‑51 (MQ=255)          aGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCAcc                                   <  1:5412985/51‑1 (MQ=255)           gcgcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCa                                  <  1:6378189/51‑1 (MQ=255)            cgcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAg                                 <  1:1244452/51‑1 (MQ=255)            cgcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAg                                 >  1:3659574/1‑51 (MQ=255)             gcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAgg                                >  1:2865077/1‑51 (MQ=255)             gcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAgg                                >  1:1636164/1‑51 (MQ=255)             gcgcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAgg                                >  1:5601741/1‑51 (MQ=255)               gcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCa                              >  1:4035737/1‑51 (MQ=255)               gcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCa                              >  1:5265597/1‑51 (MQ=255)               gcTTTACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCa                              >  1:3399450/1‑51 (MQ=255)                 tttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACt                            >  1:1506062/1‑51 (MQ=255)                 tttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACt                            <  1:4231741/51‑1 (MQ=255)                 tttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACt                            >  1:3672650/1‑51 (MQ=255)                  ttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTg                           <  1:7276265/51‑1 (MQ=255)                  ttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTg                           <  1:4224723/51‑1 (MQ=255)                  ttACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTg                           >  1:65478/1‑51 (MQ=255)                   tACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGa                          <  1:1995031/51‑1 (MQ=255)                   tACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGa                          >  1:131221/1‑51 (MQ=255)                   tACAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGa                          <  1:5226810/51‑1 (MQ=255)                    aCAACTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAg                         <  1:2577701/51‑1 (MQ=255)                        cTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCAcc                     >  1:364555/1‑51 (MQ=255)                        cTCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCAcc                     >  1:4181563/1‑51 (MQ=255)                         tCCAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCa                    >  1:7134091/1‑51 (MQ=255)                          ccAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAg                   >  1:6468835/1‑51 (MQ=255)                          ccAGTCAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAg                   >  1:5467168/1‑51 (MQ=255)                               cAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGc              >  1:1607552/1‑51 (MQ=255)                               cAATGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGc              >  1:5939573/1‑51 (MQ=255)                                aaTGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCa             >  1:1568554/1‑51 (MQ=255)                                 aTGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCaa            <  1:3414303/51‑1 (MQ=255)                                 aTGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCa             >  1:6472319/1‑50 (MQ=255)                                  tGAAAAAAGAGAGCCGTCATCTGCACCCGGCACTGAGCACCAGCTGGCAAt           >  1:4298355/1‑51 (MQ=255)                                  tGAAAAAAGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAAt           >  1:4045934/1‑51 (MQ=255)                                    aaaaaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATgg         >  1:4303236/1‑51 (MQ=255)                                     aaaaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATggg        >  1:559969/1‑51 (MQ=255)                                      aaaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATgggg       >  1:4461163/1‑51 (MQ=255)                                      aaaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATgggg       >  1:4702143/1‑51 (MQ=255)                                       aaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATggggg      <  1:4898926/51‑1 (MQ=255)                                       aaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATggggg      >  1:1566791/1‑51 (MQ=255)                                       aaaGAGAGCCGTCATCTGCACCAGGCACTGAGCACCAGCGGGCAATggggg      >  1:311685/1‑51 (MQ=255)                                         agagagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGca    >  1:4915317/1‑51 (MQ=255)                                         agagagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGca    <  1:7127687/51‑1 (MQ=255)                                         agagagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGca    >  1:5710889/1‑51 (MQ=255)                                           agagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGcaca  >  1:526845/1‑51 (MQ=255)                                           agagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGcaca  >  1:7097758/1‑51 (MQ=255)                                           agagCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGcaca  >  1:1884037/1‑51 (MQ=255)                                              |                                                TAATCGTTCAGCGCGCTTTACAACTCCAGTCAATGAAAAAAGAGATCCGTCATCTGCACCAGGCACTGAGCACCAGCTGGCAATGGGGGCACA  >  REL606/2268906‑2268998

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