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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,3540GA11.5% 138.9 / 6.8 52noncoding (48/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base A (4/2);  total (30/22)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTG  >  NC_000913/4606263‑4606501
                                                                                           |                                                                                                                                                   
gTCTACCAATTCCGCCACCTTCGCCAACAGAGAGAAATATAAAAGTGTTTTTTCGTGCGAGGCGGGGTCGTTGGACCCACACGTCCGTAAGGACACTAACACCTGAAGCTAGAGCGTCTACCAGTTCCGACAGCTTCGCATACCATCaa                                                                                            >  8:44175/1‑149 (MQ=18)
        aTTCCGCCACCTTCGCCCAGTGCGAGGTATATCAACGTGGTTGTTGGTGCGAGGGGGGGGACTTGCAGCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaa                                                                                     >  8:354186/1‑148 (MQ=21)
                  cTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc                                                                           >  2:170276/1‑148 (MQ=255)
                   ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGAACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCt                                                                                                                          <  5:295618/100‑1 (MQ=255)
                   ttCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGAACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCt                                                                                                                          >  6:295618/1‑100 (MQ=255)
                           gTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg                                                                 >  8:69621/1‑149 (MQ=255)
                             ccGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCTCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgc                                                               <  7:44175/148‑1 (MQ=255)
                               gAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCAGCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                                                              >  2:152718/1‑148 (MQ=255)
                                   aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc                                                         <  1:170276/149‑1 (MQ=255)
                                    atatCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCATGGAGGCGCATTc                                                         >  2:94617/1‑148 (MQ=255)
                                       tCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                                                                     <  6:400247/101‑1 (MQ=255)
                                       tCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                                                                     >  5:400247/1‑101 (MQ=255)
                                       tCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                                                            <  7:331442/126‑1 (MQ=255)
                                           cggcgggtttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtt                                                 <  2:158052/142‑1 (MQ=255)
                                               gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtttt                                               >  2:89002/1‑147 (MQ=255)
                                               gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGtttt                                               >  3:272237/1‑147 (MQ=255)
                                                 ttttGGTGCGAGGGGGGGGACTTGAACCCCCACGTACGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGAATACCATCAATTCTTAAAAAGAATTGCAACCACGGAGGCGCATTCTAGTGTttttaatc                                           >  6:371397/1‑145 (MQ=21)
                                                                 gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                                                                                  <  7:230217/94‑1 (MQ=255)
                                                                 gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                                                            >  7:72771/1‑100 (MQ=255)
                                                                 gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                                                                            <  8:72771/100‑1 (MQ=255)
                                                                 gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGtt                            >  4:295628/1‑148 (MQ=255)
                                                                 gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACAATCAAATATTaaaaa                                                                                  >  8:230217/1‑94 (MQ=14)
                                                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCTTCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTaa                          >  6:295324/1‑149 (MQ=255)
                                                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAATAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTaa                          <  1:238809/149‑1 (MQ=255)
                                                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTc                                      >  1:256797/1‑137 (MQ=255)
                                                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTc                                      <  2:256797/137‑1 (MQ=255)
                                                                  ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTaa                          <  1:89002/149‑1 (MQ=255)
                                                                    aCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAAt                         >  6:172245/1‑148 (MQ=255)
                                                                      ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCTCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTAt                      >  1:289915/1‑149 (MQ=255)
                                                                         aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                   <  5:172245/149‑1 (MQ=255)
                                                                           cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg                                                      <  3:179211/112‑1 (MQ=255)
                                                                           cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg                                                      >  4:179211/1‑112 (MQ=255)
                                                                             cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaa                                                                                   >  4:240887/1‑81 (MQ=255)
                                                                             cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaa                                                                                   <  3:240887/81‑1 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGaa                                                                               <  3:171136/84‑1 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc                                                                           <  8:315407/88‑1 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc                                                                           >  7:315407/1‑88 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                               <  4:238239/132‑1 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATa                               >  3:238239/1‑132 (MQ=255)
                                                                              ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAAGTATCGACagag              >  4:241117/1‑149 (MQ=255)
                                                                              ccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCATAACTTTAAAAGTATTGATACCACGGAGGCGCATTCATTTCGCTGACAGCTGTTACTCAAAAGTTAATGATCGagctag              >  6:242634/1‑144 (MQ=12)
                                                                              ccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGaa                                                                               >  4:171136/1‑84 (MQ=255)
                                                                               cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                                                                                  <  6:13520/80‑1 (MQ=255)
                                                                               cACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                                                                                  >  5:13520/1‑80 (MQ=255)
                                                                                aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGgt            <  7:69621/149‑1 (MQ=255)
                                                                                 cGTCCGTAAGGTAAGTAACTCCTGAAGATAGAGCGTTTACCAATTGCGCTTACTTCTCATACCTTAAATTATTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTAGTCAATTGTTAATTa                       <  8:165627/137‑1 (MQ=255)
                                                                                 cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTa                       >  7:165627/1‑137 (MQ=255)
                                                                                 cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGgtg           <  6:170900/149‑1 (MQ=255)
                                                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                   >  7:404408/1‑137 (MQ=255)
                                                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGa                   <  8:404408/137‑1 (MQ=255)
                                                                                     cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAAt       >  6:268188/1‑149 (MQ=255)
                                                                                       tAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTg     >  5:379319/1‑149 (MQ=255)
                                                                                          ggACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTg  >  2:410228/1‑149 (MQ=255)
                                                                                           gACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTg  >  8:33380/1‑148 (MQ=255)
                                                                                           |                                                                                                                                                   
GTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGCTTTTCGTCAATAGTTAATTATCGACAGAGGTGTAATTGCTG  >  NC_000913/4606263‑4606501

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: ‑