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%	257.9 / 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 (35/30);  total (35/30)

AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090236‑3090327                                               |                                              aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                            <  1:1209336/50‑1 (MQ=255)  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                           <  1:3263020/50‑1 (MQ=255)   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                          <  1:3680939/50‑1 (MQ=255)     ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                        >  1:3066237/1‑50 (MQ=255)      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                         <  1:3930363/48‑1 (MQ=255)       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:3114789/50‑1 (MQ=255)       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:2560050/50‑1 (MQ=255)        tGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                     >  1:4905802/1‑50 (MQ=255)        tGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                     >  1:204107/1‑50 (MQ=255)          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:4109388/1‑50 (MQ=255)          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:546317/1‑50 (MQ=255)          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:2217525/1‑50 (MQ=255)           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:3078772/1‑49 (MQ=255)            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                 <  1:2296030/50‑1 (MQ=255)            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                 <  1:3969948/50‑1 (MQ=255)             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                >  1:4026279/1‑50 (MQ=255)             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                <  1:656280/50‑1 (MQ=255)               tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                               >  1:414681/1‑49 (MQ=255)               tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                              <  1:3615654/50‑1 (MQ=255)                aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                              >  1:3924440/1‑49 (MQ=255)                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:3503462/50‑1 (MQ=255)                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:1032538/50‑1 (MQ=255)                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:272344/50‑1 (MQ=255)                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                             >  1:4929652/1‑49 (MQ=255)                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                             <  1:938159/49‑1 (MQ=255)                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:3934017/1‑50 (MQ=255)                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:5297329/1‑49 (MQ=255)                     tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                        >  1:4516945/1‑50 (MQ=255)                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                       >  1:683302/1‑50 (MQ=255)                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                       >  1:2081613/1‑50 (MQ=255)                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                       <  1:4842585/47‑1 (MQ=255)                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   <  1:1744464/50‑1 (MQ=255)                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   >  1:2196372/1‑50 (MQ=255)                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   >  1:5139108/1‑50 (MQ=255)                           cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  >  1:3371078/1‑50 (MQ=255)                           cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  <  1:1877202/50‑1 (MQ=255)                           cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  <  1:4082064/50‑1 (MQ=255)                            tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:2652472/1‑50 (MQ=255)                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:1652597/1‑49 (MQ=255)                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:2203804/1‑49 (MQ=255)                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:3665394/1‑49 (MQ=255)                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 <  1:1131888/48‑1 (MQ=255)                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                <  1:15249/49‑1 (MQ=255)                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                <  1:612776/49‑1 (MQ=255)                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:1417762/49‑1 (MQ=255)                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:2639756/49‑1 (MQ=255)                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:5531143/49‑1 (MQ=255)                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:4303514/49‑1 (MQ=255)                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg             >  1:2894751/1‑50 (MQ=255)                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt            <  1:2104355/50‑1 (MQ=255)                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt            >  1:2735818/1‑50 (MQ=255)                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt            >  1:1924062/1‑50 (MQ=255)                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt           >  1:3593720/1‑50 (MQ=255)                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg          >  1:3507884/1‑50 (MQ=255)                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt        >  1:2315351/1‑50 (MQ=255)                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:4650929/1‑50 (MQ=255)                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:3881788/1‑50 (MQ=255)                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:77556/1‑50 (MQ=255)                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       <  1:3928790/48‑1 (MQ=255)                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       <  1:1328932/48‑1 (MQ=255)                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc     >  1:3660716/1‑50 (MQ=255)                                         ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt    >  1:3760075/1‑50 (MQ=255)                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt    <  1:978811/49‑1 (MQ=255)                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc   <  1:2075304/50‑1 (MQ=255)                                           tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg  >  1:134504/1‑50 (MQ=255)                                               |                                              AGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090236‑3090327

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