Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,239 A→G 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,2390AG96.9% 104.1 / ‑5.3 32noncoding (48/87 nt)leuPtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/14);  minor base T (0/1);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  NC_000913/4606106‑4606368
                                                                                                                                     |                                                                                                                                  
cACGTCCGTAAGGACACAAACACCTGAAGCTAGCGCGTCTACCAATTCCGTCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                                               >  2:615868/1‑139 (MQ=255)
  cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGTCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                                             <  1:615868/139‑1 (MQ=255)
                           aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAccga                                                                                                    >  2:783696/1‑137 (MQ=255)
                           aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                                                    >  2:660214/1‑139 (MQ=255)
                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                                                         <  2:429968/108‑1 (MQ=255)
                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                                                         >  1:429968/1‑108 (MQ=255)
                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                                                 >  1:152452/1‑116 (MQ=255)
                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                                                 <  2:152452/116‑1 (MQ=255)
                                     tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                                                          >  1:506588/1‑139 (MQ=255)
                                           aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGTACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                                                    <  2:506588/139‑1 (MQ=255)
                                            aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                                                   <  2:94251/139‑1 (MQ=255)
                                                cGCCACCTTCGCACAGTTATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                                               <  1:660214/139‑1 (MQ=255)
                                                cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCGATTCCGCCACCTTCGCCCAGTg                                                                               <  1:783696/139‑1 (MQ=255)
                                                cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                                              <  1:341134/124‑1 (MQ=255)
                                                cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                                              >  2:341134/1‑124 (MQ=255)
                                                  ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                                             <  1:44915/139‑1 (MQ=255)
                                                    atcataaaacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                                           <  2:1030228/131‑1 (MQ=255)
                                                        tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                                                   <  1:94574/127‑1 (MQ=255)
                                                        tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                                                   >  2:94574/1‑127 (MQ=255)
                                                                tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAACAATAACAACGt                                                               >  2:748195/1‑139 (MQ=255)
                                                                   tCTTACTTTTTTTGATATCGCCTCGTTTGGTTCGAGGGGGGGGACCTGAACCCCCTCGTTCGTAATGACACTCACAGCTACAGCTTGCGCGTCTACCACTTCCGCCACCTTCGCCCAGTGCGAGCAATATGAACGTGGt                                                            >  2:450632/1‑139 (MQ=21)
                                                                                 atatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                                              >  2:1054000/1‑91 (MQ=255)
                                                                                 atatCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                                              <  1:1054000/91‑1 (MQ=255)
                                                                                                     aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                          >  1:317445/1‑139 (MQ=255)
                                                                                                           gggCCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCGAGCGCGTCTACCAGTTCCGCCACCTTCGCCCAGTGCGTGAAATATCAACCTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGt                    <  1:450632/139‑1 (MQ=21)
                                                                                                           gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGTCTTGAACCTCAACTTCCgg                    >  2:190038/1‑138 (MQ=18)
                                                                                                                  gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCa                           >  2:227568/1‑125 (MQ=255)
                                                                                                                  gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCa                           <  1:227568/125‑1 (MQ=255)
                                                                                                                  gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGt                    <  2:298656/132‑1 (MQ=255)
                                                                                                                  gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGt                    >  1:298656/1‑132 (MQ=255)
                                                                                                                     cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTa          >  2:986674/1‑139 (MQ=255)
                                                                                                                             tCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTGTTTGGTGCGAGGGGGGGGGCCTTGAACCCCCCCGTCCGGAAGGACACTaataacaga  >  1:487818/1‑132 (MQ=25)
                                                                                                                                     |                                                                                                                                  
CACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  NC_000913/4606106‑4606368

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