Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 696,793 T→A 13.0% noncoding (38/75 nt) glnW ← tRNA‑Gln

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913696,7930TA13.0% 66.3 / 2.4 30noncoding (38/75 nt)glnWtRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base T (21/5);  new base A (2/2);  total (23/7)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCC  >  NC_000913/696644‑696937
                                                                                                                                                     |                                                                                                                                                
tatTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                                                                                                  >  2:1895271/1‑150 (MQ=255)
    cGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCa                                                                                                                                              >  2:726943/1‑150 (MQ=255)
     gATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaa                                                                                                                                             >  1:1200553/1‑150 (MQ=255)
       tACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaa                                                                                                                                           >  1:1679744/1‑150 (MQ=255)
          aTTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                                                                                                        >  2:226433/1‑150 (MQ=255)
             cAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                                                                                     >  1:1823370/1‑150 (MQ=255)
                       tggCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc                                                                                                                           >  2:824595/1‑150 (MQ=255)
                           tACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTgg                                                                                                                       >  1:1691382/1‑150 (MQ=255)
                             cgacgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCg                                                                                                                     >  1:2009713/1‑150 (MQ=255)
                             cgacgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCg                                                                                                                     >  2:1006613/1‑150 (MQ=255)
                               acgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGAt                                                                                                                   >  1:1221366/1‑150 (MQ=255)
                               acgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGAt                                                                                                                   <  2:1988861/150‑1 (MQ=255)
                                    aTTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATAcccc                                                                                                              >  2:1157452/1‑150 (MQ=255)
                                        gAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACAAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGAGATACCCCaaaa                                                                                                          >  1:123891/1‑148 (MQ=255)
                                           ccTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGTATTCGAACCAGGGAATGCCGGTATCAAAACCCGTGCCATAACCCCTGGCCAAAACCCAAtcagcc                                                                                                      >  2:778162/1‑145 (MQ=25)
                                                 aCCCCTTCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCgg                                                                                                 <  1:1702654/150‑1 (MQ=255)
                                                                gTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCg                                                                                  >  1:1930084/1‑150 (MQ=255)
                                                                       gCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAACCGGGCGGTGACCCGCTTACTCGAagaaga                                                                           >  1:11913/1‑150 (MQ=255)
                                                                             aCCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCATAAGACGGGCGGTGAAACGTTTAGTCGACGATGCTCGCTg                                                                     >  1:679874/1‑150 (MQ=18)
                                                                              ccAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAAACGAGCGGTGAAACGCTTAATCGAAGAAGATCGATgg                                                                    >  1:612558/1‑150 (MQ=21)
                                                                                    gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATACACGGCCGGTTAACCGATAACTCGTAGAAGATGGATGGGGTAcc                                                              >  1:1646118/1‑150 (MQ=21)
                                                                                     aGCTACGTAGCCGGATTGTTTCTTCGATGGCTGGGTTACCTTGATGCGATCCAGGGATTGCCGGTATCAAAACCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGTGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCt                                                             <  2:1288963/150‑1 (MQ=25)
                                                                                      gCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTg                                                            <  2:1930084/150‑1 (MQ=255)
                                                                                          cGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAGAGACGGGCGGTGACCCGAGTAGTCGAAGATGATGGCTGGGGTACCtcgctg                                                        >  1:955801/1‑145 (MQ=18)
                                                                                                                    tGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCaaaa                              <  2:1691382/150‑1 (MQ=255)
                                                                                                                             aggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaag                                                                               >  2:281735/2‑92 (MQ=14)
                                                                                                                             aggATTCGAACCTCGGAATGCCGGAATCAGAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaag                                                                               <  1:281735/91‑1 (MQ=17)
                                                                                                                              ggATTCGAACCTCGGAATGCCGGAATCATAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaagaa                                                                             <  1:774539/93‑1 (MQ=17)
                                                                                                                              ggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCgaagaa                                                                             >  2:774539/1‑93 (MQ=17)
                                                                                                                                    gAACCAGGGAATGCCGGTATCAAAAACCGGTACCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGc              >  1:1693612/1‑150 (MQ=255)
                                                                                                                                                gCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAAAACCGGGCGGTGAAACGATTACTCGAAGAAGATGGCTGGGGTACCTGGCTACCAACCAGGCAATGCCGGTATCAACAACAGGTGCCTAACAGATTGGCGATAccc  >  2:241139/1‑150 (MQ=37)
                                                                                                                                                     |                                                                                                                                                
TATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCC  >  NC_000913/696644‑696937

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: