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%	254.3 / NA	65	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 (34/31);  total (34/31)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328                                                |                                               cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:4865903/50‑1 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:4804451/50‑1 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:4469209/50‑1 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:3652480/50‑1 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:1825644/50‑1 (MQ=255)    cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:881259/1‑50 (MQ=255)    cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:788752/1‑50 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:199967/1‑49 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:4172266/1‑50 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:4014260/1‑50 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:1052026/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       >  1:3878371/1‑50 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:1056208/1‑49 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:1137778/1‑49 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:3439847/49‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:690477/50‑1 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:1010793/1‑50 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 <  1:5056904/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:1864082/1‑49 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:5082126/50‑1 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                             <  1:4001212/50‑1 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:4140854/1‑49 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:4474098/1‑49 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:1119716/50‑1 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                          <  1:4099439/50‑1 (MQ=255)                      tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  1:2263292/1‑50 (MQ=255)                      tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  1:3073290/1‑50 (MQ=255)                          atcTATCTGGCGCTTTCCTGATTTCGCTGATTGCATACCAAAACAGCt                       <  1:1558640/46‑1 (MQ=19)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:446263/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:92798/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  1:3443915/49‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  1:4304558/49‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  1:3085793/48‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  1:5092209/1‑50 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  1:5107630/1‑50 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  1:992164/1‑50 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   >  1:1465211/1‑50 (MQ=255)                             tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:3738957/1‑50 (MQ=255)                             tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:3004388/1‑50 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:1805909/1‑49 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:1606294/50‑1 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:4672306/1‑50 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:2145955/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:1702556/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:1378064/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              <  1:2399037/50‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              <  1:1766467/50‑1 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             <  1:3450719/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:2347261/50‑1 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           >  1:4239327/1‑50 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  1:4910208/50‑1 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          >  1:1788771/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         >  1:4171612/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:1711374/50‑1 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         >  1:3819193/1‑50 (MQ=255)                                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        >  1:2661370/1‑50 (MQ=255)                                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        >  1:1815548/1‑50 (MQ=255)                                        ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  1:352362/1‑50 (MQ=255)                                        ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  1:654725/1‑50 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:4482991/48‑1 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:4193043/1‑50 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:1277803/50‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    >  1:4150877/1‑50 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:3448593/50‑1 (MQ=255)                                             gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  1:665747/1‑50 (MQ=255)                                                |                                               CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328

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