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

GGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGA  >  NC_000913/4606093‑4606355
                                                                                                                                                  |                                                                                                                    
ggATTTGAACCCCCACGTCCGTAAGGACATTAACTCCTGCAGCTAGCGCGTATACCAATTCCGCGACCTTATCACATTCATCTTACTTTTTTTGATCTCTCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                                                    <  4:187352/149‑1 (MQ=25)
                              tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                                                                     >  6:31097/1‑118 (MQ=255)
                              tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGaa                                                                                                                     <  5:31097/118‑1 (MQ=255)
                                acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                               >  3:21828/1‑122 (MQ=255)
                                acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                               <  4:21828/122‑1 (MQ=255)
                                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                                  >  3:342673/1‑99 (MQ=255)
                                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                                  <  4:342673/99‑1 (MQ=255)
                                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                                                                     <  2:161958/112‑1 (MQ=255)
                                                    tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                                                                     >  1:161958/1‑112 (MQ=255)
                                                          ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGACTCGTTTGGTGCGAGGGGGGGGACGTGCCCCCCCACCGCCGTAAGAACACTAACACATGAAGATAGCGAGTCTACCAATTGCGCCACCTTCGCCCAGGGCGAGCaa                                                          >  2:89609/1‑149 (MQ=25)
                                                              gCCACCTTCGCACAGTCATCTAAATTATTTAGAAATCGCGTCGTTTGGTGCTAGGGTGGTCACTTGCGCCCACAAGTCCGTAAGAACACTAACACATGAATCTAGAGCGTCTAGGAATTCCGCCCCCTTCGCGCGGTGCGAGCAatag                                                       >  2:95445/1‑147 (MQ=255)
                                                              gCCACCTTCGCACAGTCATCATACTTTATTTTATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGCACCCCCACGTCCGTAAGAACACAAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat                                                         <  1:95445/146‑1 (MQ=21)
                                                                 aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                                              <  7:134814/138‑1 (MQ=255)
                                                                 aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGACACCTTCGCCCAGTGCGa                                                              >  8:134814/1‑138 (MQ=255)
                                                                                               tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                >  5:109866/1‑106 (MQ=255)
                                                                                               tatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                <  6:109866/106‑1 (MQ=255)
                                                                                                atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                                  >  4:232450/1‑87 (MQ=255)
                                                                                                atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                                  <  3:232450/87‑1 (MQ=255)
                                                                                                atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                         <  7:397089/128‑1 (MQ=255)
                                                                                                atCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                         >  8:397089/1‑128 (MQ=255)
                                                                                                   gCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACAAACTCATAAAGCAAGCGAGTCAACCAATTACGACACCTTCGCACAGTGCGAGCAATATCCACGTGGTTTATGATGCGAGGGGGGGGCCTTGGAccc                 >  6:147597/1‑149 (MQ=18)
                                                                                                          ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                                   >  7:404547/1‑44 (MQ=38)
                                                                                                          ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                                   <  8:404547/44‑1 (MQ=38)
                                                                                                                  aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACAAACACCCGAAGCTAGCGCGTCGACCAATTCCGCCCCCTTCGCCCAGTGCGAGCAATATCAACGTGTGTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGa  <  5:147597/149‑1 (MQ=25)
                                                                                                                               gAACCCCCACGTCCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                <  4:361798/74‑1 (MQ=255)
                                                                                                                                   ccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                <  5:283770/70‑1 (MQ=18)
                                                                                                                                   ccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                                >  6:283770/1‑70 (MQ=18)
                                                                                                                                             gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAgg   <  4:235286/121‑1 (MQ=255)
                                                                                                                                             gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAgg   >  3:235286/1‑121 (MQ=255)
                                                                                                                                               aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                               >  7:300710/1‑59 (MQ=255)
                                                                                                                                               aaGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                               <  8:300710/59‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                    
GGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGA  >  NC_000913/4606093‑4606355

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: