Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,193 T→C 9.1% F19F (TTT→TTC)  tgt → tRNA‑guanine transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913426,1930TC 9.1% 111.0 / 2.2 44F19F (TTT→TTC) tgttRNA‑guanine transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (13/27);  new base C (2/2);  total (15/29)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

ATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGA  >  NC_000913/426047‑426339
                                                                                                                                                  |                                                                                                                                                  
aTCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGa                                                                                                                                                  <  2:652105/149‑1 (MQ=255)
     cAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTg                                                                                                                                             <  2:425651/149‑1 (MQ=255)
     cAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAAATTGATCTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTg                                                                                                                                             <  1:354928/149‑1 (MQ=255)
          aaTTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTa                                                                                                                                        <  2:312753/149‑1 (MQ=255)
            ttAATGAGCGCGGCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGt                                                                                                                                      <  1:77947/149‑1 (MQ=255)
                  aGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCCGGTCTTTGATCGTGGCGTAGTGGAAAc                                                                                                                                >  2:278184/1‑149 (MQ=255)
                      cgTCGGGGAGTAAATCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCt                                                                                                                            >  1:547725/1‑149 (MQ=255)
                              aGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTAt                                                                                                                    <  1:96888/149‑1 (MQ=255)
                              aGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTAt                                                                                                                    <  2:464301/149‑1 (MQ=255)
                              aGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTAt                                                                                                                    <  1:53347/149‑1 (MQ=255)
                                  aTTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTCGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCt                                                                                                                >  2:15483/1‑149 (MQ=255)
                                  aTTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTCGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCt                                                                                                                >  2:770027/1‑149 (MQ=255)
                                    tCCGCGGCGCTGGTTCAAAACGTTGGACTGTTTCTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGt                                                                                                              <  1:269026/149‑1 (MQ=255)
                                      cgcgGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTg                                                                                                            <  2:496109/149‑1 (MQ=255)
                                      cgcgGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTg                                                                                                            <  2:107015/149‑1 (MQ=255)
                                      cccggcgCTGGTTTAAAACGTTGGACTGTGTTTCTTACGTAGTGGAGAAAAAATGAAATTTGAACGGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTg                                                                                                            <  2:256935/147‑1 (MQ=255)
                                       gcggcgCTGGTGTAAAGCGGTAGACTGTTTTGCTGACGGAGGGGTGAAAAAACGAAAGTTGAACTGGACACCACCGACGGTCGCGGACGCCGTCGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGCTgg                                                                                                           <  2:307253/149‑1 (MQ=255)
                                         ggcgCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAACAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCCGGTCTGTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCa                                                                                                         >  1:262517/1‑149 (MQ=255)
                                            gcTGGTTTAAAACGTTGGACTGTTTTTCTGACGGAGTTGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCt                                                                                                      <  1:255472/149‑1 (MQ=255)
                                                    aaaaCGTTGGACTGTCTTTCTGACGTAGTGGAGAAAAAAAGAAATTTGAACTGGACACCACAGACGGACGCGCAAGCCGTGGCCGCCTGGCCATTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGGCTGTTGACACCTACGGCAcc                                                                                              >  2:325675/1‑149 (MQ=255)
                                                              aCTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGa                                                                                    <  1:277900/149‑1 (MQ=255)
                                                              aCTGGTGTTCTGCCGTAGGGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGa                                                                                    <  2:6143/149‑1 (MQ=255)
                                                               cTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCATAAAAGGGAt                                                                                   <  1:300101/149‑1 (MQ=255)
                                                                     ttCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGcc                                                                             >  2:699271/1‑149 (MQ=255)
                                                                     ttCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGcc                                                                             >  2:748999/1‑149 (MQ=255)
                                                                        tGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGga                                                                          <  1:134960/149‑1 (MQ=255)
                                                                                 ggAGAAAAAATGAAATTTCAACTGGACACCACCGACGGTCGCGCACGCCGCTGCCGGCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGa                                                                 <  1:783236/149‑1 (MQ=255)
                                                                                  gagaAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTCGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGaa                                                                <  1:15483/149‑1 (MQ=255)
                                                                                  gagaAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTCGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGaa                                                                <  1:770027/149‑1 (MQ=255)
                                                                                   agaAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCACACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGACCCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAg                                                               >  2:365458/1‑149 (MQ=255)
                                                                                     aaaaaaTGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGcc                                                             >  1:298132/1‑149 (MQ=255)
                                                                                             aaaTTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGgcgc                                                     <  1:752466/149‑1 (MQ=255)
                                                                                              aagtttaaCCGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGgcgcg                                                    <  2:342232/143‑1 (MQ=255)
                                                                                                tttGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGAGGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCa                                                  >  1:786821/1‑149 (MQ=255)
                                                                                                       tGGACACCACGGACGGTCGCGCGCGCCGTGGCGGCCTGGTCTTTGACCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTACAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATc                                           <  2:402714/149‑1 (MQ=255)
                                                                                                       tGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGCTGACGCCGGAAGAAGTTGAAGCCACTCGCGAGCAAATTATc                                           >  2:500714/1‑149 (MQ=255)
                                                                                                            accaccGACGGTCGCGCACGCCGTGGCCGCCTGGGCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGCAAAATGGATGACGTCGGAGGCAGGTGAAGCCACTGGCGCGCGAATTATCCCCgg                                      >  1:457949/1‑149 (MQ=255)
                                                                                                               accGACGGTCGCGCACGCCGTGGCCGACTGGAGCATTTTCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCaa                                   >  2:772441/1‑149 (MQ=255)
                                                                                                                      gTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTc                            <  2:252291/149‑1 (MQ=255)
                                                                                                                               gCCGTGGCCGCCTGGTCTTTGCTCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTgg                   <  1:259872/149‑1 (MQ=255)
                                                                                                                               gCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTgg                   <  1:396949/149‑1 (MQ=255)
                                                                                                                                           tGGTCCTTGCACGTCGCCTGGGGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAATGTGGATGACGCCTGAACAAGTTGAAGCCACTTGCGCGCAAAGTATCCTCGGCAACACCTTCCACCTGTCGCTGCGCCCGGGc       <  2:580828/149‑1 (MQ=255)
                                                                                                                                             gTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCa     <  1:563353/149‑1 (MQ=255)
                                                                                                                                             gTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCa     >  2:458982/1‑149 (MQ=255)
                                                                                                                                                tttGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGa  >  1:580687/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAATGAAATTTGAACTGGACACCACCGACGGTCGCGCACGCCGTGGCCGCCTGGTCTTTGATCGTGGCGTAGTGGAAACGCCTTGTTTTATGCCTGTTGGCACCTACGGCACCGTAAAAGGGATGACGCCGGAAGAAGTTGAAGCCACTGGCGCGCAAATTATCCTCGGCAACACCTTCCACCTGTGGCTGCGCCCGGGCCAGGA  >  NC_000913/426047‑426339

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑