Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	382,261	A→G	I639T (ATT→ACT)	sbcC ←	exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	382,261	0	A	G	100.0%	253.3 / NA	59	I639T (ATT→ACT)	sbcC	exonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/32);  total (27/32)

TTCTTCATCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAATTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCTGA  >  REL606/382188‑382333                                                                          |                                                                         ttcttcATCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGtt                                                                        <  1:1302402/76‑1 (MQ=255)  tcttcATCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGttt                                                                       <  2:1575614/76‑1 (MQ=255)  tcttcATCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGttt                                                                       <  2:678369/76‑1 (MQ=255)    ttcATCTACCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGc                                                                     >  1:1294131/1‑76 (MQ=255)       aTCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGt                                                                  <  1:325016/76‑1 (MQ=255)       aTCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGt                                                                  <  1:1886140/76‑1 (MQ=255)           tCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGAt                                                              <  1:1132386/76‑1 (MQ=255)           tCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGAt                                                              <  1:933167/76‑1 (MQ=255)           tCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGAt                                                              >  2:585396/1‑76 (MQ=255)             cTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATAc                                                            >  2:1950270/1‑76 (MQ=255)             cTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTACTGTTGATAc                                                            <  1:1050319/76‑1 (MQ=255)              tgtgGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACt                                                           >  1:1232135/1‑76 (MQ=255)              tgtgGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACt                                                           <  2:1373773/76‑1 (MQ=255)               gtgGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTg                                                          >  2:1863903/1‑76 (MQ=255)               gtgGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTg                                                          <  2:977742/76‑1 (MQ=255)                tgGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTgg                                                         <  1:642125/76‑1 (MQ=255)                    aaTGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATaa                                                     <  2:614580/76‑1 (MQ=255)                    aaTGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATaa                                                     >  2:1489475/1‑76 (MQ=255)                     aTGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAAt                                                    <  2:863335/76‑1 (MQ=255)                        tCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTg                                                 <  1:1241233/76‑1 (MQ=255)                           gTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATAGTGGATAATTTGCTg                                              <  1:1073225/76‑1 (MQ=255)                           gTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTg                                              <  2:1646093/76‑1 (MQ=255)                            tGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGa                                             >  1:1499919/1‑76 (MQ=255)                              cATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGAtt                                           >  2:987770/1‑76 (MQ=255)                              cATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGAtt                                           <  2:1288343/76‑1 (MQ=255)                               aTAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTa                                          <  1:1188746/76‑1 (MQ=255)                               aTAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTa                                          >  1:712626/1‑76 (MQ=255)                                tagcccGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTAt                                         <  1:487290/73‑1 (MQ=255)                                taccccGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTAt                                         <  1:253097/73‑1 (MQ=255)                                tAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTAt                                         >  2:665889/1‑76 (MQ=255)                                tAACCCGTCAATGTCGTTAAAAGTAGTAGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTAt                                         <  1:1756849/76‑1 (MQ=255)                                  aCCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATgc                                       <  2:1781256/76‑1 (MQ=255)                                   cccGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATgcg                                      <  2:924947/76‑1 (MQ=255)                                        cAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCa                                 <  2:510445/76‑1 (MQ=255)                                           tGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATc                              >  2:686684/1‑76 (MQ=255)                                           tGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATc                              <  1:1680165/76‑1 (MQ=255)                                            gTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCt                             <  1:1189801/76‑1 (MQ=255)                                            gTCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCt                             >  1:374615/1‑76 (MQ=255)                                             tCGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTg                            >  2:1657620/1‑76 (MQ=255)                                              cGTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGc                           >  1:27569/1‑76 (MQ=255)                                               gTTAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGcc                          >  1:1238075/1‑76 (MQ=255)                                                ttAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTCATTATGCGCGGCAATCTGccc                         >  1:696032/1‑76 (MQ=255)                                                 tAAAAGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCt                        >  1:58013/1‑76 (MQ=255)                                                   aaaGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTACGCGCGGCAATCTGCCCTTg                      <  2:662565/76‑1 (MQ=255)                                                    aaGTAGTTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGt                     >  2:773639/1‑76 (MQ=255)                                                     agtagtTGCTGGCGTTGTTCAGCTTGCTGTTGGTACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTa                    >  2:750133/1‑76 (MQ=255)                                                      gtagtTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTaa                   <  1:1922567/76‑1 (MQ=255)                                                      gtagtTGCCGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTaa                   >  1:952200/1‑76 (MQ=255)                                                        agtTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAAtt                 >  2:188445/1‑76 (MQ=255)                                                        agtTGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAAtt                 >  1:572506/1‑76 (MQ=255)                                                           tGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCAt              <  1:839164/76‑1 (MQ=255)                                                           tGCTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCAt              >  2:295698/1‑76 (MQ=255)                                                             cTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGc            <  2:1426046/76‑1 (MQ=255)                                                             cTGGCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGc            >  2:907492/1‑76 (MQ=255)                                                               ggCGTTGTTCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCg          >  1:45374/1‑76 (MQ=255)                                                                     gttCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCt    >  1:1291438/1‑76 (MQ=255)                                                                     gttCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCt    >  1:1303947/1‑76 (MQ=255)                                                                     gttCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCt    <  1:826023/76‑1 (MQ=255)                                                                      ttCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCTg   >  2:1277693/1‑76 (MQ=255)                                                                       tCAGTTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCTGa  <  2:341307/76‑1 (MQ=255)                                                                          |                                                                         TTCTTCATCTTCCTGTGGCAATGTCAGTGCATAACCCGTCAATGTCGTTAAAAGTAGTTGCTGGCGTTGTTCAATTTGCTGTTGATACTGGATAATTTGCTGATTATGCGCGGCAATCTGCCCTTGTAATTCATGCCGTTGGCTGA  >  REL606/382188‑382333

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