Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,967,964	G→A	A101T (GCC→ACC)	yggG →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,967,964	0	G	A	100.0%	230.8 / NA	60	A101T (GCC→ACC)	yggG	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/25);  total (35/25)

CCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCG  >  REL606/2967893‑2968036                                                                        |                                                                         ccAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCtt                                                                      <  2:244312/76‑1 (MQ=255) ccAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCtt                                                                      >  1:876349/1‑76 (MQ=255)  cgaTGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCttt                                                                     <  2:492604/74‑1 (MQ=255)  cAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTATATTACAAAGTGTATATGGCGAAGGATGTGAACACCttt                                                                     >  1:1264689/1‑76 (MQ=255)  cAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCttt                                                                     >  1:1264674/1‑76 (MQ=255)    aTGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGc                                                                   <  2:98511/76‑1 (MQ=255)    aTGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGc                                                                   >  2:513307/1‑76 (MQ=255)    aTGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGc                                                                   <  1:1140960/76‑1 (MQ=255)    aTGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGc                                                                   >  1:1711078/1‑76 (MQ=255)        gcTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGCACACCTTTGCAATg                                                               >  1:800427/1‑76 (MQ=255)         cTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATgg                                                              >  1:1318363/1‑76 (MQ=255)         cTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATgg                                                              >  1:799360/1‑76 (MQ=255)          tAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGc                                                             >  2:1735609/1‑76 (MQ=255)          tAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGc                                                             >  2:1303558/1‑76 (MQ=255)          tAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGc                                                             >  1:886742/1‑76 (MQ=255)            ggCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTa                                                           <  1:103528/76‑1 (MQ=255)             gCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTaa                                                          <  2:1266985/76‑1 (MQ=255)                 caaTATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGc                                                      >  1:1959756/1‑76 (MQ=255)                      tCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTAt                                                 >  1:1695242/1‑76 (MQ=255)                        aaCGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATcc                                               >  1:796412/1‑76 (MQ=255)                         aCGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCg                                              >  2:334204/1‑76 (MQ=255)                          cGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCgt                                             <  2:1318363/76‑1 (MQ=255)                               aGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCta                                        <  1:274635/76‑1 (MQ=255)                               aGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCta                                        >  1:1126976/1‑76 (MQ=255)                                gCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCtat                                       <  2:1979699/76‑1 (MQ=255)                                 ccGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTTTCCGTGTCtata                                      <  1:978007/76‑1 (MQ=255)                                 ccGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCtata                                      <  1:765196/76‑1 (MQ=255)                                  cGGTAAATTACAACGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAg                                     <  1:1886599/76‑1 (MQ=255)                                  cGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAg                                     <  1:899089/76‑1 (MQ=255)                                   ggTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGc                                    <  1:747998/76‑1 (MQ=255)                                    gTAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCg                                   <  2:417389/76‑1 (MQ=255)                                    gTAAATTACAAAGTGCATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCg                                   >  2:1271120/1‑76 (MQ=255)                                     tAAATTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCgg                                  >  1:212687/1‑76 (MQ=255)                                      aaaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCggg                                 >  2:1300/1‑76 (MQ=255)                                      aaaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCggg                                 >  1:867369/1‑76 (MQ=255)                                      aaaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCggg                                 <  1:1930556/76‑1 (MQ=255)                                      aaaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCggg                                 >  1:1507383/1‑76 (MQ=255)                                       naTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCAGGc                                >  1:791510/2‑76 (MQ=255)                                       aaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGc                                >  1:1052992/1‑76 (MQ=255)                                       aaTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGc                                >  1:1052980/1‑76 (MQ=255)                                        aTTACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCt                               >  1:1001632/1‑76 (MQ=255)                                         ttACAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTg                              >  2:973381/1‑76 (MQ=255)                                           aCAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGAt                            >  1:966271/1‑76 (MQ=255)                                           aCAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGAt                            <  1:669027/76‑1 (MQ=255)                                            cAAAGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATg                           >  1:740581/1‑76 (MQ=255)                                             aaaGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATgg                          <  1:1900661/76‑1 (MQ=255)                                             aaaGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATgg                          <  1:1703375/76‑1 (MQ=255)                                               aGTGTATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGat                        >  2:864616/1‑76 (MQ=255)                                                gtgtATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGata                       >  2:198154/1‑76 (MQ=255)                                                 tgtATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGatat                      >  1:868466/1‑76 (MQ=255)                                                 tgtATATGGCGAAGGATGTGAACACCTTCGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGatat                      <  1:689097/76‑1 (MQ=255)                                                  gtATATGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATatg                     <  1:1172926/76‑1 (MQ=255)                                                   tataTGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATatga                    >  1:553251/1‑76 (MQ=255)                                                    atatGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATatgat                   <  2:1634467/76‑1 (MQ=255)                                                     tatGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATatgatg                  >  2:1585083/1‑76 (MQ=255)                                                       tGGCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGAc                <  2:732234/76‑1 (MQ=255)                                                        ggCGAAGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACg               >  1:379620/1‑76 (MQ=255)                                                             aGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATaa          <  1:1300/76‑1 (MQ=255)                                                             aGGATGTGAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATaa          >  1:153598/1‑76 (MQ=255)                                                                   tgAACACCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATAACGGAGt    <  2:286227/76‑1 (MQ=255)                                                                     aaCGCCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCg  <  2:1221107/76‑1 (MQ=255)                                                                        |                                                                         CCAATGCGCTAGGCAACAATATCAACGGTCAGCCGGTAAATTACAAAGTGTATATGGCGAAGGATGTGAACGCCTTTGCAATGGCTAACGGCTGTATCCGTGTCTATAGCGGGCTGATGGATATGATGACGGATAACGAAGTCG  >  REL606/2967893‑2968036

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