Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	925,899	(G)7→8	noncoding (73/88 nt)	serW ←	tRNA‑Ser

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,892	1	.	G	97.1%	115.3 / ‑3.8	36	noncoding (80/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (3/31);  minor base T (0/1);  total (3/33)
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.85e-01

TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGT  >  NC_000913/925754‑926026                                                                                                                                             |                                                                                                                                         ttCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑g                                                                                                                                        <  1:535795/140‑1 (MQ=255)   cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑gtg                                                                                                                                      >  1:803569/1‑138 (MQ=255)   cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑ggg                                                                                                                                      >  1:767409/1‑140 (MQ=255)   cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑ggg                                                                                                                                      >  2:34135/1‑140 (MQ=255)   cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑ggg                                                                                                                                      >  2:646820/1‑140 (MQ=255)     cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑Gg‑‑ggggg                                                                                                                                    >  1:690383/1‑140 (MQ=255)     cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑ggggg                                                                                                                                    <  2:829427/140‑1 (MQ=255)      tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GC‑‑gggggg                                                                                                                                   >  2:683248/1‑140 (MQ=255)             cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCg                                                                                                                             <  1:1284506/140‑1 (MQ=255)               ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGaa                                                                                                                           <  2:1053779/140‑1 (MQ=255)               ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGaa                                                                                                                           <  2:624982/140‑1 (MQ=255)                       aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGAt                                                                                                                   <  1:494715/140‑1 (MQ=255)                        gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAGGGGGGGGGATTCGAACCCCCGAt                                                                                                                   <  2:767409/140‑1 (MQ=255)                          aaTTCACAGGTTACAAACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCCTCTTTAAATATGGCGGG‑GAG‑GGGGGGGATTCGAACCCCCGATACg                                                                                                                <  2:1309235/140‑1 (MQ=25)                          aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACg                                                                                                                <  1:333781/140‑1 (MQ=255)                          aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACg                                                                                                                <  1:702776/140‑1 (MQ=255)                                   gTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGtat                                                                                                       <  2:240585/140‑1 (MQ=255)                                      aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAT‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                    <  1:841690/140‑1 (MQ=255)                                      aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATACGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                    <  1:1210857/140‑1 (MQ=255)                                        aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATacaca                                                                                                  <  1:390129/140‑1 (MQ=255)                                           cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTTGGACCCCCCCTACGTTGCCCGATACACA‑Ctt                                                                                              >  2:885307/1‑140 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑‑GGGGGGGGGTTGAACCCCCCCTACGTTGCCGGCTAAACA‑Ccttt                                                                                            >  2:1294939/1‑139 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGA‑‑GGGGGGGGGTTGCACCCCCCCCACGTTGCCCGATTCACA‑CTTTc                                                                                            >  2:841690/1‑140 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTTTAAACCCCCCTACCTTGCCGTGTacacaacttc                                                                                             >  2:885315/1‑135 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTCGAACCCCCCCTACGTTGCCGTATACAAA‑CTTTc                                                                                            >  1:872448/1‑140 (MQ=21)                                             ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTc                                                                                            <  2:489789/140‑1 (MQ=255)                                              cACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTCGGACCCCCCCCACGTTGCCGTATACACA‑CTTTcc                                                                                           >  2:390129/1‑140 (MQ=21)                                                cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTCGAAACCCCCCGACCTTGCCGTATACAAACCTTTccg                                                                                          >  2:743671/1‑139 (MQ=21)                                                cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAg                                                                                         <  1:885315/140‑1 (MQ=255)                                                 ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGG‑GAG‑GGGGGGGGTTCGAACCCCCCCCACCTTTCCCTATACAAC‑CTTTCCAgg                                                                                        >  2:1090327/1‑140 (MQ=21)                                                     ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTg                                                                                    <  2:1172857/140‑1 (MQ=255)                                                     ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTg                                                                                    <  2:506130/140‑1 (MQ=255)                                                     cagccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTg                                                                                    <  1:709294/137‑1 (MQ=255)                                                               cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGcc                                                                         <  1:646820/140‑1 (MQ=255)                                                                    ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTGGGACCCCCCCTACGTTTCCGTGTTTACA‑CTTTTCAGGCGGGGTCCCTCCGCCCCTc                                                                     >  2:497915/1‑140 (MQ=12)                                                                       tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGG‑‑GGGGGGGGGTGGAAACCCCCCCACGTTTCCGGATACACA‑CTTTCCCGGCGGGGTCCCTTCGCCAccccgg                                                                  >  2:322515/1‑135 (MQ=12)                                                                       tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGGGTTGAAACCCCCCTACAGTGGCGTATACACA‑CATTTCACGCGCGCGCCCTCCGCCAccccgg                                                                  >  1:1026824/1‑135 (MQ=12)                                                                       tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGa                                                                  <  1:34135/140‑1 (MQ=255)                                                                         cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                <  1:534538/140‑1 (MQ=255)                                                                              agaTAAAATGCGAAAAAAAAGCCCGCGCTTTCGTACGAACTCTTCTTTAAATATGGAGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACT‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCa                                                           <  2:1001417/140‑1 (MQ=11)                                                                                                 aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCa                                        <  1:583438/140‑1 (MQ=255)                                                                                                  gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCaa                                       <  1:388760/140‑1 (MQ=255)                                                                                                       cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTc                                  <  2:589426/140‑1 (MQ=255)                                                                                                        aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                 <  1:743671/140‑1 (MQ=255)                                                                                                               tACGAGCTCTTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTgg                          <  1:497915/140‑1 (MQ=255)                                                                                                                     ctctTCTTTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgg                    <  1:1035738/140‑1 (MQ=255)                                                                                                                       cttcttTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgggg                  <  2:1026824/140‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg                <  1:322515/140‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg                <  1:1181169/140‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg                <  1:1178370/140‑1 (MQ=255)                                                                                                                         tcttTAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCACGGGTGGCAACGGGgcg                <  2:1234063/140‑1 (MQ=255)                                                                                                                            ttAAATATGGCGGT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGcta             <  1:683248/140‑1 (MQ=255)                                                                                                                                     gtggT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGa    >  2:450673/3‑140 (MQ=255)                                                                                                                                       ggT‑GAG‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCGAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt  <  1:865175/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                         TTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGT‑GA‑‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGT  >  NC_000913/925754‑926026

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