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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,2071.G90.0% 82.9 / 3.6 30noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (13/14);  total (15/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 = 9.93e-01
*NC_0009134,606,2072.G86.7% 80.4 / 7.2 30noncoding (80/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (12/14);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑G‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGC  >  NC_000913/4606083‑4606320
                                                                                                                               ||                                                                                                                        
gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑g‑‑ggggg‑‑ggggacttg                                                                                                      >  1:818494/1‑133 (MQ=255)
        ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGG‑‑GGGGG‑‑GGACTTGAACCCCCAc                                                                                             >  1:653909/1‑140 (MQ=255)
         gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTTCGCCACCTTCGCACCGTCATCTTACTTTTTTTGATATCCGCTTGTTTGGTGC‑‑GA‑‑GCGGGGAG‑‑GGACTTGAACCCCCAc                                                                                             >  2:373448/1‑140 (MQ=18)
         gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCAc                                                                                             >  2:520174/1‑140 (MQ=255)
               tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGGGCGTCTACCAAGTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGt                                                                                       <  2:172828/140‑1 (MQ=25)
               ngAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTAGCACAGTCATCTTACTTTTTTTGATATCGCCTCATTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGt                                                                                       <  2:499206/139‑1 (MQ=25)
                  aCCCCCACGTCCGTAAGGACACTAATACCCGAAGCTATCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAg                                                                                    <  1:373448/140‑1 (MQ=37)
                  aCCCCCACGTCCGTAAGCACACTAACACCTGGAGCGAGCACGTCTACCCATTCCGCCGCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAg                                                                                    <  2:302467/140‑1 (MQ=25)
                      ccAGGTCCGTAAGGAAACCAACACCGCAAGCTAACACAACTACCAAGTCCGCCACCTTAACACAGTAATATTACTTATTTTGAGATGGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAaca                                                                                <  2:523031/140‑1 (MQ=19)
                        aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACt                                                                              <  1:520174/140‑1 (MQ=255)
                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGATCCTACTTTTTTTGTTATCTCTGCGGTTGGTGC‑‑GG‑‑G‑‑GCGCG‑‑GGAGCCTGAAATCCCCACGGCCATGAGAACACTAACACCGCAAGcga                                                                >  2:348777/1‑138 (MQ=12)
                                           ccccTGAAGCTAGCGCGTCTACAAATTCCCCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGGG‑‑G‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                           <  1:296470/138‑1 (MQ=25)
                                             ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGGGG‑‑G‑‑GGGGG‑‑GGAGTTGAACCCGCACGTCCATAAGAACATTAACACCTGAAGCTAGCGCGTc                                                         >  2:296470/1‑140 (MQ=37)
                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                                      >  1:554593/1‑140 (MQ=255)
                                                        cgcgTCTACCAATTCCGCCACCTTCGCCCAGTCCTCTCACTTTTTTTGATATCCCCCCGACGCGGGC‑‑GA‑‑G‑‑GGGCGACGGTGTGTAACCCCCGCGTCCGGACGAACCCTTAAACCAGAAGCTAGCGCGTCTACCAATTCCg                                              >  2:745349/1‑140 (MQ=255)
                                                         gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGG‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                             <  2:554593/140‑1 (MQ=25)
                                                                  aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                    >  1:653660/1‑140 (MQ=255)
                                                                  aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCATCGTTTGGTGC‑‑GAGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   <  2:653909/140‑1 (MQ=255)
                                                                   aTTCCGCCACCTTCGCACAGTCATCTTTCTTTTTTTGATATGGCCTCGGTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAGCCCCCGCGTCCGTAAGAACGCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   <  1:745349/140‑1 (MQ=21)
                                                                     tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGCCTACCAATTCCGCCACCTTCGCCCa                                 >  1:150185/1‑140 (MQ=255)
                                                                       cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                               <  2:653660/140‑1 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   <  1:440553/134‑1 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   >  2:440553/1‑134 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                               >  1:510802/1‑138 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                               <  2:510802/138‑1 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                             >  2:723828/1‑140 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                             >  2:208678/1‑140 (MQ=255)
                                                                         ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑GGCGGGGG‑‑GGACTTGGACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCACGTCAACCACTTCCGCCACCTCCGCCCaccgc                             >  1:503757/1‑136 (MQ=21)
                                                                            ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                          <  1:593240/140‑1 (MQ=255)
                                                                              ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTTAACCCCCACGTCCCTAAAAACACTAACACCCGAAGCTATGTTGTCCACCAAATCCGCCCCCTTCCCCCCGTGGGGGCa                        >  2:642008/1‑140 (MQ=255)
                                                                                cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   <  1:771796/127‑1 (MQ=255)
                                                                                cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                   >  2:771796/1‑127 (MQ=255)
                                                                                  aacaGTCATCTTACTTCTTTTGGTAACGCCTAGTTCGGTCC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat                    <  2:395587/139‑1 (MQ=18)
                                                                                                 ttttttGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg     >  1:434966/1‑140 (MQ=255)
                                                                                                    tttGATATCGCCTCGTTTGGTGC‑‑GAGGG‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGTGTTTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc  >  2:127067/1‑140 (MQ=21)
                                                                                                                               ||                                                                                                                        
GCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑‑GA‑‑G‑‑GGGGG‑‑GGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGC  >  NC_000913/4606083‑4606320

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