Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,392,506 C→G 10.6% noncoding (35/76 nt) glyX → tRNA‑Gly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,392,5060CG10.6% 92.7 / 6.7 38noncoding (35/76 nt)glyXtRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base C (16/18);  new base G (2/2);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

AGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTT  >  NC_000913/4392373‑4392649
                                                                                                                                     |                                                                                                                                               
acttGGTAGAGCACGACCTTGCCAAGCGCGGGGTCGCGAGATCGAGTCGCGTTTCCCGCTCCAAAATTAGAAAAGTGCGGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                  <  1:102589/147‑1 (MQ=14)
aGTTTGTAGCGCACGACCTTGCCAACGGCGGGGTCGCGAGTTCGAGGCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCGGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTcgc                                                                                                                                  <  5:91647/149‑1 (MQ=14)
                ccTTGGCAAGGTCGTGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGtt                                                                                                                  <  7:347172/149‑1 (MQ=255)
                      cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                                                              >  8:48425/1‑147 (MQ=39)
                      cAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTccc                                                                                                              <  7:48425/147‑1 (MQ=39)
                                                           tCCAAAATTTGAAAATTGATGCAAAGCACAGCCCACCCAAGCCGGAATAGGTAAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGGGAGTTCGAGTCTAGTTTACAGCTGCAAAATTTGAAAGTGCTGTAGGGCAgagagagagc                                                                       >  2:403742/1‑143 (MQ=17)
                                                           tCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGa                                                                             <  1:403742/143‑1 (MQ=255)
                                                             ctaaaTTTGAAAATTACTGCAACTCAAAGACGAAGCAAGCGGGAATAGCTCAGTTGGTAGAGCACGTCCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCaa                                                                     <  6:52984/147‑1 (MQ=17)
                                                             cAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCGCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCGccc                                                                       >  7:71535/1‑147 (MQ=255)
                                                              aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                                                                                 >  7:163095/1‑120 (MQ=255)
                                                              aaaaTTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTg                                                                                                 <  8:163095/120‑1 (MQ=255)
                                                                       aaaGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                                                          >  8:128192/1‑86 (MQ=38)
                                                                       aaaGTGCTGAAAAGCACAGACCACCCATGCGTGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAg                                                                              <  7:241991/130‑1 (MQ=255)
                                                                           tgctgcAAAGCACAGACCACCCATGTTTGAATAGCTCAGTTGCTTGAGCACTACCTTGCCAATGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCAGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                       <  8:71535/149‑1 (MQ=17)
                                                                               gcAATGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                       <  2:16476/65‑1 (MQ=18)
                                                                               gcAACGCACAGACCACCCACGCGGCAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTg                                                   <  5:333623/149‑1 (MQ=255)
                                                                               gcAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                                                       >  1:16476/1‑65 (MQ=39)
                                                                                  aaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                         <  2:337719/60‑1 (MQ=33)
                                                                                  aaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                         >  5:145926/1‑60 (MQ=33)
                                                                                  aaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                         >  1:337719/1‑60 (MQ=33)
                                                                                  aaGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCg                                                                                                                                         <  6:145926/60‑1 (MQ=33)
                                                                                              cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGCTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                       <  7:413781/130‑1 (MQ=255)
                                                                                              cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCAAAGGTCGGGGTCGCGAGTTCGAGTCTCGCTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCa                                                       >  8:413781/1‑130 (MQ=255)
                                                                                                 aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGc                                                          <  1:103504/124‑1 (MQ=255)
                                                                                                 aaGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGc                                                          >  2:103504/1‑124 (MQ=255)
                                                                                                    cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGa                                                                             >  2:55779/1‑102 (MQ=255)
                                                                                                    cGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGa                                                                             <  1:55779/102‑1 (MQ=255)
                                                                                                          tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCt                                                                                         <  6:320298/84‑1 (MQ=37)
                                                                                                          tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCt                                                                                         >  5:320298/1‑84 (MQ=37)
                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                                                   <  5:108265/99‑1 (MQ=255)
                                                                                                                 gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGAACACCCAAGc                                                                   >  6:108265/1‑99 (MQ=255)
                                                                                                                     gTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAGTACTGAAAGGAACAGACCACCCAAGCGGGAATAGCTGAGGTGGGAGAGCACGAACTGGCCAAGGTCGGGGTAGAGAGt               >  8:1543/1‑147 (MQ=11)
                                                                                                                      tAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAgg                                                                                   >  4:132376/1‑78 (MQ=255)
                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAcc                                      >  7:22892/3‑118 (MQ=255)
                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAcc                                      <  8:22892/116‑1 (MQ=255)
                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTc                           <  2:131292/127‑1 (MQ=255)
                                                                                                                           cccGACCTTGGCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTc                           >  1:131292/3‑129 (MQ=255)
                                                                                                                                ccTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGTCACAGACCACCCAAGCGGGAATAGCGCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGgtt  >  8:68896/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                               
AGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTT  >  NC_000913/4392373‑4392649

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