Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	111,433	G→A	*130* (TAG→TAA)	mutT →	dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	111,433	0	G	A	100.0%	180.1 / NA	56	*130* (TAG→TAA)	mutT	dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/29);  total (27/29)

CCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAC  >  NC_000913/111323‑111567                                                                                                               |                                                                                                                                       ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCg                                                                                                            >  1:1349985/1‑139 (MQ=255)  cGTGGGGTAAAGAAGGGCAACCCGGTGATTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCATCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGa                                                                                                           <  2:907027/139‑1 (MQ=255)             gaagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                <  2:269704/139‑1 (MQ=255)             gaagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                <  2:1349985/139‑1 (MQ=255)                 ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATg                                                                                            >  2:583613/1‑139 (MQ=255)                 ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATg                                                                                            >  1:738242/1‑139 (MQ=255)                  gCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGa                                                                                                           <  2:98327/123‑1 (MQ=255)                  gCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGa                                                                                                           >  1:98327/1‑123 (MQ=255)                   cAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGcc                                                                                          >  2:308289/1‑139 (MQ=255)                     aCCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTg                                                                                        >  2:73659/1‑139 (MQ=255)                      cccGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTGCGCCAGGCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                                               >  1:1038702/1‑131 (MQ=255)                      cccGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                                               <  2:1038702/131‑1 (MQ=255)                        cGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATg                                                                                     <  1:583613/139‑1 (MQ=255)                                 gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  2:1358752/139‑1 (MQ=255)                                         tGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGTACACg                                                                                               >  2:1076656/1‑112 (MQ=255)                                         tGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGTACACg                                                                                               <  1:1076656/112‑1 (MQ=255)                                           gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  <  1:769294/139‑1 (MQ=255)                                            tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                 <  1:1631057/139‑1 (MQ=255)                                            tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                 <  1:1786497/139‑1 (MQ=255)                                            tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                 >  1:609080/1‑139 (MQ=255)                                            tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                                                 >  1:26604/1‑139 (MQ=255)                                                tcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                             <  2:609080/139‑1 (MQ=255)                                                tcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                                             >  2:626432/1‑139 (MQ=255)                                                   tAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                          <  1:308289/139‑1 (MQ=255)                                                        ccGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                                                     <  1:1240390/139‑1 (MQ=255)                                                          gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               >  2:948462/1‑63 (MQ=255)                                                          gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               <  1:948462/63‑1 (MQ=255)                                                          gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAt                                                   >  1:447224/1‑139 (MQ=255)                                                          gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAt                                                   >  1:1402276/1‑139 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                                               <  2:442152/93‑1 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                                               >  1:442152/1‑93 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            >  2:1009384/1‑112 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  1:832935/112‑1 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            >  2:832935/1‑112 (MQ=255)                                                            tgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  1:1009384/112‑1 (MQ=255)                                                              atTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGTACACGATGCCTGATGCGACGCTGgcgc                                                                         <  1:783471/113‑1 (MQ=255)                                                              atTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGTACACGATGCCTGATGCGACGCTGgcgc                                                                         >  2:783471/1‑113 (MQ=255)                                                                         gccaATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGAt                                    >  1:1523842/1‑139 (MQ=255)                                                                                 aCCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCa                                         >  2:1188494/1‑126 (MQ=255)                                                                                 aCCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCa                                         <  1:1188494/126‑1 (MQ=255)                                                                                   cGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                                       <  2:1685156/94‑1 (MQ=255)                                                                                   cGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                                       >  1:1685156/1‑94 (MQ=255)                                                                                     gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAgg                                                      >  2:560502/1‑109 (MQ=255)                                                                                     gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAgg                                                      <  1:560502/109‑1 (MQ=255)                                                                                     gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGttttt                        <  2:145668/139‑1 (MQ=255)                                                                                         ttGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAg                    <  2:26604/139‑1 (MQ=255)                                                                                          tGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAgg                   <  1:626432/139‑1 (MQ=255)                                                                                            cGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                                                     <  2:675483/103‑1 (MQ=255)                                                                                            cGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                                                     >  1:675483/1‑103 (MQ=255)                                                                                             gAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  2:385857/79‑1 (MQ=255)                                                                                             gAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            >  1:385857/1‑79 (MQ=255)                                                                                                           tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGcaca                                                                                                 >  2:1590353/1‑44 (MQ=255)                                                                                                           tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGcaca                                                                                                 <  1:1590353/44‑1 (MQ=255)                                                                                                           tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGAc                                                                                 >  2:189824/1‑60 (MQ=255)                                                                                                           tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGAc                                                                                 <  1:189824/60‑1 (MQ=255)                                                                                                           tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAc  >  1:829928/1‑139 (MQ=255)                                                                                                               |                                                                                                                                       CCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAC  >  NC_000913/111323‑111567

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