Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,399,985	(C)7→8	coding (252/516 nt)	ogt ←	O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,399,978	1	.	C	100.0%	95.6 / NA	28	T87A (ACG→GCG)	ogt	O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (15/13);  total (15/13)

TGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT  >  NC_000913/1399840‑1400100                                                                                                                                            |                                                                                                                            tGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGt                                                                                                                              <  2:368672/139‑1 (MQ=255)     aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                        <  1:85644/139‑1 (MQ=255)     aTTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccc                                                                                                                        >  1:141614/1‑139 (MQ=255)       tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccccc                                                                                                                      >  1:467361/1‑139 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  1:129443/139‑3 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  2:141614/139‑3 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  1:803271/139‑3 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  1:683983/139‑3 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  2:384706/139‑3 (MQ=255)          aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  2:559222/139‑3 (MQ=255)                     gcgccAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCCGTAGCAGTggg                                                                                                          <  1:668538/136‑1 (MQ=255)                       aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAg                                                                                                       <  1:173321/139‑1 (MQ=255)                            cGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTAt                                                                                                  <  2:467361/139‑1 (MQ=255)                               cACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                               <  2:267802/139‑1 (MQ=255)                                              gACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTAc                                                                                >  2:57113/1‑139 (MQ=255)                                                  gCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGc                                                                            >  1:548697/1‑139 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGa                                                                                                         >  1:543962/1‑109 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGa                                                                                                         <  2:543962/109‑1 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                           <  1:689447/139‑1 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   >  1:104565/1‑96 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑ccccccccg                                                                                                                   <  2:104565/98‑3 (MQ=255)                                                   cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTCTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                           >  2:689447/1‑139 (MQ=255)                                                                    ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                          <  2:887701/139‑1 (MQ=255)                                                                    ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCtt                                                          >  1:686434/1‑139 (MQ=255)                                                                                       tGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                       <  1:57113/139‑1 (MQ=255)                                                                                          ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtg                                    >  2:946127/1‑139 (MQ=255)                                                                                                   aTAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                       >  2:285255/1‑127 (MQ=255)                                                                                                   aTAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGa                                       <  1:285255/127‑1 (MQ=255)                                                                                                     aGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAAt                         >  2:151762/1‑139 (MQ=255)                                                                                                      gTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATg                        <  1:516823/139‑1 (MQ=255)                                                                                                       tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccccccccg                                                                                                                  <  2:986824/47‑4 (MQ=255)                                                                                                       tGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cccccccccg                                                                                                                  >  1:986824/1‑44 (MQ=255)                                                                                                          gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCt                    >  1:284074/1‑139 (MQ=255)                                                                                                          gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                          <  2:832833/36‑1 (MQ=255)                                                                                                          gTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑cc                                                                                                                          >  1:832833/1‑36 (MQ=255)                                                                                                           tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTc                   >  2:107163/1‑139 (MQ=255)                                                                                                            aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCa                  <  1:151762/139‑1 (MQ=255)                                                                                                                ttttCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc              >  2:940496/1‑139 (MQ=255)                                                                                                                 tttCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc             <  2:686434/139‑1 (MQ=255)                                                                                                                 tttCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGcc             >  2:378496/1‑139 (MQ=255)                                                                                                                   tCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                   <  2:814919/115‑1 (MQ=255)                                                                                                                   tCCAGACTTCGCGCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtgg                                   >  1:814919/1‑115 (MQ=255)                                                                                                                            cgcgCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt  >  1:204892/1‑139 (MQ=255)                                                                                                                            cgcgCTGAAATGGCGTC‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAt  >  2:985975/1‑139 (MQ=255)                                                                                                                                            |                                                                                                                            TGGGATTCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGAT  >  NC_000913/1399840‑1400100

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