Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,142,600	G→A	A254V (GCG→GTG)	ygjD ←	O‑sialoglycoprotein endopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,142,600	0	G	A	100.0%	278.4 / NA	80	A254V (GCG→GTG)	ygjD	O‑sialoglycoprotein endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (49/31);  total (49/31)

CTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCT  >  REL606/3142527‑3142672                                                                          |                                                                         cTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGCTCCAGCAcc                                                                        >  2:2217704/1‑76 (MQ=255) cTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCAGGACCAGTCGCTTAAAGCCCGTCTGATCCAGCAcc                                                                        >  2:1124868/1‑76 (MQ=255)   tCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCg                                                                      >  2:2459721/1‑76 (MQ=255)     gCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCt                                                                    >  1:5740157/1‑76 (MQ=255)         gTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCa                                                                >  1:1049255/1‑76 (MQ=255)          tAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCAc                                                               >  1:6137046/1‑76 (MQ=255)           aaCGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACt                                                              >  2:3323927/1‑76 (MQ=255)             cGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTa                                                            >  2:4824748/1‑76 (MQ=255)             cGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTa                                                            >  1:3047167/1‑76 (MQ=255)              gTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTaa                                                           >  1:3811669/1‑76 (MQ=255)               tACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAat                                                          >  1:5994810/1‑76 (MQ=255)                aCGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAatc                                                         >  1:2940282/1‑76 (MQ=255)                    ttAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAAtcatca                                                     >  2:1544579/1‑76 (MQ=255)                    ttAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAAtcatca                                                     >  2:3470822/1‑76 (MQ=255)                       gCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCg                                                  >  2:3351750/1‑76 (MQ=255)                         aCTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTa                                                >  1:1316567/1‑76 (MQ=255)                         aCTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTa                                                >  1:5861694/1‑76 (MQ=255)                          cTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTAt                                               >  1:350195/1‑76 (MQ=255)                            cACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCg                                             >  2:5167539/1‑76 (MQ=255)                             aCGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGa                                            >  2:1469037/1‑76 (MQ=255)                               gccgccCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGacc                                          >  1:892723/1‑76 (MQ=255)                               cccgcccgcCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGacc                                          <  1:4776507/75‑1 (MQ=255)                                 gccccgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGaccac                                        <  2:482952/74‑1 (MQ=255)                                 cgcccgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGaccac                                        >  1:1296188/1‑76 (MQ=255)                                   cccgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCg                                      <  1:335006/76‑1 (MQ=255)                                    ccgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCCaccccccc                                     >  1:1836213/1‑74 (MQ=255)                                     cgccATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCa                                    >  1:1663525/1‑76 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt                                >  2:4607261/1‑76 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt                                <  1:5907980/76‑1 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt                                >  2:4187637/1‑76 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt                                <  2:341526/76‑1 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCtt                                <  2:1754183/76‑1 (MQ=255)                                         aTGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCACCtt                                >  2:1403226/1‑76 (MQ=255)                                          tGACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTc                               <  1:4397760/76‑1 (MQ=255)                                           gACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCa                              <  1:1504414/76‑1 (MQ=255)                                           gACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCa                              <  1:4485820/76‑1 (MQ=255)                                           gACCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCa                              <  2:3105136/76‑1 (MQ=255)                                            aCCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaa                             <  1:396847/76‑1 (MQ=255)                                            aCCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaa                             >  2:1799750/1‑76 (MQ=255)                                            aCCAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaa                             <  2:3047167/76‑1 (MQ=255)                                             ccAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCaaa                            >  1:2548339/1‑76 (MQ=255)                                             ccAGTCCCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTACCGACGACCGCATCTCCaaa                            >  2:1410568/1‑76 (MQ=255)                                              cAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAg                           <  2:1023806/76‑1 (MQ=255)                                              cAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAg                           <  2:3074022/76‑1 (MQ=255)                                              cAGTCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAg                           <  1:1198382/76‑1 (MQ=255)                                                 tCGCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGgcg                        >  1:4145893/1‑76 (MQ=255)                                                   gCTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGcgcg                      >  2:771310/1‑76 (MQ=255)                                                    cTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCgg                     >  1:5429405/1‑76 (MQ=255)                                                    cTTAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCgg                     >  2:2049901/1‑76 (MQ=255)                                                     ttAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCggg                    >  2:4835573/1‑76 (MQ=255)                                                     ttAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCggg                    >  2:3796721/1‑76 (MQ=255)                                                     ttAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCggg                    >  1:2914764/1‑76 (MQ=255)                                                     ttAAAGCCCGGCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCggg                    <  1:5950825/76‑1 (MQ=255)                                                      tAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGc                   <  1:609981/76‑1 (MQ=255)                                                      tAAAGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGc                   <  1:4747963/76‑1 (MQ=255)                                                       aaaGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCg                  <  1:1905636/76‑1 (MQ=255)                                                       aaaGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCg                  >  2:3733050/1‑76 (MQ=255)                                                       aaaGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCg                  >  2:2329590/1‑76 (MQ=255)                                                        aaGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGa                 >  2:598923/1‑76 (MQ=255)                                                        aaGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGa                 <  1:3237632/76‑1 (MQ=255)                                                         aGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGat                >  1:6031714/1‑76 (MQ=255)                                                         aGCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGat                <  1:4354185/76‑1 (MQ=255)                                                          gCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGata               >  2:5116140/1‑76 (MQ=255)                                                          gCCCGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGata               <  1:1749864/76‑1 (MQ=255)                                                          gCCCGTCTGACCCAGCCCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGata               <  1:4130236/76‑1 (MQ=255)                                                           cccGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGatat              <  1:4150003/76‑1 (MQ=255)                                                            ccGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGCTATc             >  1:2797199/1‑76 (MQ=255)                                                            ccGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATc             <  1:1593811/76‑1 (MQ=255)                                                            ccGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATc             >  1:3022876/1‑76 (MQ=255)                                                            ccGTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATc             <  1:5445145/76‑1 (MQ=255)                                                            ccGTCTGACCCCGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATc             <  2:1196180/76‑1 (MQ=255)                                                              gtgtGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAg           <  2:1368989/73‑1 (MQ=255)                                                              gTCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAg           <  1:5905219/76‑1 (MQ=255)                                                               tCTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGc          >  2:701189/1‑76 (MQ=255)                                                                cTGATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCa         >  1:3040384/1‑76 (MQ=255)                                                                 tGATTCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAAGCGCGAGGGATATCAGCAc        >  2:3991557/1‑76 (MQ=255)                                                                  gATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCAcg       <  2:4161165/76‑1 (MQ=255)                                                                  gATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCAcg       <  1:5571375/76‑1 (MQ=255)                                                                  gATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCAcg       >  1:3844570/1‑76 (MQ=255)                                                                  gATCCAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGACATATAAGCAcg       >  1:1702773/1‑76 (MQ=255)                                                                      cAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGTGATATCAGCACGCGTc   <  2:449110/76‑1 (MQ=255)                                                                      cAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTc   <  2:2013797/76‑1 (MQ=255)                                                                      cAGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTc   >  2:6022350/1‑76 (MQ=255)                                                                       aGCACCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCt  <  1:3015123/76‑1 (MQ=255)                                                                       aGCACCCGCTTGCACTTAACCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCt  <  1:3603496/76‑1 (MQ=255)                                                                          |                                                                         CTTCGCCCGTAACGTACGGTTAGCACTCACGCCGCCCGCCATGACCAGTCGCTTAAAGCCCGTCTGATCCAGCGCCCGCTTGCACTTAATCATCAGCGTATCGACCACCGCATCTTCAAAGGCGCGGGCGATATCAGCACGCGTCT  >  REL606/3142527‑3142672

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