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_000913	4,606,322	1	.	G	95.0%	127.3 / 0.3	40	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (22/16);  total (23/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 = 4.26e-01

TTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606183‑4606449                                                                                                                                             |                                                                                                                                    tttGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                                                                                                                                      <  1:1021667/140‑1 (MQ=255)            cTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg‑gg‑g‑ggggac                                                                                                                      <  2:605238/140‑4 (MQ=255)                 ttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGaa                                                                                                                  <  2:650128/140‑1 (MQ=255)                  tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAAc                                                                                                                 <  1:143208/140‑1 (MQ=255)                        ggggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGG‑GG‑G‑GGGACTTGAACCCCCa                                                                                                            <  1:909286/138‑1 (MQ=255)                         aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGGACCCCCACg                                                                                                          >  1:963766/1‑140 (MQ=255)                         aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACg                                                                                                          <  2:197309/140‑1 (MQ=255)                         aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg‑gg‑g‑g                                                                                                                           <  2:427351/122‑1 (MQ=255)                         aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gg‑gg‑g‑g                                                                                                                           >  1:427351/1‑122 (MQ=255)                              ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGt                                                                                                     <  2:588612/140‑1 (MQ=255)                                    ttGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGac                                                                                               >  1:1113574/1‑140 (MQ=255)                                      gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCa                                                                                                            >  1:1091678/1‑125 (MQ=255)                                      gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCa                                                                                                            <  2:1091678/125‑1 (MQ=255)                                      gAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGacac                                                                                             >  1:1034510/1‑140 (MQ=255)                                       aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                            <  2:1034510/140‑1 (MQ=255)                                            ccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGAGCGA‑‑GG‑GGCG‑GGGTCTTGAACCCCCACGTCCGGAAAGTCACTAatac                                                                                       >  2:338966/1‑140 (MQ=25)                                                         aaCACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                                                          >  2:93136/1‑140 (MQ=255)                                                           cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGt                                                                        >  1:300686/1‑140 (MQ=255)                                                              tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTa                                                                     >  2:312130/1‑140 (MQ=255)                                                                acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                   >  1:173981/1‑140 (MQ=255)                                                                     tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAAGACACTAACATCTCAAGCTAGCGCGTCTACCAATTc                                                              >  2:728290/1‑140 (MQ=25)                                                                      gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGCGG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                             >  2:76254/1‑140 (MQ=255)                                                                          cTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGCTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                         >  2:987541/1‑140 (MQ=255)                                                                           tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                        <  1:312130/140‑1 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                       >  1:1043552/1‑140 (MQ=255)                                                                                cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                   <  1:987541/140‑1 (MQ=255)                                                                                 gTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGTCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                  >  2:806661/1‑140 (MQ=255)                                                                                  tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                 <  2:1043552/140‑1 (MQ=255)                                                                                     aCCAATTCCGCCACCTTCGCACAGTGCGAGCAATAACAACGTGGTTTTTGGAGCGA‑‑GG‑GG‑GAGGCCCTTGAACCCCAACTTCTGCAGGAAAACTAATACCAGAAGCTAGAGCGACTGCTAGTTGCGCGTCCTTCACAAAc                                              >  2:1019226/1‑140 (MQ=255)                                                                                         aTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                          <  2:1108999/124‑1 (MQ=255)                                                                                         aTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                          >  1:1108999/1‑124 (MQ=255)                                                                                               ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                   >  1:90222/1‑126 (MQ=255)                                                                                               ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCt                                    >  1:93269/1‑140 (MQ=255)                                                                                               ccACCTTCGCCCAGGGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCGCCACCTTCg                                                   <  2:90222/126‑1 (MQ=18)                                                                                                      cGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAg                             <  1:325274/140‑1 (MQ=255)                                                                                                         ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                          <  1:914015/140‑1 (MQ=255)                                                                                                             tGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCt                      >  2:365802/1‑140 (MQ=255)                                                                                                                 aGCAATATCAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCa                  >  1:473923/1‑140 (MQ=255)                                                                                                                         cAACGTGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg          <  2:473923/140‑1 (MQ=255)                                                                                                                              tGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc     >  1:531311/1‑140 (MQ=255)                                                                                                                              tGGTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc     >  2:700430/1‑140 (MQ=255)                                                                                                                                gTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa   >  2:915526/1‑140 (MQ=255)                                                                                                                                gTTTTTGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa   >  2:249135/1‑140 (MQ=255)                                                                                                                                 tttttGGTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  1:909020/140‑1 (MQ=255)                                                                                                                                 tttttGGTGCGA‑‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg               >  2:1008693/1‑126 (MQ=255)                                                                                                                                 tttttGGTGCGA‑‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg               <  1:1008693/126‑1 (MQ=255)                                                                                                                                      ggTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg             <  1:382589/124‑1 (MQ=255)                                                                                                                                      ggTGCGAG‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg             >  2:382589/1‑124 (MQ=255)                                                                                                                                             |                                                                                                                                    TTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GG‑GG‑G‑GGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606183‑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: ‑