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%	207.1 / NA	52	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 (19/33);  total (19/33)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328                                                |                                               cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              >  1:1618494/1‑50 (MQ=255)    cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           <  1:1799272/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                          <  1:1230188/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                          <  1:1669639/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:3863481/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:1665044/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:228764/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:1862741/50‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      <  1:3152073/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:3352656/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:2771407/1‑50 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:3487966/1‑49 (MQ=255)             aaaTAACCATCTGAACTATCAGAAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:2360985/50‑1 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:3810529/1‑50 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                                <  1:285343/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:2672904/1‑49 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:543010/50‑1 (MQ=255)                 aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               >  1:1655369/1‑49 (MQ=255)                   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  1:2162714/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:3477579/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:1072747/50‑1 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:1217069/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:533215/1‑48 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         <  1:1025140/46‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        <  1:4274341/47‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        <  1:3538274/47‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:3307639/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:4024276/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:518517/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  1:2997535/48‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  1:2771291/50‑1 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   <  1:2140774/50‑1 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:2003021/1‑49 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                 <  1:2101468/49‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:3160153/50‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  1:303763/1‑50 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:2952303/50‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  1:522317/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:2891871/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:3856602/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:2484866/1‑50 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             >  1:3545599/1‑50 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:3576586/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:1139899/50‑1 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           >  1:1270010/1‑50 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          >  1:2745457/1‑50 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     <  1:203302/49‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:2916064/50‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:1368196/50‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    >  1:781610/1‑50 (MQ=255)                                            tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg   >  1:4094499/1‑50 (MQ=255)                                             gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  1:521557/1‑50 (MQ=255)                                                |                                               CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328

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