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_000913111,4330GA100.0% 178.0 / NA 55*130* (TAG→TAA) mutTdGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/26);  total (29/26)

GTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACG  >  NC_000913/111302‑111568
                                                                                                                                    |                                                                                                                                       
gTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGAt                                                                                                                                  >  2:19739/1‑139 (MQ=255)
               ggggAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGTTTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccg                                                                                                                   <  1:19739/139‑1 (MQ=255)
                   aGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATc                                                                                                               >  2:1495940/1‑139 (MQ=255)
                     ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccgc                                                                                                                  >  1:732706/1‑134 (MQ=255)
                     ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccgc                                                                                                                  <  2:732706/134‑1 (MQ=255)
                     ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                     <  1:697233/131‑1 (MQ=255)
                     ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                     >  2:697233/1‑131 (MQ=255)
                     ccGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCg                                                                                                             >  2:340183/1‑139 (MQ=255)
                        tGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACa                                                                                                          <  1:1495940/139‑1 (MQ=255)
                               aagaagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGcac                                                                                                   >  1:1069757/1‑139 (MQ=255)
                                 gaagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 <  2:1455551/139‑1 (MQ=255)
                                 gaagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 <  1:1137939/139‑1 (MQ=255)
                                  aagGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACg                                                                                                >  1:556066/1‑139 (MQ=255)
                                    gggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGAt                                                                                              >  2:418754/1‑139 (MQ=255)
                                      gCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGc                                                                                            >  2:1308025/1‑139 (MQ=255)
                                       cAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGcc                                                                                           <  1:1007102/139‑1 (MQ=255)
                                         aCCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 <  2:385468/131‑1 (MQ=255)
                                         aCCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 >  1:385468/1‑131 (MQ=255)
                                          cccGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCATAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCACACACGATGCCTGa                                                                                        >  2:1079333/1‑139 (MQ=255)
                                                     gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                             <  1:418754/139‑1 (MQ=255)
                                                     gATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                             <  2:1069757/139‑1 (MQ=255)
                                                        gTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCccac                                                                                                 <  1:1024410/117‑1 (MQ=39)
                                                        gTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 >  2:1024410/1‑117 (MQ=255)
                                                            cTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCt                                                                      >  2:488839/1‑139 (MQ=255)
                                                               gtcggtcTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 >  1:400944/1‑109 (MQ=255)
                                                               gtcggtcTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                 <  2:400944/109‑1 (MQ=255)
                                                               gtcggtcTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   <  1:340183/139‑1 (MQ=255)
                                                               gtcggtcTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCACACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   <  1:1079333/139‑1 (MQ=255)
                                                                              gatgatTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGcgcg                                                                         <  1:319845/118‑1 (MQ=255)
                                                                              gatgatTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGcgcg                                                                         >  2:319845/1‑118 (MQ=255)
                                                                              gatgatTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAt                                                    >  1:17829/1‑139 (MQ=255)
                                                                              gatgatTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAt                                                    >  2:352859/1‑139 (MQ=255)
                                                                                tgatTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAttt                                                  >  1:614583/1‑139 (MQ=255)
                                                                                         gccag‑ccaATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCAt                                         <  1:352859/139‑1 (MQ=255)
                                                                                            ag‑ccaATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGa                                      <  2:556066/139‑1 (MQ=255)
                                                                                             g‑ccaATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGAt                                     >  2:776270/1‑139 (MQ=255)
                                                                                                  aTGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACg                                                                                 >  2:1195496/1‑91 (MQ=255)
                                                                                                  aTGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACg                                                                                 <  1:1195496/91‑1 (MQ=255)
                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGtttt                          >  2:338793/1‑139 (MQ=255)
                                                                                                         ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAA‑GGATTTCTAACTCATTGATAAATTTGttttt                         >  1:99305/1‑139 (MQ=255)
                                                                                                          gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   <  1:1100261/97‑1 (MQ=255)
                                                                                                          gTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   >  2:1100261/1‑97 (MQ=255)
                                                                                                           tAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTg                        <  1:1308025/139‑1 (MQ=255)
                                                                                                                gCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   >  2:1301575/1‑91 (MQ=255)
                                                                                                                gCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                   <  1:1301575/91‑1 (MQ=255)
                                                                                                                  gAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGt                       >  1:665908/1‑133 (MQ=255)
                                                                                                                  gAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGt                       <  2:665908/133‑1 (MQ=255)
                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccgc                                                                                                                  >  2:358304/1‑41 (MQ=255)
                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccgc                                                                                                                  <  1:358304/41‑1 (MQ=255)
                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGTGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGt                   >  1:143287/1‑136 (MQ=255)
                                                                                                                   aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGTGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGt                   <  2:143287/136‑1 (MQ=255)
                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGgc                                                                            <  2:1088977/66‑1 (MQ=255)
                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGgc                                                                            >  1:1088977/1‑66 (MQ=255)
                                                                                                                                tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAc   >  1:607170/1‑139 (MQ=255)
                                                                                                                                 gTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAcg  <  2:17829/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAG‑CCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACG  >  NC_000913/111302‑111568

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