Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,841,780 A→G E97E (GAA→GAG nudG → 5‑hydroxy‑CTP diphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,841,7800AG97.3% 122.6 / ‑4.1 37E97E (GAA→GAGnudG5‑hydroxy‑CTP diphosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/16);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCAT  >  NC_000913/1841647‑1841915
                                                                                                                                     |                                                                                                                                       
tACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAg                                                                                                                                   >  1:96269/1‑140 (MQ=255)
  cgcgAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGgc                                                                                                                                 >  1:565918/1‑140 (MQ=255)
    cgAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGgcgc                                                                                                                               <  2:96269/140‑1 (MQ=255)
         aaCTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGctggtc                                                                                                                          <  1:265760/140‑1 (MQ=255)
              ggCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATc                                                                                                                                     <  2:232161/124‑1 (MQ=255)
              ggCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATc                                                                                                                                     >  1:232161/1‑124 (MQ=255)
                   cGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCAc                                                                                                                >  2:149127/1‑140 (MQ=255)
                             gtgGGTGAATATGTTGCCAGCCATCAGGGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGgc                                                                                                      >  1:113750/1‑140 (MQ=255)
                             gtgGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGgc                                                                                                      >  2:278869/1‑140 (MQ=255)
                             gtgGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGgc                                                                                                      >  2:746764/1‑140 (MQ=255)
                                          ttGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCg                                                                                         <  2:227610/140‑1 (MQ=255)
                                          ttGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCg                                                                                         >  1:227610/1‑140 (MQ=255)
                                           tGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGc                                                                                        >  2:534221/1‑140 (MQ=255)
                                               agccaTCAGGGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGc                                                                                    <  2:113750/140‑1 (MQ=255)
                                                    tCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCg                                                                                         >  1:152975/1‑130 (MQ=255)
                                                    tCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCg                                                                                         <  2:152975/130‑1 (MQ=255)
                                                           gaAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACAt                                                                        <  1:3584/140‑1 (MQ=255)
                                                           gaAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACAt                                                                        <  2:822244/140‑1 (MQ=255)
                                                            aaGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACAtt                                                                       <  2:174745/140‑1 (MQ=255)
                                                              gTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTcc                                                                     <  1:746764/140‑1 (MQ=255)
                                                               tttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACAt                                                                        <  1:500994/136‑1 (MQ=255)
                                                               tttCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACAt                                                                        >  2:500994/1‑136 (MQ=255)
                                                                   ggggCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCattat                                                                <  1:149127/140‑1 (MQ=255)
                                                                                          ccTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTAcgc                                         <  2:395574/140‑1 (MQ=255)
                                                                                            tGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACgcgc                                       >  2:140233/1‑140 (MQ=255)
                                                                                              gCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGcgccg                                     >  1:505549/1‑140 (MQ=255)
                                                                                                   tACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGa                                >  1:543901/1‑140 (MQ=255)
                                                                                                   tACCCGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGa                                >  1:325499/1‑140 (MQ=255)
                                                                                                       cGACTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAg                            <  2:325499/140‑1 (MQ=255)
                                                                                                         aCTTCCACGGGACGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCg                          <  2:543901/140‑1 (MQ=255)
                                                                                                                    aCGTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTa               <  2:235869/140‑1 (MQ=255)
                                                                                                                      gTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAg             >  2:124965/1‑140 (MQ=255)
                                                                                                                      gTTACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAg             >  2:279790/1‑140 (MQ=255)
                                                                                                                        tACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAggg           >  2:243261/1‑140 (MQ=255)
                                                                                                                        tACAGGCACATGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAggg           >  1:219153/1‑140 (MQ=255)
                                                                                                                             gCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTg      >  2:550462/1‑140 (MQ=255)
                                                                                                                                 aTGAGCATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGtcat  >  1:755096/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TACGCGAAGAACTGGGCATCGAAGCAACTGTGGGTGAATATGTTGCCAGCCATCAGCGAGAAGTTTCGGGGCGGATTATCCATCTTCATGCCTGGCACGTACCCGACTTCCACGGGACGTTACAGGCACATGAACATCAGGCGCTGGTCTGGTGCTCACCTGAAGAGGCGCTGCAATATCCGCTGGCCCCTGCTGACATTCCATTATTAGAGGCGTTTATGGCTTTACGCGCCGCCAGACCAGCGGATTAGTGCTAAGGGTTTTGTCAT  >  NC_000913/1841647‑1841915

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