Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,520,976	(G)5→6	noncoding (46/76 nt)	valU →	tRNA‑Val

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,520,971	1	.	G	97.8%	149.1 / ‑3.5	45	noncoding (41/76 nt)	valU	tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (17/27);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGATCATCTTGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑A‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGT  >  NC_000913/2520833‑2521097                                                                                                                                             |                                                                                                                               gCGATCATCTTGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg‑a‑g                                                                                                                               >  1:836195/1‑140 (MQ=255)       atcTTGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑A‑ggggggt                                                                                                                         >  2:606304/1‑138 (MQ=255)         cTTGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑A‑ggggggtcg                                                                                                                       >  1:212434/1‑136 (MQ=255)          ttGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑A‑ggggggtcgg                                                                                                                      <  1:541537/140‑6 (MQ=255)                      acgacgAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGACGGCGGTTCGATccc                                                                                                           >  2:247338/1‑140 (MQ=255)                            aacgaacAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCAACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGACCACCTCCCTTACAAGG‑ACGGGGGTCGGCGGTTCGGTCCCGtcatc                                                                                                     <  2:786002/140‑1 (MQ=37)                              cgaacAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATcac                                                                                                   <  1:606304/140‑1 (MQ=255)                                aacAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATcaccc                                                                                                 <  1:941831/140‑1 (MQ=255)                                aacAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATcaccc                                                                                                 <  2:691367/140‑1 (MQ=255)                                     tttCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCaa                                                                                            <  1:317215/140‑1 (MQ=255)                                     tttCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCaa                                                                                            <  1:756202/140‑1 (MQ=255)                                      ttCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCaa                                                                                            >  1:627172/1‑139 (MQ=255)                                      ttCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCaa                                                                                            <  2:627172/139‑1 (MQ=255)                                       tCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCTGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAact                                                                                          <  1:720207/140‑1 (MQ=255)                                          ttAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGGTCCCGTCATCACCCACCAactact                                                                                       <  2:836195/140‑1 (MQ=255)                                                               tttGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTa                                                                  >  1:904522/1‑140 (MQ=255)                                                                    aCTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAg                                                             <  2:262741/140‑1 (MQ=255)                                                                       ctcCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGCCATCACCCACCAACTACTTTCTGTAGTCTCCGCCGTGTAGCAAGaaa                                                          >  1:633097/1‑140 (MQ=255)                                                                        tcCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAAt                                                         >  2:366326/1‑140 (MQ=255)                                                                                  gCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATTCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTgg                                               <  2:316412/140‑1 (MQ=255)                                                                                  gCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTgg                                               >  1:608690/1‑140 (MQ=255)                                                                                  gCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGAgg                                               >  1:145768/1‑140 (MQ=255)                                                                                  gCGACTCCACACAGCGGGGGTGAGTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCATCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTgg                                               <  1:462177/140‑1 (MQ=255)                                                                                    gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTTGATCGCGTCATTACCCCCCAAACACTTTATTTAGTCTCCCCCGGGTAGCACGGAATTGAGAAGAGGGt                                             >  2:462177/1‑140 (MQ=18)                                                                                    gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTgg                                               >  2:130453/1‑138 (MQ=255)                                                                                    gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTgg                                               <  1:130453/138‑1 (MQ=255)                                                                                    gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGt                                             >  1:296689/1‑140 (MQ=255)                                                                                    gACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAACAAGAAATTGACaaaggggg                                             >  2:867025/1‑134 (MQ=21)                                                                                      cTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGCG‑GGGGGTCGGCGGCTCCATCCCCTCCTCACCCACCAACTACTTTATGTAGTCTCCGCCGTGGAGCAAGAAATTGAGAAGgggggca                                           >  2:644539/1‑137 (MQ=25)                                                                                          acacaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAg                                       <  2:741447/140‑1 (MQ=255)                                                                                            acaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCt                                     <  1:914182/140‑1 (MQ=255)                                                                                            acaGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCt                                     <  2:608690/140‑1 (MQ=255)                                                                                             caGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTc                                    <  2:212434/140‑1 (MQ=255)                                                                                              aGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCa                                   >  2:437288/1‑140 (MQ=255)                                                                                              aGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCa                                   <  2:904522/140‑1 (MQ=255)                                                                                                ggggggTGATTAGCTCAGCCTGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCCTCATCACCCACCAACTTCTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGc                                 <  2:487028/139‑1 (MQ=25)                                                                                                cGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGt                                                 <  1:700711/124‑1 (MQ=255)                                                                                                cGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGt                                                 >  2:700711/1‑124 (MQ=255)                                                                                                cGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGc                                 >  1:802191/1‑140 (MQ=255)                                                                                                cGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGCCGGCGGTTCGATCCCGCCATCACCCACCAACTAGGTTCTGTATTTTTTCCCGTGTAGCAAGGAATTCAGAAGAGGGGGATTAGATCAGc                                 >  2:586782/1‑140 (MQ=17)                                                                                                 gggggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCt                                <  2:296689/140‑1 (MQ=255)                                                                                                     gTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGga                            <  1:437288/140‑1 (MQ=255)                                                                                                        aTTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGagag                         <  2:802191/140‑1 (MQ=255)                                                                                                           aGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAc                      <  2:402697/140‑1 (MQ=255)                                                                                                            gCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc                     <  2:933353/140‑1 (MQ=255)                                                                                                            gCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCAcc                     <  1:501921/140‑1 (MQ=255)                                                                                                             cTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                    <  1:644539/140‑1 (MQ=255)                                                                                                             cTCAGCTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCt                    <  1:366326/140‑1 (MQ=255)                                                                                                                  cTGGGAGAGCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCCACTACTTTATGGAGTCTCCGCCGGGTAGCAAGGAATTGGGAAGTGGGTGATTAACTCAGCTGGGAGAGCACCTCCCtt               >  2:501921/1‑140 (MQ=25)                                                                                                                       agagCACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGGGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAag          >  1:873707/1‑140 (MQ=255)                                                                                                                           cACCTCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTCTGTACTCTTCCCCGTGAAACAAGAAATTGACAAAAGGGGGGTTAGCTTCTCTTGGAGAACACCTCCCTTAAAaggagg      >  2:76444/1‑140 (MQ=16)                                                                                                                               tCCCTTACAAGG‑AGGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGt  >  1:666278/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                               GCGATCATCTTGCTTAAATTTACGACGAACGAACAATTTCTTTAGAAAATGCGTTGACTCATTTTGAACTCTCCCTATAATGCGACTCCACACAGCGGGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGG‑A‑GGGGGTCGGCGGTTCGATCCCGTCATCACCCACCAACTACTTTATGTAGTCTCCGCCGTGTAGCAAGAAATTGAGAAGTGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGT  >  NC_000913/2520833‑2521097

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