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%	189.0 / NA	58	*130* (TAG→TAA)	mutT	dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/31);  total (27/31)

GTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTT  >  NC_000913/111302‑111545                                                                                                                                    |                                                                                                                 gTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGAt                                                                                                           >  1:1274300/1‑139 (MQ=255)           gggAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGtttt                                                                                                 >  2:393895/1‑139 (MQ=255)             gAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcg                                                                                               <  1:821563/139‑1 (MQ=255)                    aGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATc                                                                                        >  1:1401037/1‑139 (MQ=255)                      ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCg                                                                                      >  1:129708/1‑139 (MQ=255)                               aaaGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGca                                                                             <  2:824619/139‑1 (MQ=255)                                   aagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                         <  1:388495/139‑1 (MQ=255)                                      ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATg                                                                      >  2:489862/1‑139 (MQ=255)                                              ggTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGc                                                              >  1:1026785/1‑139 (MQ=255)                                                 gAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccg                                                                                            <  1:1187553/106‑1 (MQ=255)                                                 gAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccg                                                                                            >  2:1187553/1‑106 (MQ=255)                                                  aGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACg                                                          <  2:1401037/139‑1 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      <  2:108344/139‑1 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      >  1:1094160/1‑139 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      <  1:359635/139‑1 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      >  1:108344/1‑139 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      <  2:129708/139‑1 (MQ=255)                                                      gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                      <  1:489862/139‑1 (MQ=255)                                                         gTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGgcgc                                                   >  2:559834/1‑139 (MQ=255)                                                          tCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTCAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGcgcg                                                  <  2:1026785/139‑1 (MQ=255)                                                            gCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTc                                                >  1:52325/1‑139 (MQ=255)                                                             cTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCt                                               >  2:1351848/1‑139 (MQ=255)                                                                gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                            <  1:1351848/139‑1 (MQ=255)                                                                 tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                           >  1:1136097/1‑139 (MQ=255)                                                                     tcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGc                                       <  2:1094160/139‑1 (MQ=255)                                                                        tAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                    <  2:1274300/139‑1 (MQ=255)                                                                           tGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                           >  1:948634/1‑129 (MQ=255)                                                                           tGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATc                                           <  2:948634/129‑1 (MQ=255)                                                                           tGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAg                                 <  2:1136097/139‑1 (MQ=255)                                                                           tGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAg                                 <  1:559834/139‑1 (MQ=255)                                                                            gCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAgg                                <  2:1116924/139‑1 (MQ=255)                                                                             ccGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                               <  1:393895/139‑1 (MQ=255)                                                                               gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                               <  1:1394295/137‑1 (MQ=255)                                                                               gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                               >  2:1394295/1‑137 (MQ=255)                                                                               gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTT‑AACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                 >  1:765036/1‑118 (MQ=255)                                                                               gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTT‑AACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                 <  2:765036/118‑1 (MQ=255)                                                                                atgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                          <  2:310928/93‑1 (MQ=255)                                                                                atgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                          >  1:310928/1‑93 (MQ=255)                                                                                                     aaCCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTg       <  2:52325/139‑1 (MQ=255)                                                                                                        cGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                         >  2:1429224/1‑70 (MQ=255)                                                                                                        cGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                         <  1:1429224/70‑1 (MQ=255)                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                          <  2:715819/68‑1 (MQ=255)                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                          >  1:715819/1‑68 (MQ=255)                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAg                                 >  2:594450/1‑109 (MQ=255)                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAg                                 <  1:594450/109‑1 (MQ=255)                                                                                                          gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                            >  2:1380660/1‑97 (MQ=255)                                                                                                          gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                            <  1:1380660/97‑1 (MQ=255)                                                                                                          gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGttttt  >  2:142442/1‑139 (MQ=255)                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGa                                                                                     >  1:679201/1‑47 (MQ=255)                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGa                                                                                     <  2:679201/47‑1 (MQ=255)                                                                                                                     gCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaaa                                  >  2:590251/1‑96 (MQ=255)                                                                                                                     gCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaaa                                  <  1:590251/96‑1 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGc                                                                              <  1:23817/41‑1 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGc                                                                              >  2:23817/1‑41 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCg                                                             >  2:804941/1‑58 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCg                                                             <  1:804941/58‑1 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTc                    <  2:1027979/99‑1 (MQ=255)                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTc                    >  1:1027979/1‑99 (MQ=255)                                                                                                                                    |                                                                                                                 GTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTT  >  NC_000913/111302‑111545

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