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%	265.2 / NA	64	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 (33/31);  total (33/31)

AACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGC  >  NC_000913/3090251‑3090312                                |                               aaCCATCTGAACTATCAGGAACTTTCCTGATTTgg                             >  1:1444758/1‑35 (MQ=39) aaCCATCTGAACTATCAGGAACTTTCCTGATTTgg                             <  1:6618847/35‑1 (MQ=39) aaCCATCTGAACTATCAGGAACTTTCCTGATTTgg                             <  1:6618924/35‑1 (MQ=39)  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGc                            >  1:870336/1‑35 (MQ=39)  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGc                            <  1:4647442/35‑1 (MQ=39)  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGc                            >  1:831201/1‑35 (MQ=255)   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCt                           <  1:476971/35‑1 (MQ=255)     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                         <  1:6469545/35‑1 (MQ=39)     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                         >  1:458021/1‑35 (MQ=39)      tCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                        <  1:4867026/35‑1 (MQ=39)      tCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                        <  1:338236/35‑1 (MQ=255)       cTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                       <  1:2139841/35‑1 (MQ=255)       cTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                       <  1:2807394/35‑1 (MQ=39)        tGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                      >  1:7892039/1‑35 (MQ=39)        tGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                      <  1:6930345/35‑1 (MQ=255)         gAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                     >  1:4269428/1‑35 (MQ=255)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    >  1:6795001/1‑35 (MQ=255)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    <  1:832546/35‑1 (MQ=255)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    >  1:5360494/1‑35 (MQ=255)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    >  1:1118040/1‑35 (MQ=255)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    >  1:3729374/1‑35 (MQ=38)          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCa                    >  1:3137538/1‑35 (MQ=38)           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                   <  1:742388/35‑1 (MQ=255)           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                   >  1:7260858/1‑35 (MQ=38)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  >  1:6097542/1‑35 (MQ=38)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  >  1:7520919/1‑35 (MQ=38)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  >  1:2649878/1‑35 (MQ=38)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  <  1:2734395/35‑1 (MQ=255)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  <  1:6581954/35‑1 (MQ=38)            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                  <  1:5087878/35‑1 (MQ=38)             tATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                 >  1:4265656/1‑35 (MQ=38)             tATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                 >  1:2466708/1‑35 (MQ=38)             tATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                 <  1:5746301/35‑1 (MQ=38)             tATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                 <  1:579200/35‑1 (MQ=38)             tATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                 >  1:6193201/1‑35 (MQ=255)              aTCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                <  1:4324952/35‑1 (MQ=38)              aTCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                <  1:1587404/35‑1 (MQ=38)               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCa               <  1:6195547/35‑1 (MQ=255)                cAGGAACTTTCCTGATTTGGCTGATTGCATACCaa              >  1:7150615/1‑35 (MQ=38)                cAGGAACTTTCCTGATTTGGCTGATTGCATACCaa              <  1:7053892/35‑1 (MQ=38)                cAGGAACTTTCCTGATTTGGCTGATTGCATACCaa              <  1:349060/35‑1 (MQ=38)                 aGGAACTTTCCTGATTTGGCTGATTGCATACCaaa             >  1:486363/1‑35 (MQ=38)                 aGGAACTTTCCTGATTTGGCTGATTGCATACCaaa             <  1:774330/35‑1 (MQ=38)                  ggAACTTTCCTGATTTGGCTGATTGCATACCaaaa            >  1:5890304/1‑35 (MQ=38)                  ggAACTTTCCTGATTTGGCTGATTGCATACCaaaa            >  1:5444173/1‑35 (MQ=38)                  ggAACTTTCCTGATTTGGCTGATTGCATACCaaaa            >  1:714337/1‑35 (MQ=255)                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAAc           >  1:1088105/1‑35 (MQ=38)                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACa          >  1:5420310/1‑35 (MQ=255)                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAg         >  1:3459625/1‑35 (MQ=255)                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAg         >  1:348342/1‑35 (MQ=255)                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAg         >  1:5494154/1‑35 (MQ=38)                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCt       <  1:5865660/35‑1 (MQ=255)                        ttCCTGATTTGGCTGATTGCATACCAAAACAGCtt      <  1:5929792/35‑1 (MQ=255)                        ttCCTGATTTGGCTGATTGCATACCAAAACAGCtt      >  1:6802300/1‑35 (MQ=38)                         tCCTGATTTGGCTGATTGCATACCAAAACAGCttt     >  1:4948994/1‑35 (MQ=255)                         tCCTGATTTGGCTGATTGCATACCAAAACAGCttt     <  1:3735097/35‑1 (MQ=255)                         tCCTGATTTGGCTGATTACATACCAAAACAGCttt     <  1:3698123/35‑1 (MQ=25)                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTc    >  1:7078687/1‑35 (MQ=255)                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTc    <  1:4221944/35‑1 (MQ=255)                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCg   >  1:2774455/1‑35 (MQ=38)                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCg   >  1:2183186/1‑35 (MQ=38)                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCg   >  1:5966262/1‑35 (MQ=38)                            tGATTTGGCTGATTGCATACCAAAACAGCTTTCGc  <  1:7101082/35‑1 (MQ=38)                            tGATTTGGCTGATTGCATACCAAAACAGCTTTCGc  <  1:1658816/35‑1 (MQ=255)                            tGATTTGGCTGATTGCATACCAAAACAGCTTTCGc  <  1:2894147/35‑1 (MQ=38)                                |                               AACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGC  >  NC_000913/3090251‑3090312

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