Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,287,294 C→G noncoding (35/85 nt) tyrV ← tRNA‑Tyr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,2940CG86.7% 36.5 / ‑4.7 15noncoding (35/85 nt)tyrVtRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (5/8);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGA  >  NC_000913/1287050‑1287533
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGTGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg                                                                                                                                                                                                                                           >  1:32818/1‑251 (MQ=255)
           aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGc                                                                                                                                                                                                                                <  1:217996/251‑1 (MQ=255)
              gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGc                                                                                                                                                                                                                             <  1:149215/251‑1 (MQ=255)
                                               gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCCGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGGACTTCAAATTTTg                                                                                                                                                                                            >  2:232137/1‑251 (MQ=255)
                                                            cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTga                                                                                                                                                                               <  1:94257/251‑1 (MQ=255)
                                                                                                                gTCGAAGCTTCGCGTCCGTCCCCGCATGGGCAGAATAGTTGATTGCGGACTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGc                                                                                                                            <  2:149075/250‑1 (MQ=38)
                                                                                                                                                 aaTATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGAAGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGTGGTAATGCTTGAGATGGGGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAGCGTTCTCTCGCTTTCGCtc                                                                                          >  1:170200/1‑251 (MQ=255)
                                                                                                                                                         aTTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcgaa                                                                                  <  1:232137/251‑1 (MQ=255)
                                                                                                                                                           tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGATCCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAc                                                                                 <  2:157582/250‑1 (MQ=255)
                                                                                                                                                                aTTCGTTGGGAAGTTCAGGGACTTTTGAACGCGATGTTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTCGGCCGCTCGGGAACCCCACCAGGGGTAATTAAAATTTTGAGGTTAGGCTGGGGATGGTGGTGGGGGAAGGATTATTCGGCGTTTCGCGCCTCACCCTTCGGGCAGTTGCCTGTGGCAACGTTGTCTCGCTTTCGCtcgatgagaagcttagt                                                                           >  1:198631/1‑238 (MQ=25)
                                                                                                                                                                            gTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCa                                                               <  1:268777/251‑1 (MQ=255)
                                                                                                                                                                                          gAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATg                                                  >  1:111619/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gttgtcGGCAGATCTAGACTCTGCTCCCTTTGGCCGCTCGGGAACCTGCCGATGGGTAATTCAAATTTTGAGGTAATGATTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCtgt       <  2:261658/246‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGCGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGa  >  1:149957/1‑250 (MQ=255)
                                                                                                                                                                                                                                          ggCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGGTGCCTGTGGCAACGATCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACCTCGCAGTGGTGa  >  2:243116/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               
GGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGA  >  NC_000913/1287050‑1287533

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: