Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,728,796 C→T 19.0% D150D (GAC→GAT rnt → RNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,728,7960CT19.0% 46.9 / 5.1 21D150D (GAC→GATrntRNase T; exoribonuclease T; structured DNA 3' exonuclease; RNA processing; DNA repair
Reads supporting (aligned to +/- strand):  ref base C (8/9);  new base T (2/2);  total (10/11)
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

TACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGT  >  NC_000913/1728654‑1728939
                                                                                                                                              |                                                                                                                                               
tACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACActgctg                                                                                                                                          <  1:136167/150‑1 (MQ=255)
tACGTAAAGGTATTAAAGCGAGCGGCTGGAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACActgctg                                                                                                                                          <  1:78211/150‑1 (MQ=255)
            ttAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGATACTGCTGCACTGGCCGGGc                                                                                                                              >  2:9915/1‑150 (MQ=255)
            ttAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGATACTGCTGCACTGGCCGGGc                                                                                                                              >  2:9897/1‑150 (MQ=255)
                agcgagcgGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGc                                                                                                                          >  2:109364/1‑150 (MQ=255)
                agcgagcgGCTGTAACCGCGCCAGTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACGTTTGACACTGCTGCACTGGCCGGGCTGGc                                                                                                                          >  2:236226/1‑150 (MQ=255)
                                cgcgCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTa                                                                                                          >  1:78907/1‑150 (MQ=255)
                                     cATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTc                                                                                                     >  2:202411/1‑150 (MQ=255)
                                           ggtggCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGc                                                                                               <  2:179871/150‑1 (MQ=255)
                                                                tGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGATACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGa                                                                          <  1:9915/150‑1 (MQ=255)
                                                                tGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGATACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGa                                                                          <  1:9897/150‑1 (MQ=255)
                                                                tGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGa                                                                          <  2:100946/150‑1 (MQ=255)
                                                                       aGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGAc                                                                   <  1:109364/150‑1 (MQ=255)
                                                                       aGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGAc                                                                   <  1:236226/150‑1 (MQ=255)
                                                                       aGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGAc                                                                   <  1:139124/150‑1 (MQ=255)
                                                                                             aaCGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCgcgc                                             >  1:58130/1‑150 (MQ=255)
                                                                                                 cgcCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTa                                         <  2:123773/150‑1 (MQ=255)
                                                                                                 cgcCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTa                                         <  2:123777/150‑1 (MQ=255)
                                                                                                                   cccGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActgc                       >  1:113351/1‑150 (MQ=255)
                                                                                                                   cccGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCActgc                       >  1:113349/1‑150 (MQ=255)
                                                                                                                                        ttttGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCACTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGt  >  2:84384/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TACGTAAAGGTATTAAAGCGAGCGGCTGTAACCGCGCCATTATGGTGGCGCACAATGCCAATTTTGATCACAGCTTTATGATGGCCGCCGCAGAACGCGCCTCACTGAAACGTAACCCGTTCCACCCTTTCGCCACTTTTGACACTGCTGCACTGGCCGGGCTGGCACTCGGACAAACCGTATTGTCAAAGGCTTGCCAGACCGCTGGCATGGACTTCGACAGCACCCAGGCGCACTCCGCGCTGTACGACACCGAACGCACTGCTGTGCTGTTTTGTGAAATCGT  >  NC_000913/1728654‑1728939

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