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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,207	1	.	G	90.9%	127.8 / 5.8	44	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (19/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTA  >  NC_000913/4606057‑4606351                                                                                                                                                          |                                                                                                                                                   cAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGa                                                                                                                                                     <  2:37733/151‑1 (MQ=255)      aaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑ggg                                                                                                                                             <  1:1007589/151‑1 (MQ=255)      aaGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑ggg                                                                                                                                             <  1:861496/151‑1 (MQ=255)        ggAAGTACAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑ggggg                                                                                                                                           >  1:52375/1‑151 (MQ=255)         gAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑gggggg                                                                                                                                          >  1:56265/1‑151 (MQ=255)          aaGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑ggggggg                                                                                                                                         >  1:872426/1‑150 (MQ=255)              aaaGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑gg‑gggggggact                                                                                                                                      <  2:212927/151‑5 (MQ=255)               aaGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTg                                                                                                                                     >  1:974085/1‑151 (MQ=255)                       tGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACACTCACCTTACTTTTAATGATTTCCTCTCGTGGAGGGGCAT‑‑GG‑GGGGGGACTTCAAc                                                                                                                                  >  2:370520/1‑146 (MQ=14)                        ggTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCa                                                                                                                             <  1:598884/150‑1 (MQ=255)                            cGA‑GGGGGGAGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTCTTACTTTTATTGATATCGCTTCGT‑TTGGTGCGA‑‑GGCTGGGGGACTTGAACCCCCCCGTcc                                                                                                                        >  2:641365/1‑151 (MQ=21)                              a‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGt                                                                                                                      >  1:187063/1‑151 (MQ=255)                                   gggggACTTGAACCCCCACGTCCGTAAGGACACTAACGCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGa                                                                                                                  <  2:775281/151‑1 (MQ=255)                                       gACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                              >  1:417931/1‑151 (MQ=255)                                          ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGATATCTCCTCGT‑TTCGTGCGAG‑GG‑GGGGGGACTTTAACCCCGACGTCCGTACGAAGACTa                                                                                                            >  2:358386/1‑150 (MQ=21)                                           tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                           <  1:648072/150‑1 (MQ=255)                                            gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAac                                                                                                          <  1:1055907/150‑1 (MQ=255)                                               cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTATTGATATCGTCTCGT‑TGGGTGCGA‑‑GGCGGGGGGACTTCAACTCCTACTTCCGAAATaattctcacccgt                                                                                                      >  2:580374/1‑140 (MQ=18)                                                ccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                                                                    <  1:1031074/151‑1 (MQ=255)                                                  ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                                                                   <  1:85474/151‑1 (MQ=255)                                                   cacctacGTAAGGAAACTACCACCTGAAGCTAGCGTGACTACCAATTCCGCCACCTGCGCACAGTCATCTTACTTTTTGTGATATCGCCTCGT‑TTGGTGCGAGGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCGGaa                                                                                                   <  2:747763/145‑1 (MQ=37)                                                            aaGGACACTAACACCTGAAGCTAGCGAGTCTACCAATTCCGCCGCCGACGAACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                                          <  2:500805/150‑1 (MQ=37)                                                                     aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                                 <  1:894263/150‑1 (MQ=255)                                                                           tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                                           >  1:820833/1‑150 (MQ=255)                                                                           tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                                           >  1:85680/1‑150 (MQ=255)                                                                                      cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                                 >  2:227580/1‑149 (MQ=255)                                                                                       gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                              <  1:803076/151‑1 (MQ=255)                                                                                          tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTCAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCCCCCa                                                                >  2:176826/1‑146 (MQ=255)                                                                                              aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                         <  2:650465/149‑1 (MQ=255)                                                                                              aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                         >  1:650483/1‑149 (MQ=255)                                                                                               aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                                                      <  1:1060247/151‑1 (MQ=255)                                                                                                    gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                 <  1:226666/151‑1 (MQ=255)                                                                                                        ccTCCGCACGGTCATCGTACCTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACGAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCTCCTTCGCCCAGTGCGAGCAATATCAAc                                               <  2:1066008/149‑1 (MQ=21)                                                                                                         cgtagaacaGTCAACTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                 <  2:353009/140‑1 (MQ=255)                                                                                                           tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                          >  1:830953/1‑151 (MQ=255)                                                                                                           tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                          >  1:667939/1‑151 (MQ=255)                                                                                                              cacaGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                                       <  2:820802/151‑1 (MQ=255)                                                                                                               acaGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                      >  2:462446/1‑151 (MQ=255)                                                                                                               acaGTCATCTTACTTTTTTTGATATCGCATCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                              <  2:542268/143‑1 (MQ=255)                                                                                                                caGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                     >  1:728028/1‑151 (MQ=255)                                                                                                                  gTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTgg                                    <  2:1029500/150‑1 (MQ=255)                                                                                                                    cATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTg                                  >  2:469137/1‑150 (MQ=255)                                                                                                                         tACTTTTTTTGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                <  1:30686/147‑1 (MQ=255)                                                                                                                         tACTTTTTTTGATATCACCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgg                             <  2:509037/150‑1 (MQ=255)                                                                                                                          aCTATTTTTGATATCGCCTCGT‑TTGGTTGGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACTCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa                               <  2:316370/147‑1 (MQ=21)                                                                                                                                tttGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGAc                     >  1:84988/1‑151 (MQ=255)                                                                                                                                 ttGATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACTGCGAGCAATATCAACGAGGTTTTTGGTGCGAGGCGGGGGAc                     >  2:53830/1‑150 (MQ=255)                                                                                                                                   gATATCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTg                  >  1:220824/1‑151 (MQ=255)                                                                                                                                     tatCGCCTCGT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGTGCGAGGGGGGGGACtt                   >  2:405046/1‑148 (MQ=255)                                                                                                                                              gT‑TTGGTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTCCGAGGCGCGGGACTTGAACCCCCACGt       >  2:426076/1‑151 (MQ=255)                                                                                                                                                    gTGCGAG‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTa  <  2:6452/151‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                   CAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGT‑TTGGTGCGA‑‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTA  >  NC_000913/4606057‑4606351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑