Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,090,282	C→T	intergenic (+30/‑65)	yggI → / → endA	conserved protein YggI/DNA‑specific endonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,090,282	0	C	T	100.0%	354.3 / NA	86	intergenic (+30/‑65)	yggI/endA	conserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/50);  total (36/50)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090235‑3090325                                                |                                            cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           >  1:4230270/1‑50 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           <  1:6698626/50‑1 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           >  1:3087442/1‑50 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                           <  1:2886329/50‑1 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                          >  1:4294919/1‑50 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                          >  1:5165216/1‑50 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                          >  1:2434811/1‑50 (MQ=255)   gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                         <  1:678652/50‑1 (MQ=255)   gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                         <  1:4715114/50‑1 (MQ=255)    cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                        <  1:2055422/50‑1 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                        >  1:2288697/1‑49 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                        >  1:6559833/1‑49 (MQ=255)      ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                      >  1:3788134/1‑50 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                       <  1:773609/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                       <  1:4141871/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                       <  1:5008182/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                     >  1:1865936/1‑50 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                    <  1:5981406/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                    <  1:1800869/50‑1 (MQ=255)         tGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   >  1:5374572/1‑50 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   <  1:6132823/49‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   <  1:3124032/49‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                   <  1:3212251/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                  <  1:2978168/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                  <  1:548084/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                 <  1:2038801/50‑1 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                 >  1:5957549/1‑49 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                <  1:2206029/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                               <  1:3852925/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                               <  1:3744364/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                               <  1:4545740/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                               >  1:4434731/1‑50 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                               <  1:5848685/50‑1 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:3822724/1‑50 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:4114193/1‑50 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                              >  1:2408888/1‑50 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                             <  1:3074739/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                             <  1:4663824/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                             <  1:5123357/50‑1 (MQ=255)                 aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                            >  1:1276359/1‑49 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                          <  1:813208/50‑1 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                          <  1:31476/50‑1 (MQ=255)                   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:2109624/1‑50 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:6551031/1‑49 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:3212944/1‑49 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                        <  1:3860983/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                        <  1:1756920/50‑1 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:100295/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:622880/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                         >  1:5222165/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                       <  1:2639137/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                      <  1:5879405/46‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:2495076/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  >  1:5935016/1‑50 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:6059173/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:5413982/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:6622314/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                  <  1:1378214/50‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                 <  1:4016553/50‑1 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                <  1:2651812/50‑1 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                <  1:3489250/50‑1 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                >  1:5908404/1‑50 (MQ=255)                             tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc               >  1:4334403/1‑50 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc            >  1:4025856/1‑50 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc            >  1:1435713/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc            <  1:4916966/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg           <  1:6420652/50‑1 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          >  1:375617/1‑50 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          <  1:4346343/50‑1 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          <  1:335103/50‑1 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt          <  1:2156577/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt         <  1:6478579/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt         >  1:6264245/1‑50 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg        >  1:1730535/1‑50 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg        <  1:6548169/50‑1 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg        >  1:1583424/1‑50 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc       >  1:1888188/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt      <  1:5911681/50‑1 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt      >  1:2631179/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt      >  1:4306193/1‑50 (MQ=255)                                        ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg    >  1:1504893/1‑50 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg    <  1:3033610/49‑1 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg     <  1:2146712/48‑1 (MQ=255)                                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt  >  1:6411640/1‑50 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt  <  1:3481695/49‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt  <  1:3099269/49‑1 (MQ=255)                                                |                                            CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCT  >  NC_000913/3090235‑3090325

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