Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,841,897 G→A *136* (TAG→TAA
T80I (ACT→ATT) 
nudG →
ynjH ←
CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
DUF1496 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,841,8970GA100.0% 137.3 / NA 43*136* (TAG→TAA
T80I (ACT→ATT) 
nudG
ynjH
CTP pyrophosphohydrolase, also hydrolyzes 2‑hydroxy‑dATP, 8‑hydroxy‑dGTP, 5‑hydroxy‑CTP, dCTP and 5‑methyl‑dCTP
DUF1496 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/21);  total (22/21)

GTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG  >  NC_000913/1841765‑1842018
                                                                                                                                    |                                                                                                                         
gTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTaa                                                                                                                     >  2:337282/1‑139 (MQ=255)
 ttACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAg                                                                                                                    >  2:676127/1‑139 (MQ=255)
       gCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGtttt                                                                                                              >  1:854519/1‑139 (MQ=255)
            tGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAgg                                                                                                                   <  1:607513/129‑1 (MQ=255)
            tGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAgg                                                                                                                   >  2:607513/1‑129 (MQ=255)
                             tctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcagc                                                                                        <  2:232173/139‑1 (MQ=255)
                                      cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 >  1:498131/1‑121 (MQ=255)
                                      cACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 <  2:498131/121‑1 (MQ=255)
                                                 gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 <  1:308589/110‑1 (MQ=255)
                                                 gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 >  2:308589/1‑110 (MQ=255)
                                                 gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         <  1:1137583/118‑1 (MQ=255)
                                                 gcgcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         >  2:1137583/1‑118 (MQ=255)
                                                                 tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCCTTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaa                                                    >  2:502367/1‑139 (MQ=255)
                                                                 tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                           >  1:222574/1‑100 (MQ=255)
                                                                 tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                           <  2:222574/100‑1 (MQ=255)
                                                                 tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt                                                         <  1:499948/134‑1 (MQ=255)
                                                                 tGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCtt                                                         >  2:499948/1‑134 (MQ=255)
                                                                    ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         <  1:525781/99‑1 (MQ=255)
                                                                    ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         >  2:525781/1‑99 (MQ=255)
                                                                    ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGcagcag                                                                                      <  1:699156/102‑1 (MQ=255)
                                                                    ccccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGcagcag                                                                                      >  2:699156/1‑102 (MQ=255)
                                                                    ccccTACTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcagc                                                                                        >  1:101869/1‑100 (MQ=255)
                                                                    ccccTACTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcagc                                                                                        <  2:101869/100‑1 (MQ=255)
                                                                      ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                      <  2:1271822/68‑1 (MQ=255)
                                                                      ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                      >  1:1271822/1‑68 (MQ=255)
                                                                      ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         >  2:56890/1‑97 (MQ=255)
                                                                      ccTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         <  1:56890/97‑1 (MQ=255)
                                                                            gACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                         <  1:676127/139‑1 (MQ=255)
                                                                                           gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                             >  2:170888/1‑136 (MQ=255)
                                                                                           gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACgca                             <  1:170888/136‑1 (MQ=255)
                                                                                             ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              >  1:78486/1‑69 (MQ=255)
                                                                                             ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              <  2:78486/69‑1 (MQ=255)
                                                                                                 tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                              <  2:1096981/81‑1 (MQ=255)
                                                                                                 tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                              >  1:1096981/1‑81 (MQ=255)
                                                                                                  ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTa                   >  1:953689/1‑139 (MQ=255)
                                                                                                          ttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTg           <  1:502367/139‑1 (MQ=255)
                                                                                                           tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                         <  1:937303/108‑1 (MQ=255)
                                                                                                           tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                         >  2:937303/1‑108 (MQ=255)
                                                                                                            aCGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            <  2:277203/56‑1 (MQ=255)
                                                                                                            aCGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            >  1:277203/1‑56 (MQ=255)
                                                                                                                   gccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg  <  1:1299589/139‑1 (MQ=255)
                                                                                                                            gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            >  2:308046/1‑40 (MQ=38)
                                                                                                                            gCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            <  1:308046/40‑1 (MQ=255)
                                                                                                                                    |                                                                                                                         
GTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG  >  NC_000913/1841765‑1842018

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: