Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,330 (G)8→10 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.G96.8% 94.9 / ‑3.1 31noncoding (80/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (18/12);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01
*NC_0009134,606,3222.G83.9% 78.0 / 11.7 31noncoding (80/87 nt)leuQtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (15/11);  total (18/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 = 4.17e-01

ATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606189‑4606449
                                                                                                                                         ||                                                                                                                                 
atCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCCAACAATATCAACGAGGTTTTTGGAGCGA‑‑gggg                                                                                                                              >  2:160282/1‑140 (MQ=18)
   gCCTCGTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggggggg                                                                                                                         <  1:543871/140‑1 (MQ=255)
      tCGGTTGGTGCGA‑‑GGGGTGGGGGATTGAACCCCCCCGTCCGTAAGAACACTAACCCCTTAAGCCAGCGCGTCTACCAATTCCGCCACCCTCGCCCAGTGCGAGTAATATCAACGTGGTTTTTGGTgggg‑‑gggggggggga                                                                                                                      <  2:287510/140‑4 (MQ=18)
           tGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGaa                                                                                                                  <  2:644673/140‑1 (MQ=255)
                 ggggggggggggCCTTGTACCCTCACGTCCGTAAGTACATTAACACCTTAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAAcccc                                                                                                              <  2:383378/136‑1 (MQ=21)
                 gA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACg                                                                                                          <  1:657169/140‑1 (MQ=255)
                   agggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCa                                                                                                            <  1:702359/138‑1 (MQ=255)
                              aCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                  <  1:86330/140‑1 (MQ=255)
                                  gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGaca                                                                                              >  1:572079/1‑140 (MQ=255)
                                   aaCCCCCACGTCCGTAAGAACACTAGCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCACCGTGGTTTTTGGCTCGAGCGGGGGGGGAATTGAACCCCCACGTCCGGAAGGacac                                                                                             >  2:691349/1‑140 (MQ=21)
                                        ccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCCAGCAATAACAACGTGGTTCTTGTTGCGA‑‑GGGCGGTGGTAGTTGAACCCACCCGTCCGTAAGTACACTAtca                                                                                        >  2:396587/1‑140 (MQ=18)
                                           cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATGAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                                                     <  1:244311/140‑1 (MQ=255)
                                                gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCGACGGCCTTAAGGAAACTAACACCTGAAGc                                                                                >  1:290877/1‑140 (MQ=37)
                                                gTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                                                                >  1:106252/1‑140 (MQ=255)
                                                                 tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                                              >  1:708050/1‑140 (MQ=255)
                                                                   aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                             >  1:888708/1‑140 (MQ=255)
                                                                    aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCCCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTccc                                                            >  1:15663/1‑139 (MQ=25)
                                                                       tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                         <  1:64826/140‑1 (MQ=255)
                                                                        aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                        <  2:888708/140‑1 (MQ=255)
                                                                          cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGTGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                      >  2:987902/1‑140 (MQ=255)
                                                                          cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                      >  2:541866/1‑140 (MQ=255)
                                                                              tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                  >  2:136735/1‑140 (MQ=255)
                                                                               cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                 >  1:959189/1‑140 (MQ=255)
                                                                                tACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAGCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCAt                                                <  2:106252/140‑1 (MQ=255)
                                                                                 aCCAATTCCGCCCCCTTCGCCCAGTGGGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                               <  1:691349/140‑1 (MQ=21)
                                                                                                  cGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                              >  2:720480/1‑140 (MQ=255)
                                                                                                     ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGaa                           <  2:572079/140‑1 (MQ=255)
                                                                                                        gTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTg                        >  1:218189/1‑140 (MQ=255)
                                                                                                         tGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTAGCATACCATCAATTCTTAAAAAGAATTGc                       >  1:453286/1‑140 (MQ=21)
                                                                                                             aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAcc                   >  1:760672/1‑140 (MQ=255)
                                                                                                                   atCAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg             >  1:182699/1‑140 (MQ=255)
                                                                                                                   atCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg            >  2:77610/1‑140 (MQ=255)
                                                                                                                     cAACGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgc           >  2:73625/1‑140 (MQ=255)
                                                                                                                      aaCGTGGTTTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg          <  2:959189/140‑1 (MQ=255)
                                                                                                                         gTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCCTCGCATACTATCAATTATTAAAAAGGATTGATACCACGGAGGCGCAtt      >  2:451030/1‑140 (MQ=18)
                                                                                                                              ttttGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  2:182699/140‑1 (MQ=255)
                                                                                                                              ttttGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  1:160282/140‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                                 
ATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606189‑4606449

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