Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,194,365	C→T	T214M (ACG→ATG)	metA →	homoserine O‑succinyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,194,365	0	C	T	100.0%	264.9 / NA	67	T214M (ACG→ATG)	metA	homoserine O‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/35);  total (32/35)

GCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATA  >  REL606/4194294‑4194437                                                                        |                                                                         gCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGgaa                                                                      <  1:985097/76‑1 (MQ=255)  cATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGgaag                                                                     <  1:937111/76‑1 (MQ=255)   aTTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGgaaga                                                                    >  2:949990/1‑76 (MQ=255)   aTTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGgaaga                                                                    >  1:77864/1‑76 (MQ=255)   aTTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGgaaga                                                                    >  2:1729441/1‑76 (MQ=255)      cgcgCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAgg                                                                 >  1:1075282/1‑76 (MQ=255)      cgcgCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAgg                                                                 <  1:1264701/76‑1 (MQ=255)       gcgcTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAggg                                                                <  1:1658024/76‑1 (MQ=255)         gcTATGCTGACGTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGa                                                              <  2:809076/76‑1 (MQ=255)          cTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGAt                                                             <  1:43212/76‑1 (MQ=255)          cTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGAt                                                             >  1:1184069/1‑76 (MQ=255)          cTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGAt                                                             >  1:1651559/1‑76 (MQ=255)             tGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCa                                                          >  2:833786/1‑76 (MQ=255)              gCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCat                                                         <  1:1532896/76‑1 (MQ=255)              gCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCat                                                         >  1:304864/1‑76 (MQ=255)              gCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCat                                                         >  1:210208/1‑76 (MQ=255)                      tCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGttt                                                 <  1:1063437/76‑1 (MQ=255)                      tCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGttt                                                 <  1:548367/76‑1 (MQ=255)                       ccGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGtttn                                                <  1:225087/76‑2 (MQ=255)                        cGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGc                                               <  1:381244/76‑1 (MQ=255)                           cAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATATGTTTGCCAg                                            <  1:458739/76‑1 (MQ=255)                            aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGt                                           <  2:1513106/76‑1 (MQ=255)                            aGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGt                                           >  2:1585449/1‑76 (MQ=255)                              cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTaa                                         <  1:936457/76‑1 (MQ=255)                              cGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTaa                                         >  1:1515344/1‑76 (MQ=255)                               gTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTaaa                                        >  1:1224461/1‑76 (MQ=255)                                ttGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAg                                       <  2:1569970/76‑1 (MQ=255)                                ttGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAg                                       >  1:293489/1‑76 (MQ=255)                                ttGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAg                                       >  2:843614/1‑76 (MQ=255)                                    ttCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATaa                                   >  1:1457441/1‑76 (MQ=255)                                     tCGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAg                                  <  1:1107638/76‑1 (MQ=255)                                      cGTGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAgc                                 >  1:892414/1‑76 (MQ=255)                                        tGATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAgcgc                               >  2:676433/1‑76 (MQ=255)                                         gATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCa                              >  1:738561/1‑76 (MQ=255)                                         gATTACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCa                              >  1:8904/1‑76 (MQ=255)                                            tACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTg                           >  2:1137102/1‑76 (MQ=255)                                            tACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTg                           <  1:689912/76‑1 (MQ=255)                                            tACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTg                           <  1:1244896/76‑1 (MQ=255)                                            tACACCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTg                           <  1:315951/76‑1 (MQ=255)                                             acacCGATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGc                          >  1:1092300/1‑76 (MQ=255)                                                  gATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                     <  1:1053654/76‑1 (MQ=255)                                                  gATCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtg                     <  1:1368343/76‑1 (MQ=255)                                                   aTCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtgt                    <  1:746523/76‑1 (MQ=255)                                                    tCTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTtgtg                   <  2:892414/76‑1 (MQ=255)                                                     cTGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGa                  >  2:1052132/1‑76 (MQ=255)                                                      tGGAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGAc                 <  1:801877/76‑1 (MQ=255)                                                        gAAATTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACgg               >  1:190264/1‑76 (MQ=255)                                                         aaaTTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACggg              >  1:1635082/1‑76 (MQ=255)                                                           aTTCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGcc            >  2:596690/1‑76 (MQ=255)                                                            ttCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCa           >  1:819045/1‑76 (MQ=255)                                                            ttCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCa           <  2:1420316/76‑1 (MQ=255)                                                             tCTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCAt          <  2:1742393/76‑1 (MQ=255)                                                              cTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATc         <  1:1090437/76‑1 (MQ=255)                                                              cTGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATc         >  2:1217689/1‑76 (MQ=255)                                                               tGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        <  2:1077730/76‑1 (MQ=255)                                                               tGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        >  2:1548373/1‑76 (MQ=255)                                                               tGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        <  1:393428/76‑1 (MQ=255)                                                               tGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        <  1:595649/76‑1 (MQ=255)                                                               tGGCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        >  1:100595/1‑76 (MQ=255)                                                               nggCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATcc        >  1:963898/2‑76 (MQ=255)                                                                ggCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATccc       <  2:53449/76‑1 (MQ=255)                                                                 gCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg      <  1:184915/76‑1 (MQ=255)                                                                 gCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg      <  1:1630109/76‑1 (MQ=255)                                                                 gCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg      >  1:314388/1‑76 (MQ=255)                                                                 gCAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCg      <  2:679605/76‑1 (MQ=255)                                                                  cAGAGATGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGa     >  2:1358390/1‑76 (MQ=255)                                                                     agaTGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAata  <  2:1321242/76‑1 (MQ=255)                                                                        |                                                                         GCATTCGCGCTATGCTGACTTTCCGGCAGCGTTGATTCGTGATTACACCGATCTGGAAATTCTGGCAGAGACGGAAGAAGGGGATGCATATCTGTTTGCCAGTAAAGATAAGCGCATTGCCTTTGTGACGGGCCATCCCGAATA  >  REL606/4194294‑4194437

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