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%	342.1 / NA	90	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 (44/46);  total (44/46)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTC  >  NC_000913/3090235‑3090326                                                |                                             cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                            >  1:1117849/1‑50 (MQ=255) cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                            <  1:3430002/50‑1 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                           <  1:6480483/50‑1 (MQ=255)   gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                          <  1:5499724/50‑1 (MQ=255)    cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                         >  1:5682595/1‑50 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                         >  1:6100862/1‑49 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                         >  1:4361080/1‑49 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                         >  1:6631304/1‑49 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                        <  1:3617090/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                        <  1:876617/48‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                       <  1:5670317/49‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                      <  1:6149685/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                      <  1:1744034/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                      >  1:3093200/1‑50 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                     <  1:1952861/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                     <  1:6552869/50‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                    <  1:6930931/49‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                    <  1:1193572/49‑1 (MQ=255)          gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                    <  1:2854448/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                   <  1:6963624/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                   <  1:5608202/49‑1 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                  >  1:7003403/1‑50 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                  >  1:3662075/1‑49 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                  >  1:353344/1‑49 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                  >  1:4821858/1‑49 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                 <  1:6091133/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                >  1:5859566/1‑50 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                <  1:3521972/50‑1 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               <  1:3088383/50‑1 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               >  1:3994955/1‑50 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               >  1:4084088/1‑50 (MQ=255)              aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               <  1:5334280/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                              <  1:4481831/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               >  1:7327381/1‑49 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                               >  1:2929636/1‑49 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                              >  1:2125527/1‑49 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                             <  1:1351284/50‑1 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                             <  1:2548977/50‑1 (MQ=255)                 aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                             >  1:4946688/1‑49 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                            >  1:1382999/1‑49 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                            <  1:5849455/49‑1 (MQ=255)                   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                          <  1:4675768/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                         <  1:843262/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                         <  1:7315908/50‑1 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                          >  1:261573/1‑48 (MQ=255)                       cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                      >  1:3991107/1‑50 (MQ=255)                         gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                    >  1:2126810/1‑50 (MQ=255)                         gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                    >  1:7164956/1‑50 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                       <  1:6306338/46‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                   <  1:4095887/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                   >  1:4868405/1‑50 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                    <  1:1607152/49‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                    <  1:147459/49‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                    <  1:2289744/49‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                     <  1:471871/48‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                  >  1:5932123/1‑50 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                 >  1:3990550/1‑50 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                 <  1:391059/50‑1 (MQ=255)                             tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                >  1:3064049/1‑50 (MQ=255)                             tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                >  1:1710415/1‑50 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                >  1:6425706/1‑49 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                >  1:5728221/1‑49 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                >  1:1672700/1‑49 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt               <  1:6271014/49‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa              <  1:598953/50‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTGTCGCTa              <  1:3912681/50‑1 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc             >  1:4597158/1‑50 (MQ=255)                                 aGGAGCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTGGCTAc             <  1:7120657/49‑1 (MQ=38)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg            >  1:7163679/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg            <  1:3622026/50‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg            <  1:3711823/50‑1 (MQ=255)                                 aGGAACTTTCCGGATTGGGCTGATTGCATACCAAAACAGCTTTCGCTAc             <  1:4920206/49‑1 (MQ=255)                                  ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt           >  1:3194780/1‑50 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt          <  1:5923614/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt          >  1:1115909/1‑50 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg         <  1:6805702/50‑1 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg         >  1:3382551/1‑50 (MQ=255)                                    aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg         <  1:1360994/50‑1 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc        >  1:2031814/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt       >  1:3179718/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt       >  1:1600332/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt       >  1:1464608/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt       >  1:896843/1‑50 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg     <  1:5727525/49‑1 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc    >  1:6147222/1‑50 (MQ=255)                                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt   >  1:6903831/1‑50 (MQ=255)                                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt   >  1:5645516/1‑50 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt   <  1:6863878/49‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc  <  1:3953553/50‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc  >  1:6933772/1‑50 (MQ=255)                                                |                                             CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTC  >  NC_000913/3090235‑3090326

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