Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,141,832 A→T D128E (GAT→GAA)  dcd ← deoxycytidine triphosphate deaminase; dCTP deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,141,8320AT93.6% 151.5 / NA 47D128E (GAT→GAA) dcddeoxycytidine triphosphate deaminase; dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (23/21);  minor base G (3/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAA  >  NC_000913/2141686‑2141975
                                                                                                                                                  |                                                                                                                                               
gCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGGTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGGTCg                                                                                                                                              >  2:727168/1‑150 (MQ=255)
 cATCTTCACGGCGGTTGTAAAGTCGCACCGCCCGGCCCGGAAGCGGCCCCAAGCCCAGCGCACCCATTAACATCCCCGGACGCACCGCCCGAGGCAACTTTCCGGAGTTGTTGAAATCCCACCCAAAGCAACCCGAACCGCCCCGGTCggt                                                                                                                                            >  1:296052/1‑149 (MQ=255)
   tCTTCACGGCGGTTGTAAGGGCGCACCGCCGGGCCGGAAAACGGCTCAAAACTCAGCGCACCAATTAAAATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGGTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATg                                                                                                                                           >  2:971107/1‑150 (MQ=255)
       cACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTg                                                                                                                                      <  2:311439/151‑1 (MQ=255)
                  tAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGt                                                                                                                           >  2:644686/1‑151 (MQ=255)
                  tAAGGACTCACCGCCGGGCCGGAATGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCAGTTCGATGCGGTGCGCGGTGAc                                                                                                                             <  1:867395/149‑1 (MQ=255)
                   aaGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTg                                                                                                                          <  2:610267/151‑1 (MQ=255)
                                cgggccggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGcccc                                                                                                             >  1:979939/1‑151 (MQ=255)
                                   gccggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGa                                                                                                          <  2:510241/151‑1 (MQ=255)
                                                   aaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAAc                                                                                          <  2:1070331/151‑1 (MQ=255)
                                                     gCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGc                                                                                        >  2:780649/1‑151 (MQ=255)
                                                         aGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGt                                                                                    >  1:885244/1‑151 (MQ=255)
                                                           cgcACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcc                                                                                  <  2:66244/151‑1 (MQ=255)
                                                           cgcACCAATTAACAGGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGAAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCGCCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGt                                                                                    <  1:849308/149‑1 (MQ=255)
                                                            gcACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGAACAATGCAACCAGACCAGCCCGGTTCGACGCGGCGCGCGGCGACGAGCACCACCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcc                                                                                  >  1:587371/1‑150 (MQ=255)
                                                               ccAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagcc                                                                              >  2:909966/1‑151 (MQ=255)
                                                               ccAATTAACATGCCCGGACGCACCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGGTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGGGAGGAACGCGCGTccagcc                                                                              >  1:399360/1‑151 (MQ=255)
                                                                           cccGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCccacc                                                                      >  2:346628/1‑147 (MQ=255)
                                                                           cccGGACGCAGCGCCACCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGAt                                                                  <  1:368786/151‑1 (MQ=255)
                                                                                aCGCATAGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCCGACCAGCCCGGTTCGATGCGGTGCGCGGAGACTTGCACCATCAGCCCCAGACGCCCCAGTGAGGAACGCCCGTCCAGCCAGCCACCCAGATCGGc                                                              >  1:358861/1‑150 (MQ=255)
                                                                                   cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTgg                                                           <  1:127657/150‑1 (MQ=255)
                                                                                   cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGc                                                          >  2:450502/1‑151 (MQ=255)
                                                                                            cGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCAcc                                                 <  2:320346/151‑1 (MQ=255)
                                                                                                aGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACt                                             >  2:620387/1‑151 (MQ=255)
                                                                                                   ttACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTcc                                           >  1:662340/1‑150 (MQ=255)
                                                                                                    tACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCTCACCAGATCGGCTGGCAGCGTCACCGACTCCa                                          <  2:1047807/150‑1 (MQ=255)
                                                                                                    tACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTcc                                           >  1:835562/1‑149 (MQ=255)
                                                                                                     aCCGGAGTTGTAGACCTCCAGCACCATGCAACCAGACCAGCCCGGTCCGTTGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTcc                                           <  1:851300/148‑1 (MQ=255)
                                                                                                     aCCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAg                                         <  2:669258/150‑1 (MQ=255)
                                                                                                     aCCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGc                                        <  2:923487/151‑1 (MQ=255)
                                                                                                             tGTAGACCTCCAGCACAATGCAACCAGACCAGCCCGGTGCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGGGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcg                                 <  1:608340/150‑1 (MQ=255)
                                                                                                             tGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcg                                 <  2:341399/150‑1 (MQ=255)
                                                                                                                   aCTCCAGCACAATGCAACCAGACAAGCCCGGTTCGATGCGGTGCGCGGAGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCAACAGCTCGGCTGACAGCGTCACCGACTCCCGCGTCACcgccagcgc                          >  1:477700/1‑151 (MQ=255)
                                                                                                                    cTCCAGCACAATGCAACCAGACCAGCCCGGTTCGCTGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGAGTCAccgcccgc                            >  1:778047/1‑148 (MQ=255)
                                                                                                                    cTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcc                         >  2:719509/1‑151 (MQ=255)
                                                                                                                      ccAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcc                         <  2:349892/149‑1 (MQ=255)
                                                                                                                        aGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCgccagcgccagc                      >  2:949806/1‑150 (MQ=255)
                                                                                                                          cacaATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctc                    >  2:725508/1‑150 (MQ=255)
                                                                                                                           acaATGTAACCAGACTAGCCCGGTTCGATGTGGTGCGAGGTGACGTGCACCATCAGCCCTAGACGCGCCAGTGAGGAAGGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTGACCGACTATAGAGTCACCGCGAGCGCCAGCcctc                  >  1:343387/1‑151 (MQ=255)
                                                                                                                            caATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTcc                 >  2:241943/1‑151 (MQ=255)
                                                                                                                                gCAACCAGACCAGCCCGGTTCGCTGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCCCCAGATCGGCTGGCAGCGTCCCCGACTCCAGCGTCACCGCCAGCGCCAGCTCCCCtcg              >  1:1026765/1‑148 (MQ=255)
                                                                                                                                  aaCCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAGGCGCCAGTGGGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCGCCGCCAGCGCCAGCTCTCCTGGGTg           <  1:907739/151‑1 (MQ=255)
                                                                                                                                     gagaccaaCCAGGTTCGATGCGGGGAGCGGGGATGTCCACGATGGGCCCCAGACGCGCCAAAGAAAAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGTCCGCGACACCGACTCCATCATCACCGCCACCGCCACCTCTCCTGGGt            >  1:370290/2‑147 (MQ=255)
                                                                                                                                     cagcccagccCGGTTCGATGCGGTGCGCGGTGCCGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAg        <  1:426818/151‑1 (MQ=255)
                                                                                                                                         ccagCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGAt      >  2:465112/1‑149 (MQ=255)
                                                                                                                                            gCCGGGTTCGAGGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCCCCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGaa  <  1:811305/150‑1 (MQ=255)
                                                                                                                                              ccGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATa     <  1:702476/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCATCTTCACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAA  >  NC_000913/2141686‑2141975

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