Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,927,357	(T)8→9	coding (74/1266 nt)	rhlB ←	ATP‑dependent RNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,927,349	1	.	T	100.0%	161.8 / NA	40	G28R (GGG→AGG)	rhlB	ATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (25/15);  total (25/15)

CTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATG  >  REL606/3927275‑3927418                                                                            |                                                                       cTGCCCGGCTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑t                                                                      <  2:367421/76‑1 (MQ=255) cTGCCCGGCTACGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑t                                                                      <  2:590541/76‑1 (MQ=255)    cccGGCTCCGTCACGCCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑tttt                                                                   <  2:925008/76‑1 (MQ=255)     ccGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttt                                                                  <  1:644947/76‑1 (MQ=255)     ccGCCTCCGGCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttt                                                                  <  2:113115/76‑1 (MQ=255)         cgacgccacgccCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑ttttttttt                                                              <  2:1338687/74‑2 (MQ=255)          tACGCCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑tttttttttc                                                             <  1:319292/76‑3 (MQ=255)              tCACGGCCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAgg                                                          <  1:1327849/76‑1 (MQ=255)              tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAgg                                                          >  2:1243834/1‑76 (MQ=255)              tCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAgg                                                          >  2:807435/1‑76 (MQ=255)                aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGc                                                        >  1:881651/1‑76 (MQ=255)                aCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGc                                                        >  1:1144320/1‑76 (MQ=255)                  gACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCTTTTTTTTTTCAAGGGCt                                                       >  2:622403/1‑76 (MQ=255)                  cccccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCtt                                                      <  1:1163201/74‑1 (MQ=255)                    ccccccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCt                                                    <  2:24617/76‑1 (MQ=255)                     ccgccgGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTa                                                   <  1:331810/76‑1 (MQ=255)                     ccGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTa                                                   <  1:1013413/76‑1 (MQ=255)                       gCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACa                                                 >  2:367116/1‑76 (MQ=255)                        ccAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACaa                                                <  1:160393/76‑1 (MQ=255)                           gCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCt                                             >  2:204897/1‑76 (MQ=255)                                  cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                      >  2:1365978/1‑76 (MQ=255)                                  cGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATg                                      >  1:1479596/1‑76 (MQ=255)                                   ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                     >  2:762026/1‑76 (MQ=255)                                   ggAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGc                                     >  1:113675/1‑76 (MQ=255)                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  <  2:1160242/76‑1 (MQ=255)                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  >  2:1051661/1‑76 (MQ=255)                                      aGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAgg                                  <  1:641408/76‑1 (MQ=255)                                         gCCAGTGCCTGGATGGGGCTACAGTTAAGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCg                               <  2:1417613/76‑1 (MQ=255)                                           cAGTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGaa                             <  1:602631/76‑1 (MQ=255)                                             gTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                           >  1:1272719/1‑76 (MQ=255)                                             gTGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGt                           <  2:377316/76‑1 (MQ=255)                                              tGCCTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTc                          <  1:858212/76‑1 (MQ=255)                                                ccTGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCgg                        <  1:1019179/76‑1 (MQ=255)                                                  tGAATGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGaa                      >  2:336174/1‑76 (MQ=255)                                                    aaTGGGCGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACATCCTTCGGATGCAGGGCGAAGTCGGaaaa                    <  2:109623/76‑1 (MQ=255)                                                          cGTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctg              >  2:452090/1‑76 (MQ=255)                                                           gTACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgt             >  1:1067054/1‑76 (MQ=255)                                                            tACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt            >  2:498052/1‑76 (MQ=255)                                                            tACAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgtt            <  1:1311823/76‑1 (MQ=255)                                                              cAGTTATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACttctgttct          >  1:1572690/1‑76 (MQ=255)                                                                 ttATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       >  1:141650/1‑76 (MQ=255)                                                                 ttATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTctgttctgtt       >  2:1010792/1‑76 (MQ=255)                                                                  tATGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTa      <  1:30117/76‑1 (MQ=255)                                                                    tGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa    >  2:1323337/1‑76 (MQ=255)                                                                    tGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa    >  1:122744/1‑76 (MQ=255)                                                                    tGAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTaaa    >  2:1255645/1‑76 (MQ=255)                                                                     gAAACCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAt   >  2:605176/1‑76 (MQ=255)                                                                      aaaCCCT‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAAtg  >  1:1131818/1‑76 (MQ=255)                                                                            |                                                                       CTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC‑‑TTTTTTTTCAAGGGCTTCTACAACCTTCGGATGCAGGGCGAAGTCGGAAAACTTCTGTTCTGTTAAATG  >  REL606/3927275‑3927418

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