Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,239 A→G 15.7% noncoding (48/87 nt) leuP ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,2390AG15.7% 51.3 / 3.6 26noncoding (48/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base G (2/2);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAA  >  NC_000913/4606100‑4606352
                                                                                                                                           |                                                                                                                 
aaCCCCCACGTCCTTAATTACACTAAAACCTGAAGATTCATCGTATACCAATTCCGCCACCATCCAACAGTCATGATACTTTTTTTGATATACCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                            <  8:10058/147‑1 (MQ=21)
                                            cTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                               >  7:55698/1‑148 (MQ=255)
                                               ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                                 <  4:181840/95‑1 (MQ=255)
                                               ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                                 >  3:181840/1‑95 (MQ=255)
                                                    tCCGCCACCTTCGCACAGTCATCTGACTATAATTGATATATAATCGGTTGGTGCGAGGGGGGGGACGTGCACCCCCACGTCCGTAAGAACAATAACACCTGACGCTAGCGCGACTACCAATTACTCAACCTTCGCCCAGTGCGGGCaa                                                       >  6:329618/1‑148 (MQ=17)
                                                            cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                           <  6:44456/120‑1 (MQ=255)
                                                            cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                           >  5:44456/1‑120 (MQ=255)
                                                                       cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTTACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                    >  1:87663/1‑148 (MQ=255)
                                                                       cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                                     <  3:336658/131‑1 (MQ=255)
                                                                       cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                                     >  4:336658/1‑131 (MQ=255)
                                                                             aCTTTTTTGGATCTCGCCTCGTTTTGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                 <  5:329618/145‑1 (MQ=255)
                                                                              cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTTACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggg                             <  2:87663/148‑1 (MQ=255)
                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                                          <  7:356903/69‑1 (MQ=255)
                                                                                ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                                          >  8:356903/1‑69 (MQ=255)
                                                                                         atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                  >  8:330202/1‑116 (MQ=255)
                                                                                         atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                  <  7:330202/116‑1 (MQ=255)
                                                                                                 gTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                               <  1:7370/63‑1 (MQ=37)
                                                                                                 gTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                               >  2:7370/1‑63 (MQ=37)
                                                                                                        gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            >  1:426161/1‑91 (MQ=18)
                                                                                                        gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGATAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            <  2:426161/91‑1 (MQ=16)
                                                                                                         cccgggtggcgACTTTAACCCCCACTTCCGTAAGAACACTAACACCCGAAGCTTCCCCGTCTACCAATTCCGCCCCCTTCGCCCAGTCCTTGAAATATCAACGTGTTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa  <  8:55698/145‑1 (MQ=21)
                                                                                                         cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            <  5:407083/90‑1 (MQ=255)
                                                                                                         cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                            >  6:407083/1‑90 (MQ=255)
                                                                                                                             cccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCACCGTGGTGTTTGGTGCGAGGGGGGGGCCTTGAACCCCCACGTCCGTaa  >  7:48046/1‑128 (MQ=255)
                                                                                                                               cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                                            <  6:304371/52‑1 (MQ=35)
                                                                                                                               cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                                            >  5:304371/1‑52 (MQ=35)
                                                                                                                                           |                                                                                                                 
AACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAA  >  NC_000913/4606100‑4606352

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: