Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,043,483	(G)5→6	noncoding (75/90 nt)	serU ←	tRNA‑Ser

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,043,478	1	.	G	88.9%	75.9 / 4.1	27	noncoding (80/90 nt)	serU	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (11/13);  total (12/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.91e-01

TGAATGTCACAGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑A‑G‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCATTTTAACAGTCCCTG  >  NC_000913/2043352‑2043608                                                                                                                                   |                                                                                                                                       tGAATGTCACAGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGaa                                                                                                                         >  1:520962/1‑140 (MQ=255)   aaTGTCACAGTATCGTACAGAAAACGGTGACGATCTAGCCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAc‑c                                                                                                                      <  1:249842/140‑1 (MQ=255)   aaTGTCACAGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAc‑c                                                                                                                      <  1:562981/140‑1 (MQ=255)           aGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAc                                                                                                                        <  2:405676/131‑1 (MQ=255)           aGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAc                                                                                                                        >  1:405676/1‑131 (MQ=255)            gTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTaga                                                                                                             <  1:472439/140‑1 (MQ=255)               tCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGAGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGtt                                                                                                          >  2:93940/1‑140 (MQ=255)                 gAACAGAAAACAGTGACGATCTGACGCTTCAAGAATAATCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGc                                                                                                        <  1:311674/140‑1 (MQ=255)                  aaCAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGcc                                                                                                       >  1:531217/1‑140 (MQ=255)                    cAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGAGGTG‑AGG‑A‑G‑G‑G‑GGGTTTGAACGCCCCGGT‑GAGTTGcccc                                                                                                     >  2:317811/1‑140 (MQ=255)                       aaaaCAGTGACGATCTAACCCGTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAATTAAGTCATTGATAAAGTGGCGGAGAA‑G‑G‑G‑C‑G‑GGGTTTGAAA‑CCCCGGTGAAGTTTTCCCTAc                                                                                                  >  2:415228/1‑140 (MQ=255)                                 cGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGGTTTGAAA‑CCCCGGTAGAGTTGCCCCTACTCCGGttttt                                                                                        >  2:548228/1‑139 (MQ=255)                                     cTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCgaga                                                                                    <  1:460285/140‑1 (MQ=255)                                             tCAAGAATATTCTAAGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGGTTTGAAA‑CCCCGGTAGAGTTGCCCCTACTCCGGGTTTCGCGACCGGTCCg                                                                            >  2:559062/1‑140 (MQ=255)                                              cAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTAGAAC‑CCCCGGTAAAGTTGCCCCTACGCCCGGTTTCGAGACCGGTCCGt                                                                           >  2:178447/1‑140 (MQ=255)                                                  aaTATTCTAAGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAg                                                                       <  1:559062/140‑1 (MQ=255)                                                       tCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGGCGGGCGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGGGTTGAAC‑CCCCGGGAGAGATGGCCCTGCCCCCGGTTTCGAGAGCTGTCCGGTCaagccgg                                                                  >  2:394723/1‑134 (MQ=255)                                                         tACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑GAA‑G‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTcc                                                                >  1:329531/1‑140 (MQ=255)                                                            gATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGc                                                             >  1:416604/1‑140 (MQ=255)                                                                gTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCAtct                                                         <  1:334640/140‑1 (MQ=255)                                                                      tttAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGGTTTCGAGACCCGTCCGTTCAGCCGCTCCGGCATCTCTccct                                                   >  2:362530/1‑138 (MQ=255)                                                                               aGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGt                                          <  2:416604/140‑1 (MQ=255)                                                                                            gTCGGGAGATAAGTCATTGATACAGAGGCAGAG‑A‑G‑A‑GCG‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCa                             <  1:317811/140‑1 (MQ=255)                                                                                              cGGGAGATAAGTCATCGATAAATTGTCGGCG‑A‑G‑A‑G‑G‑GCGGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAgg                           <  1:176094/140‑1 (MQ=255)                                                                                                    aTAAGTCATCGATAAAGTGGCGGAG‑A‑G‑A‑G‑GTG‑GGATTTGAAC‑CCACGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAAttt                     <  1:415228/140‑1 (MQ=255)                                                                                                      aaGTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTgg                   <  1:362530/140‑1 (MQ=255)                                                                                                        gTCATTGATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCa                 <  1:548228/140‑1 (MQ=255)                                                                                                        gTCATTGATAAAGTGGCAGAG‑A‑C‑AGG‑G‑G‑GGATTTGGAC‑CGCCGGTAGATTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCa                 <  2:229249/140‑1 (MQ=255)                                                                                                              gATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCATTTTaa           <  2:531217/140‑1 (MQ=255)                                                                                                              gATAAAGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCATTTTaa           >  1:49621/1‑140 (MQ=255)                                                                                                                   aGTGGCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCATTTTAACAGTc      >  2:560672/1‑140 (MQ=255)                                                                                                                       gCGGAG‑A‑G‑AGG‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTATCCGTTCAGATGGTCGCCGTGAGGCCAGGAACTTTGCCATTTTAACAGTCCTTg  >  2:257144/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                       TGAATGTCACAGTATCGAACAGAAAACAGTGACGATCTAACCCTTCAAGAATATTCTACGATTGTTCTGTTTAGGAAAAGCAAGGCGGGAAGTCGGGAGATAAGTCATTGATAAAGTGGCGGAG‑A‑G‑A‑G‑G‑G‑GGATTTGAAC‑CCCCGGTAGAGTTGCCCCTACTCCGGTTTTCGAGACCGGTCCGTTCAGCCGCTCCGGCATCTCTCCGTTCAGATGGTTGCCATGATGCCAGGAAATTTGGCATTTTAACAGTCCCTG  >  NC_000913/2043352‑2043608

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: ‑