Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,239 A→G 17.0% 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,2390AG17.0% 55.1 / 5.2 31noncoding (48/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (2/3);  total (14/17)
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

AACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  NC_000913/4606100‑4606368
                                                                                                                                           |                                                                                                                                 
aaCCCCCACGGCAGTAATGACACTAACACCTGATGCTATCCAGTATACCAATTCCGCCACCAAACCACAGGAATCAAACTTTATTTCATATCCCCTCTTTGGGGGCGCGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAacac                                                                                                                         <  8:226159/149‑1 (MQ=16)
     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                                              <  8:462477/140‑1 (MQ=255)
     ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                                              >  7:462477/1‑140 (MQ=255)
        cGTCCTGAATGACACTAACTCCTGAAGCTTGCTCGTATACCTATTCCGGCACCTTCCCACAGGCATCTTACTATTTTTGATATCGCCTCTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                                                                                  <  3:454021/149‑1 (MQ=21)
                             cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                                               <  1:202603/115‑1 (MQ=255)
                             cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATAGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                                               >  2:202603/1‑115 (MQ=255)
                                               ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                                      >  1:153167/1‑122 (MQ=255)
                                               ccAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                                      <  2:153167/122‑1 (MQ=255)
                                                    tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                                                     >  1:101652/1‑134 (MQ=255)
                                                    tCCGCCACCTTCGCACAGGCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                                                     <  2:101652/134‑1 (MQ=255)
                                                                           gtACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                                                    <  2:152925/111‑1 (MQ=255)
                                                                               ttttttttATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                                                            <  8:462915/84‑1 (MQ=255)
                                                                               tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                                      <  3:350310/90‑1 (MQ=255)
                                                                               tttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                                                                      >  4:350310/1‑90 (MQ=255)
                                                                                 tttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                                                      >  6:401322/1‑104 (MQ=255)
                                                                                      gATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAg                                               >  5:195214/1‑138 (MQ=255)
                                                                                      gATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAg                                               <  6:195214/138‑1 (MQ=255)
                                                                                              gattattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                             <  5:377867/95‑1 (MQ=255)
                                                                                              gataattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                             >  6:377867/6‑100 (MQ=0)
                                                                                                 aattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                                            >  6:87016/3‑98 (MQ=9)
                                                                                                 aattGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                                            <  5:87016/96‑1 (MQ=38)
                                                                                                                ggggACTTGAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                                               <  7:180447/80‑1 (MQ=21)
                                                                                                                        gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                                     >  2:88833/1‑98 (MQ=255)
                                                                                                                        gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                                     <  1:88833/98‑1 (MQ=255)
                                                                                                                               cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                                        <  2:176681/88‑1 (MQ=255)
                                                                                                                               cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                                        >  1:176681/1‑88 (MQ=255)
                                                                                                                                        aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa  >  7:363830/1‑133 (MQ=255)
                                                                                                                                        aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa  <  8:363830/133‑1 (MQ=255)
                                                                                                                                          gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                  >  3:444153/1‑115 (MQ=255)
                                                                                                                                          gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTaa                  <  4:444153/115‑1 (MQ=255)
                                                                                                                                          gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAacac      >  6:354244/1‑127 (MQ=255)
                                                                                                                                           |                                                                                                                                 
AACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  NC_000913/4606100‑4606368

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: