Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,775,881	+T	intergenic (‑51/+198)	insJ‑3 ← / ← cysH	IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,775,879	1	.	T	100.0%	232.6 / NA	59	intergenic (‑49/+200)	insJ‑3/cysH	IS150 hypothetical protein/phosphoadenosine phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (32/27);  total (32/27)

TCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑TTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  REL606/2775805‑2775949                                                                            |                                                                       tCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑t                                                                       >  1:926192/1‑76 (MQ=32) tCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑t                                                                       >  2:503250/1‑76 (MQ=32)   aaaaGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑ttt                                                                     <  1:1082774/76‑2 (MQ=21)   aaaaGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑ttt                                                                     <  2:1055168/76‑2 (MQ=21)   aaaaGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑ttt                                                                     >  2:683454/1‑75 (MQ=21)     aaGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑tttga                                                                   >  2:1064787/1‑73 (MQ=21)         gTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGa                                                                >  1:1617889/1‑76 (MQ=255)             tttGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAata                                                            <  1:510039/76‑1 (MQ=255)              ttGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAatat                                                           >  1:60288/1‑76 (MQ=255)              ttGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAatat                                                           >  1:797804/1‑76 (MQ=255)               tGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATc                                                          <  2:748066/76‑1 (MQ=255)               tGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATc                                                          <  2:1540289/76‑1 (MQ=255)               tGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATc                                                          <  2:1608601/76‑1 (MQ=255)               tGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATc                                                          <  2:682746/76‑1 (MQ=255)                ggCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCg                                                         <  1:1324456/76‑1 (MQ=255)                ggCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCg                                                         >  1:561222/1‑76 (MQ=255)                 gCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGc                                                        >  1:1121769/1‑76 (MQ=255)                   tttGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCt                                                      >  1:676695/1‑76 (MQ=255)                   tttGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCt                                                      >  1:128491/1‑76 (MQ=255)                     tGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTgg                                                    >  1:601507/1‑76 (MQ=255)                      gACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTggg                                                   >  1:593497/1‑76 (MQ=255)                         atatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCt                                                >  1:964409/1‑76 (MQ=255)                         atatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCt                                                >  1:81464/1‑76 (MQ=255)                          tatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCtt                                               >  1:1049718/1‑76 (MQ=255)                          tatTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCtt                                               <  2:184731/76‑1 (MQ=255)                            ttAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCtttt                                             >  1:254766/1‑76 (MQ=255)                             tAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATGGAAATATCGCCTGGGGCTTTTg                                            <  2:969631/76‑1 (MQ=18)                                 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                        >  1:1262886/1‑76 (MQ=255)                                 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                        <  2:1372700/76‑1 (MQ=255)                                 ggCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGGCCa                                        <  2:615729/76‑1 (MQ=21)                                  gCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt                                       >  1:435998/1‑76 (MQ=255)                                   cTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                      >  1:725878/1‑76 (MQ=255)                                    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                     >  1:1389739/1‑76 (MQ=255)                                    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                     >  1:985298/1‑76 (MQ=255)                                    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                     >  1:1239135/1‑76 (MQ=255)                                    tATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                     <  1:1486233/76‑1 (MQ=255)                                     aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    >  1:1728612/1‑76 (MQ=255)                                     aTTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    <  1:1233095/76‑1 (MQ=255)                                      ttCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTgg                                   <  2:486870/76‑1 (MQ=255)                                        ccATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                 <  1:1396440/76‑1 (MQ=255)                                        ccATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                 <  1:1354388/76‑1 (MQ=255)                                        ccATTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                 <  1:852277/76‑1 (MQ=255)                                          aTTTCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  1:899441/1‑76 (MQ=255)                                           tttCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                              <  2:1173381/76‑1 (MQ=255)                                            ttCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                             <  1:1489012/76‑1 (MQ=255)                                            ttCATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                             <  1:1343416/76‑1 (MQ=255)                                              cATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           >  2:557129/1‑76 (MQ=255)                                              cATCGTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           >  2:537892/1‑76 (MQ=255)                                                  gTCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                       >  1:16525/1‑76 (MQ=255)                                                   tCCAACAAAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      <  2:319203/76‑1 (MQ=255)                                                       acaaAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                  >  1:1388083/1‑76 (MQ=255)                                                        caaAATGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                 >  1:279107/1‑76 (MQ=255)                                                         aaaaTGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                <  2:1374920/76‑1 (MQ=255)                                                         aaaaTGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                >  1:603969/1‑76 (MQ=255)                                                           aaTGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              >  2:11482/1‑76 (MQ=255)                                                            aTGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc             >  2:56989/1‑76 (MQ=255)                                                            aTGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc             >  1:1191952/1‑76 (MQ=255)                                                             tGGGTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  1:1501172/76‑1 (MQ=255)                                                              gggTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg           <  1:382428/76‑1 (MQ=255)                                                               ggTGCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa          <  1:1178571/76‑1 (MQ=255)                                                                  gCAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg       <  2:228749/76‑1 (MQ=255)                                                                   cAGTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc      <  1:1640726/76‑1 (MQ=255)                                                                     gTACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt    >  1:770037/1‑76 (MQ=255)                                                                      tACAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc   <  2:632109/76‑1 (MQ=255)                                                                       aCAGGTTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCa  >  1:1403032/1‑76 (MQ=255)                                                                            |                                                                       TCAAAAGGGTATTTTGGCTTTGACATATTAGGGGCTATTCCATTTCATCGTCCAACAAAATGGGTGCAGTACAGG‑TTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  REL606/2775805‑2775949

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