Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,299,273	Δ1 bp	10.7%	intergenic (+28/‑1909)	ychE → / → oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,299,273	0	T	.	10.7%	176.6 / 15.1	56	intergenic (+28/‑1909)	ychE/oppA	UPF0056 family inner membrane protein/oligopeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base . (3/3);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299126‑1299413                                                                                                                                                     |                                                                                                                                             gTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATg                                                                                                                                             >  2:537890/1‑149 (MQ=255)     tACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTcc                                                                                                                                         <  4:140264/149‑1 (MQ=255)        gCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCA‑GCTCCGtt                                                                                                                                      <  4:202013/149‑1 (MQ=255)        gCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGtt                                                                                                                                      <  2:500472/149‑1 (MQ=255)           aGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTAc                                                                                                                                   <  4:225960/149‑1 (MQ=255)           aGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTAc                                                                                                                                   <  4:261042/149‑1 (MQ=255)                ggCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGt                                                                                                                              <  3:248397/149‑1 (MQ=255)                      aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCA‑GCTCCGtttcatttga                                                                                                                              >  3:408585/1‑136 (MQ=255)                      aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGCATGGCA‑GCTCCGtttcatttga                                                                                                                              <  4:408585/143‑8 (MQ=255)                      aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattat                                                                                                                        >  2:401304/1‑149 (MQ=255)                      aaCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattat                                                                                                                        >  1:6123/1‑149 (MQ=255)                       aCGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattatt                                                                                                                       >  2:481421/1‑149 (MQ=255)                       aCGTGATGACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCattatt                                                                                                                       >  2:399133/1‑149 (MQ=255)                           gATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGACATGCTCCGTTTACTTCGTCATTATTTTTa                                                                                                                   <  1:440355/149‑1 (MQ=255)                             ttACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACt                                                                                                                 <  1:481421/149‑1 (MQ=255)                              tACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACtt                                                                                                                <  4:245490/149‑1 (MQ=255)                                 gcgTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGttt                                                                                                                                     <  1:164358/125‑1 (MQ=255)                                 gcgTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGttt                                                                                                                                     >  2:164358/1‑125 (MQ=255)                                     attatGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTccc                                                                                                         <  3:276388/149‑1 (MQ=255)                                                              tGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGaa                                                                                <  4:265834/149‑1 (MQ=255)                                                              tGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGaa                                                                                <  3:415446/149‑1 (MQ=255)                                                               ggggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAAt                                                                               >  1:295293/1‑149 (MQ=255)                                                                 ggATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAg                                                                             <  3:9192/149‑1 (MQ=255)                                                                       aaTTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTaaaa                                                                       <  1:537890/149‑1 (MQ=255)                                                                          ttATTGTTACTGGTATTAGGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGATCCGTTTAAATCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAACGGAATAGGTAAAAata                                                                    <  3:410011/149‑1 (MQ=255)                                                                           tATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAatat                                                                   <  4:23742/149‑1 (MQ=255)                                                                                ttACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttcttc                                                              <  2:305084/149‑1 (MQ=255)                                                                                         ttAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  >  1:200475/1‑136 (MQ=255)                                                                                         ttAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAtt                                                                  <  2:200475/136‑1 (MQ=255)                                                                                             ggggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagtt                                                 >  2:578452/1‑149 (MQ=255)                                                                                              gggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagtta                                                >  3:442058/1‑149 (MQ=255)                                                                                               ggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttag                                               <  4:3792/149‑1 (MQ=255)                                                                                               ggATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttag                                               >  2:586310/1‑149 (MQ=255)                                                                                                gATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagt                                              >  1:163473/1‑149 (MQ=255)                                                                                                  ttttCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACa                                                    >  3:130120/1‑141 (MQ=255)                                                                                                  ttttCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATACTCTTCTCAAATTAcc                                                    <  4:130120/141‑2 (MQ=255)                                                                                                     tCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATaa                                         <  3:311283/149‑1 (MQ=255)                                                                                                             cTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTcc                                 <  3:407722/149‑1 (MQ=255)                                                                                                                 ttaattaaTTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTaa                             <  1:190085/149‑1 (MQ=255)                                                                                                                  taattaatTCCTCCCAAATGATACGGATC‑TGCCACGCTCCGGTTACTTCATCATCAATATTACATTGTTCCCGCGTAGATAGCAAAAGCAAAAGGAATAGGTAAAAATATTGTTCCCAAATCACAGTTAGTTATAAGGATTTCCTCAAc                            <  4:331697/149‑1 (MQ=255)                                                                                                                    attaattCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagtt                                             >  2:426897/1‑130 (MQ=255)                                                                                                                    attaattCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACagttagtt                                             <  1:426897/130‑1 (MQ=255)                                                                                                                        attCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtc                      <  3:10594/149‑1 (MQ=255)                                                                                                                        attCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTtc                      <  1:399133/149‑1 (MQ=255)                                                                                                                            cTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcct                  >  2:298952/1‑149 (MQ=255)                                                                                                                             tttCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTctcctc                 >  1:40747/1‑149 (MQ=255)                                                                                                                                   agtaaaCGGAGCATGGCA‑GCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatc           >  3:212504/4‑149 (MQ=255)                                                                                                                                    tGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatca          >  4:344056/1‑149 (MQ=255)                                                                                                                                     gAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAATTTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACcatcat         <  3:564993/149‑1 (MQ=255)                                                                                                                                      aaaCGGAGCATGGCA‑GCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 >  2:231321/1‑92 (MQ=38)                                                                                                                                      aaaCGGAGCATGGCA‑GCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATAttc                                                                 <  1:231321/92‑1 (MQ=38)                                                                                                                                          ggAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAt    >  2:23491/1‑149 (MQ=255)                                                                                                                                          ggAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAt    <  3:446198/149‑1 (MQ=255)                                                                                                                                            aGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTAttt  <  3:344056/149‑1 (MQ=255)                                                                                                                                                   cATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagcaaaag                                                                                   >  3:499501/1‑62 (MQ=255)                                                                                                                                                   cATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATcaaaagcaaaag                                                                                   <  4:499501/62‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GTTTTACGCCAGACCGGCATCAACGTGATTACGCGTATTATGGGGCTATTGCTGATGGCATTGGGGATTGAATTTATCGTTACTGGTATTAAGGGGATTTTCCCCGGCCTGCTTAATTAATTCCTTTCAAATGAAACGGAGC‑TGCCATGCTCCGTTTACTTCGTCATTATTTTTACTTTGTTCCCGCGCAGTTATCAAAAGCAAAAGGAATAGGTAAAAATATTCTTCTCAAATTACAGTTAGTTATAAGGATTTCCTTAACTGCTTCTCCTCACCATCATGTTATTT  >  NC_000913/1299126‑1299413

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