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%	129.9 / NA	39	T87A (ACG→GCG)	ogt	O‑6‑alkylguanine‑DNA:cysteine‑protein methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (17/22);  total (17/22)

TCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399846‑1400111                                                                                                                                      |                                                                                                                                      tCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccccc                                                                                                                                 <  2:756720/139‑1 (MQ=255)    aTCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccg                                                                                                                              <  2:1516228/139‑3 (MQ=255)               gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  2:1123581/139‑1 (MQ=255)               gCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    >  1:28999/1‑139 (MQ=255)                 aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  <  2:116931/139‑1 (MQ=255)                 aCCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  <  2:1146535/139‑1 (MQ=255)                   cAACGGCACGCGCCGCGCCAGGCCGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                                     >  2:1193156/1‑117 (MQ=255)                   cAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cc                                                                                                                                     <  1:1193156/117‑1 (MQ=255)                           cgcgccgcgcCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaata                                                                                                        >  1:415046/1‑139 (MQ=255)                            gcgccgcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaataa                                                                                                       <  2:415046/139‑1 (MQ=255)                                 gcgccAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTa                                                                                                  >  1:510073/1‑139 (MQ=255)                                             cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCa                                                                                      >  2:960626/1‑139 (MQ=255)                                             cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccgta                                                                                                                            <  2:77203/100‑5 (MQ=255)                                             cccAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑ccccccccgta                                                                                                                            >  1:77203/1‑96 (MQ=255)                                                              ttGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCGAAATATTCACGAAGCtt                                                                     <  1:960626/139‑1 (MQ=255)                                                               tGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTg                                                                    <  1:1336287/139‑1 (MQ=255)                                                                        tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  >  2:135095/1‑84 (MQ=255)                                                                        tGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  <  1:135095/84‑1 (MQ=255)                                                                           aTTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaata                                                                                                        <  2:1150544/91‑1 (MQ=255)                                                                           aTTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaata                                                                                                        >  1:1150544/1‑91 (MQ=255)                                                                             tACCTGCCCGCTGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGcc                                                      <  1:994906/139‑1 (MQ=255)                                                                               ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           >  1:1024611/1‑68 (MQ=255)                                                                               ccTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           <  2:1024611/68‑1 (MQ=255)                                                                                   cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    >  2:416928/1‑71 (MQ=255)                                                                                   cccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  1:416928/71‑1 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           >  2:612425/1‑63 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGc                                                                                                                           <  1:612425/63‑1 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                          >  2:39593/1‑80 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCaa                                                                                                          <  1:39593/80‑1 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   >  2:276702/1‑54 (MQ=255)                                                                                    ccGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   <  1:276702/54‑1 (MQ=255)                                                                                      gCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATtggt                                             <  2:1259448/139‑1 (MQ=255)                                                                                                  gCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCAc                                                                            <  1:657168/96‑1 (MQ=255)                                                                                                  gCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCAc                                                                            >  2:657168/1‑96 (MQ=255)                                                                                                   cGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    >  1:325829/1‑55 (MQ=255)                                                                                                   cGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGa                                                                                                                    <  2:325829/55‑1 (MQ=255)                                                                                                   acaggttttttCGAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATACACGAGGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATgcgc                                <  1:915110/132‑1 (MQ=255)                                                                                                     tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTc                              >  1:828193/1‑139 (MQ=255)                                                                                                     tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   <  2:1097161/37‑1 (MQ=255)                                                                                                     tAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   >  1:1097161/1‑37 (MQ=255)                                                                                                      aGTGTTTTCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCa                             >  2:204224/1‑139 (MQ=255)                                                                                                      aGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   >  2:1256726/1‑36 (MQ=255)                                                                                                      aGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑cccc                                                                                                                                   <  1:1256726/36‑1 (MQ=255)                                                                                                          ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  <  1:198657/50‑1 (MQ=255)                                                                                                          ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAg                                                                                                                  >  2:198657/1‑50 (MQ=255)                                                                                                          ttttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGc                         <  2:828193/139‑1 (MQ=255)                                                                                                           tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATAtt                                                                               >  1:1485950/1‑84 (MQ=255)                                                                                                           tttCCAGACTTCGCGCTGAAATGGCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATAtt                                                                               <  2:1485950/84‑1 (MQ=255)                                                                                                                  aCTTCGCGCTGAAATGGCGTCCCCCCCCATACCTGTGGGAAGCTTAGCAATAATGCTAATATTACCGGCAAAATATTCACGATGCTTGTCGCTTAAGCCACCTGGATTGGTGGCAGAAATGCGCTCATAGCctacctgt                 >  2:915110/1‑133 (MQ=255)                                                                                                                                 ggCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  2:28999/139‑1 (MQ=255)                                                                                                                                 ggCGTCCCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTcc  <  2:510073/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                      TCGATCCGTTTGCCGCACCAACGGCACGCGCCGCGCCAGGACGGCCCAATTGCTCAGCCAGTTGGCCGTAATGCATTACCTGCCCGCAGGGGATAGTGCGTAGTGTTTTCCAGACTTCGCGCTGAAATGGCGT‑CCCCCCCGTAGCAGTGGGAAGCGTATCAATAATGCTAAGATTACCGGCAAAATATTCACGAAGCTTGTCGCTTAAACCGCCTGGATTGGTGGCAGAAATGCGCTCATAGCCTTCTTTGCGATAATGGATGTCC  >  NC_000913/1399846‑1400111

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