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,8970GA98.0% 153.0 / ‑6.3 49*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);  major base A (24/24);  minor base C (0/1);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG  >  NC_000913/1841770‑1842018
                                                                                                                               |                                                                                                                         
aGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGtt                                                                                                                <  2:511034/139‑1 (MQ=255)
     caTGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCGTTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtc                                                                                                           >  2:1188462/1‑139 (MQ=255)
     caTGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGtc                                                                                                           <  1:835162/139‑1 (MQ=255)
         aaCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCCGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGTCAGACCAGCGGATTAATGCTAAGGGTTTTGtcatca                                                                                                       <  1:836781/139‑1 (MQ=255)
                ggCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGc                                                                                                >  1:310337/1‑139 (MQ=255)
                  cgcTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              >  1:298250/1‑139 (MQ=255)
                    ctggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            <  2:999135/139‑1 (MQ=255)
                     tggtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgc                                                                                           <  2:184971/139‑1 (MQ=255)
                       gtctggtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         <  2:310337/139‑1 (MQ=255)
                            gtGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTa                                                                                    >  2:734531/1‑139 (MQ=255)
                                  aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGTATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                              <  2:779284/139‑1 (MQ=255)
                                  aCCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCtt                                                                              <  2:298250/139‑1 (MQ=255)
                                      gAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGc                                                                          >  1:173104/1‑139 (MQ=255)
                                              gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCAc                                                                                                      >  2:546303/1‑103 (MQ=255)
                                              gcTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCAc                                                                                                      <  1:546303/103‑1 (MQ=255)
                                                 gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                      <  1:652725/84‑1 (MQ=255)
                                                 gCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTa                                                                                                                      >  2:652725/1‑84 (MQ=255)
                                                           cTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGa                                                     <  2:55215/139‑1 (MQ=255)
                                                              gCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGaata                                                  >  1:634205/1‑139 (MQ=255)
                                                                  ctgctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAAttt                                              >  1:1160300/1‑139 (MQ=255)
                                                                     ctgACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGa                                           <  2:173104/139‑1 (MQ=255)
                                                                       gACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATc                                         <  1:734531/139‑1 (MQ=255)
                                                                             ccATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGAGTTTGTCATCACGCTGGCCTTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATa                                   >  2:1158807/1‑139 (MQ=255)
                                                                                    taGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                              <  1:1045457/105‑1 (MQ=255)
                                                                                    taGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                              >  2:1045457/1‑105 (MQ=255)
                                                                                      gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         >  2:1021225/1‑76 (MQ=255)
                                                                                      gagGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTgcag                                                                                         <  1:1021225/76‑1 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTACTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACTGCCCCTTCAGAATAATTTTGATCCTCATAAACGCAGCACt                        <  1:799758/139‑1 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAgg                                                                                                                   <  1:544886/48‑1 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAgg                                                                                                                   >  2:544886/1‑48 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCAc                                                                                                      <  2:1021022/61‑1 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCAc                                                                                                      >  1:1021022/1‑61 (MQ=255)
                                                                                        ggCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAAGGCAGCACt                        >  2:799758/1‑139 (MQ=255)
                                                                                         gCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              <  2:990834/68‑1 (MQ=255)
                                                                                         gCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              >  1:990834/1‑68 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              <  1:606305/65‑1 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCAt                                                                                              >  2:606305/1‑65 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            <  2:602161/67‑1 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            >  1:602161/1‑67 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                              >  1:390691/1‑97 (MQ=255)
                                                                                            tttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGc                                                              <  2:390691/97‑1 (MQ=255)
                                                                                             ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 <  2:956218/61‑1 (MQ=255)
                                                                                             ttATGGCTTTACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTgg                                                                                                 >  1:956218/1‑61 (MQ=255)
                                                                                                    tttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            >  1:231655/1‑59 (MQ=255)
                                                                                                    tttACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTg                                                                                            <  2:231655/59‑1 (MQ=255)
                                                                                                      tACGCGCCGCCAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTACTTAGGCTTTAATTAGCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTAAAGGGCTgc          >  2:801697/1‑139 (MQ=255)
                                                                                                             cgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAAtt                                                                   >  2:210615/1‑75 (MQ=255)
                                                                                                             cgccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAAtt                                                                   <  1:210615/75‑1 (MQ=255)
                                                                                                              gccAGACCAGCGGATTAATGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTg  >  1:1207576/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                         
AGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCATCACGCTGGCATTGCAGCAGTATTCCTTCGGCTTTAATTACCGCCCCTTCAGAATAATTTTGATCCTGATAAACGCAGCACTGAGTACAGGGCTGCGCTGACTG  >  NC_000913/1841770‑1842018

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