Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	180,550	T→G	P438P (CCT→CCG)	dgt →	deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	180,550	0	T	G	100.0%	179.6 / NA	49	P438P (CCT→CCG)	dgt	deoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/19);  total (30/19)

ATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCTATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCG  >  NC_000913/180402‑180695                                                                                                                                                     |                                                                                                                                                  aTGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGAt                                                                                                                                                 >  2:377260/1‑151 (MQ=255)  tgtgTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTTTCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGAtt                                                                                                                                                >  2:292200/1‑151 (MQ=255)   gtgtTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGAtt                                                                                                                                                >  2:183153/1‑150 (MQ=255)    tgtgTAGCCATCCAGATGTCGAGAGGCTTGAATTGCAGCGCTATTGGGTCATTAGCGCATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCATACTTTACTGAACTGGTAGAAAAAGCACGGGTGAAACGTTTTCCGATTg                                                                                                                                               <  2:132616/150‑1 (MQ=255)     gtTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGa                                                                                                                                              >  2:129852/1‑150 (MQ=255)             aTCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCtt                                                                                                                                     >  2:589630/1‑151 (MQ=255)               ccAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTAt                                                                                                                                   >  1:132792/1‑151 (MQ=255)                   aTGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTcc                                                                                                                                >  2:497622/1‑150 (MQ=255)                         aGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAAct                                                                                                                         >  1:858389/1‑151 (MQ=255)                                    ttGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGcc                                                                                                              >  1:346675/1‑151 (MQ=255)                                     tGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGcc                                                                                                              >  1:8634/1‑150 (MQ=255)                                      gCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCAt                                                                                                            <  2:527089/151‑1 (MQ=255)                                              aTCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTg                                                                                                      <  2:137410/149‑1 (MQ=255)                                                   gTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATg                                                                                               <  1:234577/151‑1 (MQ=255)                                                   gTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATg                                                                                               >  1:326961/1‑151 (MQ=255)                                                        tAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCg                                                                                            >  2:780636/1‑149 (MQ=255)                                                          gCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGc                                                                                        <  1:784524/151‑1 (MQ=255)                                                                 attaGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAg                                                                                  <  2:267183/150‑1 (MQ=255)                                                                    agagATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTaa                                                                               <  1:589645/150‑1 (MQ=255)                                                                      agaTTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAAtt                                                                            >  1:513544/1‑151 (MQ=255)                                                                         tttATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTAcc                                                                         >  1:787179/1‑151 (MQ=255)                                                                            aTCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTc                                                                      >  1:846988/1‑151 (MQ=255)                                                                              cGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAg                                                                    <  2:745421/151‑1 (MQ=255)                                                                               gTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTc                                                                      >  1:653759/1‑148 (MQ=255)                                                                                      attaAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCt                                                            >  1:712876/1‑151 (MQ=255)                                                                                        taAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTg                                                           >  2:129349/1‑150 (MQ=255)                                                                                                  gTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTa                                                <  2:405691/151‑1 (MQ=255)                                                                                                   ttATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTAt                                               >  2:740055/1‑151 (MQ=255)                                                                                                   ttAGCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTTCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCg                                                   <  2:342072/147‑1 (MQ=255)                                                                                                     aTCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATg                                              <  2:437024/150‑1 (MQ=255)                                                                                                      tCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATggg                                            >  2:389510/1‑151 (MQ=255)                                                                                                          aCTTTACTGAACTGGTAGAAAAAGAAAGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGaa                                          >  1:730891/1‑149 (MQ=255)                                                                                                              tACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAAtatta                                     <  1:399264/150‑1 (MQ=255)                                                                                                                  gAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATAttatta                                  <  2:591091/149‑1 (MQ=255)                                                                                                                  gAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTAcc                                >  1:68667/1‑151 (MQ=255)                                                                                                                    aCTGGTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGt                              >  2:564/1‑151 (MQ=255)                                                                                                                       ggTAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGtt                             <  1:90918/149‑1 (MQ=255)                                                                                                                         tAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTgccg                         >  1:815845/1‑151 (MQ=255)                                                                                                                         tAGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTgcc                          >  2:178186/1‑150 (MQ=255)                                                                                                                          aGAAAAAGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTgccgc                        >  2:52097/1‑151 (MQ=255)                                                                                                                             aaaaGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCctg                     <  2:612078/151‑1 (MQ=255)                                                                                                                             aaaaGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCctg                     >  1:490894/1‑151 (MQ=255)                                                                                                                                aGAACGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCctgctg                  >  2:519991/1‑151 (MQ=255)                                                                                                                                   aCGGGTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAg               <  2:178323/151‑1 (MQ=255)                                                                                                                                      ggTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGAt            <  1:60794/151‑1 (MQ=255)                                                                                                                                       gTGAAACGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGAtt           <  1:695271/151‑1 (MQ=255)                                                                                                                                          aaaCGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATtat         >  2:438083/1‑150 (MQ=255)                                                                                                                                             cGTTTCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGTCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCa     <  2:458792/151‑1 (MQ=255)                                                                                                                                                ttCCCGATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCg  >  2:851172/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                  ATGTGTTTAGCCATCCAGATGTCGAGCGGCTTGAATTGCAGGGCTATCGGGTCATTAGCGGATTATTAGAGATTTATCGTCCTTTATTAAGCCTGTCGTTATCAGACTTTACTGAACTGGTAGAAAAAGAACGGGTGAAACGTTTCCCTATTGAATCGCGCTTATTCCACAAACTCTCGACGCGCCATCGGCTGGCCTATGTCGAGGCTGTCAGTAAATTACCGTCAGATTCTCCTGAGTTTCCGCTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCG  >  NC_000913/180402‑180695

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