Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,725,176 T→G E48A (GAG→GCG)  glyQ ← glycine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,725,1760TG95.8% 65.2 / NA 24E48A (GAG→GCG) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (13/10);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGC  >  NC_000913/3725020‑3725412
                                                                                                                                                            |                                                                                                                                                                                                                                            
ntgnggtntnncnncttnttatTGTCCGGCGATGGCTTANNGACCTCCTGNAACGTATAGGAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGGNGGGCGGTGANNAGNCGGAANANNNGNNNCCNCCATAGGTTCTGGCCCCAGNGCNNGCAGACNGGTCATTGGGTGAGAGGTTCCCgcgc                                                                                                                                                                                                          <  2:115289/174‑1 (MQ=255)
          gCTCCTGAATATTGTCCGGCGATGGCTTAATGGCCACCTGGAACTGATAGTCGTGCGATAAACGCTTGGGGTTTTCGGCGTTGCGACGATCGCTCGGGCGAAGAGAAGGCTGTACATAAGCTGCAGCGATCGGTTCTGGCGCCAGCGAGCGCAGACAGGTCATTGGGAGCGAGGTTCCCGCGCCGAc                                                                                                                                                                                                      >  1:288962/1‑187 (MQ=255)
              cTGAATGTTGTCCGGCAATGGCTGAATGCCCACCTGAAACTGATATTAGTGATGTAAACGGTAGGGGTTTTCACCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCAATCGGTTCTGGCCCCAGAGCGCGCAGACAGGTCATTGGGTGAGAAGATCCCGCGCCGACTTACATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCACGTAATCCTGTAAGGTCAGGATc                                                                                                                                   >  1:429412/1‑250 (MQ=255)
                                             accTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTg                                                                                                   <  1:177057/251‑1 (MQ=255)
                                                 ggAACTGATAGTAGTGATGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTTGAAGTTCCTGGTATCAAACTTTTGCa                                                                                                 <  1:137847/249‑1 (MQ=255)
                                                                 cTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCACCCGCCATCGGTTCTGGCCCCCGCGCGCGCAGACAGGCCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGAt                                                                               >  1:290457/1‑251 (MQ=255)
                                                                        cGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCTAGCCCAGTAATCCTGGAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGGATCAAACTTGTGCATATTATTTCGTGCTGGATACGtgtg                                                                        <  1:343123/251‑1 (MQ=255)
                                                                             ggggTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATAGGTTCTGGCCCAAGAGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATAAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAttt                                                                   >  1:367102/1‑251 (MQ=255)
                                                                              gggTTTTCGCCGTAGCGACCATCGGACGGGCGCCGAGAAGGCTGAACATAAGCAGCCGCCCTCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTg                                                                                                                           >  1:157461/1‑194 (MQ=255)
                                                                                         gTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATATGCAGCCGCCATCGGTTCTGCCCCCAGCGCACGCAGACAGGTCATTGGGTGAGAGGTTCCAGCGATGACTTCTATGTCCACAGGTTGAACAATGGTGCAGCCATGGCGAGCCCAGTAATCCTGTAAGGTCCGGATCAAGCCCTGGACGGTCCTGGTCTCAAACTTTTGCATATTATTTCGt                                                                                      >  1:341059/1‑220 (MQ=255)
                                                                                            gCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAggg                                                         <  1:365860/246‑1 (MQ=255)
                                                                                                     ggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCCGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGc                                           >  1:158524/1‑251 (MQ=255)
                                                                                                     ggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCaa                                                                                                                                 >  1:169172/1‑165 (MQ=255)
                                                                                                            ngacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggatgg                                    >  1:317079/2‑251 (MQ=255)
                                                                                                             gacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGATAGGTTCCCGCGCCGACGTCCATg                                                                                                                                                                                               <  1:420869/95‑1 (MQ=255)
                                                                                                                    aGGCTGAACATAAGCAGCCGCCATCGGTTNTGGCCNCNGNNCGNNCAGACAgg                                                                                                                                                                                                                                  <  2:322298/53‑1 (MQ=255)
                                                                                                                    aGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAgg                                                                                                                                                                                                                                  >  1:322298/1‑53 (MQ=255)
                                                                                                                      ncNGANNATAANCAGCCGCCATCGGTTCTGGCCCCANCGCGCGCAGACAGGTCATTGGGTGAGAGGTTNCNGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGANCCCAGTAATCNTGTNNGGTCNGNANCAAGCCCTGGAAGGTCCTGGTATCAAACTNTTGCATATTATTTCGTGCTGGATNNGTNNNGANNTAaag                                                               <  2:367102/213‑1 (MQ=255)
                                                                                                                       cTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCAtn                                                                                                                                                                                               <  1:300250/85‑2 (MQ=255)
                                                                                                                       cTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATg                                                                                                                                                                                               >  2:300250/1‑85 (MQ=255)
                                                                                                                          aaCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAAGTGCCCGCGCCGACTTACATGTCCAATggg                                                                                                                                                                                      <  1:144617/91‑2 (MQ=255)
                                                                                                                                     ccgccATCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCgcgc                                                                                                                                                                                                          >  1:272171/1‑60 (MQ=255)
                                                                                                                                           tCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATg                                                                                                                                                                             <  1:29590/83‑1 (MQ=255)
                                                                                                                                           tCGGTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATg                                                                                                                                                                             >  2:29590/1‑83 (MQ=255)
                                                                                                                                              gTTCTGGCCCCAGCGCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGc  <  1:71325/251‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                                                                                                            
CCGAGGTACAGCTCCTGAATATTGTCCGGCGATGGCTTAATGACCACCTGGAACTGATAGTAGTGCTGTAAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGAAGATATACAGTACGAAACGGGAAAAAGCAGGGC  >  NC_000913/3725020‑3725412

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