Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,330 (G)8→9 noncoding (72/87 nt) leuQ ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,3221.G86.1% 99.3 / 11.0 36noncoding (80/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (16/15);  total (19/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 = 9.23e-01

TTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTC  >  NC_000913/4606184‑4606446
                                                                                                                                           |                                                                                                                             
ttGATATCGCCTCGTTTGGTGGGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                                                                                                                <  1:584977/139‑1 (MQ=255)
ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCTAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                                                                                                                >  1:713687/1‑139 (MQ=255)
ttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                                                                                                                >  1:137611/1‑139 (MQ=255)
              tttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggggggggcttg                                                                                                                >  1:155310/1‑133 (MQ=255)
                             gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                  >  2:520528/1‑139 (MQ=255)
                             gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGCGGGGTCTTGAACCCCCACGTCCGt                                                                                                  >  1:216414/1‑139 (MQ=255)
                                   tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGTGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                            <  2:713687/139‑1 (MQ=255)
                                    gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAAcccc                                                                                                           <  1:179664/123‑1 (MQ=255)
                                    gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAAcccc                                                                                                           >  2:179664/1‑123 (MQ=255)
                                    gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggggac                                                                                                                   >  1:765370/1‑111 (MQ=255)
                                    gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggggac                                                                                                                   <  2:765370/114‑4 (MQ=255)
                                    gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAACAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGTCTTGAACCCCCACGTCCGTAAGTaca                                                                                           >  2:926571/1‑139 (MQ=21)
                                      aCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                         <  1:520528/139‑1 (MQ=255)
                                          ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                     >  1:711668/1‑139 (MQ=255)
                                                       gacacTAACTCCTGAAGCTAGCGCCTCTACCAATTCCGCAACCGACGCCCAGTGCGATTATTAGAAACGGGGTTTTTGGGGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                                                        <  2:137611/138‑1 (MQ=18)
                                                          acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCa                                                                                                         >  1:797768/1‑103 (MQ=255)
                                                          acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCa                                                                                                         <  2:797768/103‑1 (MQ=255)
                                                                 ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                             >  2:86770/1‑139 (MQ=255)
                                                                      aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                         <  1:548432/139‑1 (MQ=255)
                                                                            cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                   >  2:739844/1‑139 (MQ=255)
                                                                                tCTACCAATTCCGCCACCTTCGCCCAGTGCGATCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGTCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                               >  2:548432/1‑139 (MQ=255)
                                                                                tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                     >  1:946561/1‑101 (MQ=255)
                                                                                tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                                                     <  2:946561/101‑1 (MQ=255)
                                                                                   tcctaTTCCGCCACGTTCGCTCAGTGCGAGAAACATCAACGGCGTTGTTGGGGGGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                            <  2:155310/135‑1 (MQ=14)
                                                                                             ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTc                                  >  1:235304/1‑139 (MQ=255)
                                                                                              cACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                               >  1:204193/1‑125 (MQ=255)
                                                                                              cACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                               <  2:204193/125‑1 (MQ=255)
                                                                                                           tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACg                                                                                                       >  2:957937/1‑56 (MQ=25)
                                                                                                           tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACg                                                                                                       <  1:957937/56‑1 (MQ=25)
                                                                                                               aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACg                                                                                                       <  2:526113/52‑1 (MQ=25)
                                                                                                               aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACg                                                                                                       >  1:526113/1‑52 (MQ=25)
                                                                                                               aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                     <  2:1088970/118‑1 (MQ=255)
                                                                                                               aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                     >  1:1088970/1‑118 (MQ=255)
                                                                                                                  aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg             <  1:739844/139‑1 (MQ=255)
                                                                                                                     atCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg          >  2:1069315/1‑139 (MQ=255)
                                                                                                                        aaCGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg       <  2:235304/139‑1 (MQ=255)
                                                                                                                           gTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                          >  2:504619/1‑69 (MQ=255)
                                                                                                                           gTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                          <  1:504619/69‑1 (MQ=255)
                                                                                                                            tGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc  <  1:86770/139‑1 (MQ=255)
                                                                                                                              gTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                     <  1:949784/103‑1 (MQ=255)
                                                                                                                              gTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                     >  2:949784/1‑103 (MQ=255)
                                                                                                                              gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                       >  1:83482/1‑116 (MQ=255)
                                                                                                                              gTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                       <  2:83482/116‑1 (MQ=255)
                                                                                                                                           |                                                                                                                             
TTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTC  >  NC_000913/4606184‑4606446

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: