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%	287.5 / NA	72	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/38);  total (34/38)

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090235‑3090327                                                |                                              cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                             >  1:875012/1‑50 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                            <  1:4248449/50‑1 (MQ=255)  aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                            >  1:2256705/1‑50 (MQ=255)   gCTGGTTGCTAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                           >  1:1688674/1‑50 (MQ=255)   gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                           <  1:5079821/50‑1 (MQ=255)     tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                          >  1:3208126/1‑49 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                       <  1:4430118/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                       <  1:792144/50‑1 (MQ=255)       gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                       <  1:2808780/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      >  1:1337187/1‑50 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:5398900/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                      <  1:5659822/50‑1 (MQ=255)        ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtc                                      <  1:5398779/50‑2 (MQ=255)           cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                    <  1:5432291/49‑1 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                   >  1:5502038/1‑49 (MQ=255)            gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                  <  1:1187536/50‑1 (MQ=255)             aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                 <  1:4983706/50‑1 (MQ=255)               aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                >  1:32015/1‑49 (MQ=255)                tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                               >  1:472895/1‑49 (MQ=255)                 aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                              >  1:4766724/1‑49 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:1500962/50‑1 (MQ=255)                  aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                            <  1:3662197/50‑1 (MQ=255)                   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           <  1:255903/50‑1 (MQ=255)                   ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           <  1:5762725/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:3093776/1‑49 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                          <  1:2794817/50‑1 (MQ=255)                    cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                          <  1:3645675/50‑1 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:2428812/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                           >  1:799310/1‑48 (MQ=255)                     aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                         <  1:161138/50‑1 (MQ=255)                        tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      >  1:5991670/1‑50 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                    <  1:2920058/50‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:3252053/48‑1 (MQ=255)                          aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                      <  1:5507178/48‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   <  1:3285673/50‑1 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   >  1:2417219/1‑50 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   >  1:1695874/1‑50 (MQ=255)                           aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                   <  1:4668358/50‑1 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  >  1:5485118/1‑50 (MQ=255)                            cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                  <  1:5826961/50‑1 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:5830838/1‑49 (MQ=255)                              aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                 >  1:453590/1‑49 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                <  1:4124831/49‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa               <  1:3649251/50‑1 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa               >  1:1738023/1‑50 (MQ=255)                               tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa               >  1:3533751/1‑50 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              >  1:4240690/1‑50 (MQ=255)                                cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              >  1:2179540/1‑50 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:599789/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:5612182/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc              <  1:3044488/49‑1 (MQ=255)                                 aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg             >  1:3983979/1‑50 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt           <  1:2620819/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt           <  1:2072734/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt           <  1:2001626/50‑1 (MQ=255)                                   gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt           >  1:733590/1‑50 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc         <  1:4991491/50‑1 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc         >  1:2776968/1‑50 (MQ=255)                                     aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc         >  1:718787/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt        >  1:5016179/1‑50 (MQ=255)                                      cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt        >  1:882100/1‑50 (MQ=255)                                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:3472648/1‑50 (MQ=255)                                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:1572327/1‑50 (MQ=255)                                       tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       >  1:2061587/1‑50 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg      <  1:5214319/49‑1 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg       <  1:1465569/48‑1 (MQ=255)                                         tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc     >  1:938171/1‑50 (MQ=255)                                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc     <  1:2033646/49‑1 (MQ=255)                                          ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc     <  1:5165259/49‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt    <  1:3359859/49‑1 (MQ=255)                                           cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc   >  1:1135938/1‑50 (MQ=255)                                            tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg  >  1:4750436/1‑50 (MQ=255)                                                |                                              CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCG  >  NC_000913/3090235‑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: ‑