Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,032,303	(C)5→6	coding (937/1449 nt)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,032,298	1	.	C	97.6%	150.4 / ‑3.5	42	coding (932/1449 nt)	uacT	uric acid permease
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (21/20);  total (22/20)
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.99e-01

TGCCGTGTTAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑CCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCC  >  NC_000913/3032149‑3032439                                                                                                                                                       |                                                                                                                                               tGCCGTGTTAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑c                                                                                                                                               >  2:659178/1‑151 (MQ=255)     gtgtTAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAACCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑ccc‑cc                                                                                                                                          >  1:382484/1‑151 (MQ=255)       gtTAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑ccc‑ccca                                                                                                                                        <  2:318562/151‑3 (MQ=255)         tAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑ccc‑cccac                                                                                                                                       >  1:471670/1‑147 (MQ=255)                         cGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCaaa                                                                                                                        >  1:315299/1‑150 (MQ=255)                          gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAAc                                                                                                                      <  1:134012/151‑1 (MQ=255)                          gAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAAc                                                                                                                      >  1:229166/1‑151 (MQ=255)                               caTGGGGAAGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTgg                                                                                                                 >  1:345201/2‑151 (MQ=255)                               aaTGGGGATGTTCCTGGCACTGGGTGAAATAGTCCGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGCGTCGATGTCGGAGGTACAATGATAGGAGACACGCTTAAAAGCTT‑CCCAACACGCGTCAATTACTCacaaa                                                                                                                      >  1:500210/1‑145 (MQ=255)                                aTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGc                                                                                                                >  2:530306/1‑151 (MQ=255)                                aTGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGc                                                                                                                <  2:516554/151‑1 (MQ=255)                                      aTGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGtt                                                                                                          >  1:237017/1‑151 (MQ=255)                                        gTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAg                                                                                                        <  2:501019/151‑1 (MQ=255)                                        gTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAg                                                                                                        <  2:321018/151‑1 (MQ=255)                                              ggCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTg                                                                                                    <  2:53334/149‑1 (MQ=255)                                               gCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAcg                                                                                                 >  2:121266/1‑151 (MQ=255)                                                  cTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGAcgcgc                                                                                              >  2:655178/1‑151 (MQ=255)                                                         aaaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCAt                                                                                        >  1:97663/1‑150 (MQ=255)                                                          aaTAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAg                                                                                      <  1:442154/151‑1 (MQ=255)                                                              gtcggtcgTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGc                                                                                  >  2:21008/1‑151 (MQ=255)                                                                     gTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgt                                                                           <  1:122881/151‑1 (MQ=255)                                                                      tAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGCCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                          >  2:144843/1‑151 (MQ=255)                                                                       aaaCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGgtgtg                                                                          <  1:678769/150‑1 (MQ=255)                                                                         aCTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAt                                                                       <  2:465033/151‑1 (MQ=255)                                                                             tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                   >  1:499344/1‑151 (MQ=255)                                                                             tctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttct                                                                   >  2:375623/1‑151 (MQ=255)                                                                              ctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTAttt                                                                     <  1:488499/148‑1 (MQ=255)                                                                              ctTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATttctt                                                                  >  2:605411/1‑151 (MQ=255)                                                                                 cGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCgg                                                               <  2:256408/151‑1 (MQ=255)                                                                                    aCGATATTATTTGCGGGCTGCGTGTCGCTGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGGGTGTATTTCTTCggg                                                              >  2:154496/1‑149 (MQ=255)                                                                                         attattCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCACAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATtata                                                       <  2:173023/151‑1 (MQ=255)                                                                                             ttCGCGGGCTGCGTGTCGATGGCGTAGTGACAATGATAGGCGGCTCGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCTCTGGGTGTGTATTTCTTCGGGAATTATAtt                                                     <  2:504862/149‑1 (MQ=255)                                                                                               cgcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAAt                                                  <  1:166627/150‑1 (MQ=255)                                                                                               cgcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATc                                                 >  2:58629/1‑151 (MQ=255)                                                                                               cgcgGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATc                                                 >  1:147051/1‑151 (MQ=255)                                                                                                      tGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTg                                          >  1:359211/1‑151 (MQ=255)                                                                                                       gCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGAAATTATATTAATCCTGttt                                           <  1:91175/149‑1 (MQ=255)                                                                                                           gtCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATg                                     <  1:45203/151‑1 (MQ=255)                                                                                                             cGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGt                                   <  2:456497/151‑1 (MQ=255)                                                                                                             cGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGt                                   >  2:276035/1‑151 (MQ=255)                                                                                                                   cGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCaa                              >  1:599955/1‑150 (MQ=255)                                                                                                                           cAATGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCgg                     >  1:593548/1‑151 (MQ=255)                                                                                                                            aaTGATAGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGt                    >  1:219950/1‑151 (MQ=255)                                                                                                                                  aGGCGGCACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGt              <  1:484388/151‑1 (MQ=255)                                                                                                                                       gcACGTTTAACAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCt         <  1:193133/151‑1 (MQ=255)                                                                                                                                               aaCAGCTTCCCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc  <  2:379990/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               TGCCGTGTTAATCATCGTGTTTATCGAGTCAATGGGGATGTTCCTGGCACTGGGTGAAATAGTCGGTCGTAAACTCTCTTCGCACGATATTATTCGCGGGCTGCGTGTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTT‑CCC‑CCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCC  >  NC_000913/3032149‑3032439

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