Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	179,442	A→T	H69L (CAC→CTC)	dgt →	deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	179,442	0	A	T	93.9%	207.4 / 1.9	66	H69L (CAC→CTC)	dgt	deoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base T (24/38);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 2.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGATGAA  >  NC_000913/179346‑179539                                                                                                 |                                                                                                  ggCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCtctc                                                                                                >  2:875347/1‑100 (MQ=255)      atcatcAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATgg                                                                                           <  1:1413976/100‑1 (MQ=255)       tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGa                                                                                          <  2:406404/100‑1 (MQ=255)        catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGaa                                                                                         <  2:1200374/100‑1 (MQ=255)        catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGaa                                                                                         >  1:323807/1‑100 (MQ=255)        catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGaa                                                                                         >  1:323727/1‑100 (MQ=255)        catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGaa                                                                                         >  2:1687660/1‑100 (MQ=255)         atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                        <  2:1392794/100‑1 (MQ=255)         atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                        >  1:680235/1‑100 (MQ=255)         atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                        <  1:1393376/100‑1 (MQ=255)         atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAg                                                                                        <  2:1600145/100‑1 (MQ=255)          tcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGt                                                                                       >  2:1767853/1‑100 (MQ=255)           caACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTc                                                                                      <  2:175698/100‑1 (MQ=255)             aCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCca                                                                                    >  1:1204971/1‑100 (MQ=255)                  cgggCAATAAGTCGTCAGCACCAAAAGAAGCAAGTGTTTCCACTGGAGCCGAATACCGCCGTGCGCACGCGTATTACCAACTCGAGGCTAGTCcagca                                                                                 >  1:1101221/3‑98 (MQ=255)                  ccGGTAATGCGTCTTCTGCAACAAAAGCCCCAGTTTGTTCCCCTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCcagc                                                                                  <  1:1710997/97‑1 (MQ=255)                            gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                     <  2:1204973/100‑1 (MQ=255)                            gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                     >  2:901675/1‑100 (MQ=255)                            gtcgtcTGCAACAAAAGACCCAGGTTTTTCCACCGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTa                                                                     <  1:1625644/100‑1 (MQ=255)                             tcgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTAc                                                                    <  2:614705/100‑1 (MQ=255)                              cgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTAc                                                                    <  1:41773/99‑1 (MQ=255)                              cgtcTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACa                                                                   <  1:493042/100‑1 (MQ=255)                                 cTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCg                                                                >  2:607768/1‑100 (MQ=255)                                  tGCAACAAAAGACCCGGGTTTTTCCACTGGCGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGc                                                               <  1:441617/100‑1 (MQ=255)                                   gCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGcc                                                              >  1:994074/1‑100 (MQ=255)                                   gCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGcc                                                              <  1:1867876/100‑1 (MQ=255)                                    caacaaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCa                                                             >  2:1059773/1‑100 (MQ=255)                                    caacaaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCa                                                             >  1:1046837/1‑100 (MQ=255)                                    caacaaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCACCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCa                                                             >  1:459887/1‑100 (MQ=255)                                       caaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                                          >  1:495376/1‑100 (MQ=255)                                       caaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                                          <  2:751543/100‑1 (MQ=255)                                       caaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                                          >  2:257097/1‑100 (MQ=255)                                       caaAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAg                                                          >  2:992767/1‑100 (MQ=255)                                        aaagGACCCAGGTGTTTCCCCTGGAGCGCAATGCCGCCGTGCGCACGCGACTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGa                                                         <  1:750503/100‑1 (MQ=255)                                        aaaaGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGa                                                         <  2:1179461/100‑1 (MQ=255)                                        aaaaGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGa                                                         <  2:1690373/100‑1 (MQ=255)                                         aaaGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGaa                                                        >  1:936990/1‑100 (MQ=255)                                                 aGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAg                                                >  1:1637181/1‑100 (MQ=255)                                                   gtgtttCCACTGGAGCGCAATGCCGCCGTGCTCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGc                                               <  1:1366043/96‑1 (MQ=255)                                                    tttttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCg                                             <  2:1693097/100‑1 (MQ=255)                                                      tttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTc                                           <  2:936991/100‑1 (MQ=255)                                                      tttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTc                                           <  2:1031875/100‑1 (MQ=255)                                                       ttCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCt                                          >  2:938476/1‑100 (MQ=255)                                                        tCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTc                                           <  1:860744/98‑1 (MQ=255)                                                        tCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTg                                         >  2:141117/1‑100 (MQ=255)                                                         ccACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGa                                        <  2:96637/100‑1 (MQ=255)                                                         ccACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGa                                        <  2:680236/100‑1 (MQ=255)                                                          cACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGa                                        <  1:403602/99‑1 (MQ=255)                                                            cTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAag                                     >  1:136791/1‑100 (MQ=255)                                                              ggAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                                   <  2:1288186/100‑1 (MQ=255)                                                              ggAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAagag                                   >  1:1804859/1‑100 (MQ=255)                                                               gAGCGCAATGGCGCCGTCCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGaa                                       <  1:310107/95‑1 (MQ=255)                                                                 gcgcAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCtt                                <  1:1450273/100‑1 (MQ=255)                                                                  cgcAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTa                               <  2:896513/100‑1 (MQ=255)                                                                   gcAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTaa                              >  1:924985/1‑100 (MQ=255)                                                                   gcAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTaa                              <  1:1634324/100‑1 (MQ=255)                                                                    cAATGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTaaa                             >  2:1671472/1‑100 (MQ=255)                                                                      aTGCCGCCGTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGTAAGAGCTTAAAtt                           <  2:1320209/100‑1 (MQ=255)                                                                          cgccgTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTg                       <  1:166634/100‑1 (MQ=255)                                                                             cgTGCGCACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGaa                    <  2:1637191/100‑1 (MQ=255)                                                                                 cgcACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCAt                <  2:247355/100‑1 (MQ=255)                                                                                 cgcACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCAt                <  2:247475/100‑1 (MQ=255)                                                                                   cACGCGTCTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATAc              >  2:83959/1‑100 (MQ=255)                                                                                          cTTACCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTgg       >  2:609538/1‑100 (MQ=255)                                                                                             aCCCTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGATg    <  1:271461/100‑1 (MQ=255)                                                                                              cccTCTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGATGa   <  1:295689/100‑1 (MQ=255)                                                                                               ccACTCGATGGAAATCCAGCAGGTGGGGCGCTACATCGACAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGATGaa  >  1:1203376/1‑100 (MQ=255)                                                                                                 |                                                                                                  GGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGCACGCGTCTTACCCACTCGATGGAAGTCCAGCAGGTGGGGCGCTACATCGCCAAAGAAATTTTAAGCCGTCTGAAAGAGCTTAAATTACTGGAAGCATACGGCCTGGATGAA  >  NC_000913/179346‑179539

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