Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,139,511	Δ1 bp	coding (198/663 nt)	fsaB ←	fructose‑6‑phosphate aldolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,139,511	0	C	.	97.0%	147.0 / ‑2.7	33	coding (198/663 nt)	fsaB	fructose‑6‑phosphate aldolase 2
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (13/19);  total (13/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.70e-01

GTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC  >  NC_000913/4139275‑4139748                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              gTACGAATGCCGTTTCCGCCCTGGGCATCGACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGTGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTc                                                                                                                                                                                                                                <  2:208416/251‑1 (MQ=255)  tACGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTc                                                                                                                                                                                                                                >  1:171400/1‑250 (MQ=255)   aCGAATGCCGTTTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGCCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCa                                                                                                                                                                                                                               <  1:117133/250‑1 (MQ=255)            gtttcCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGc                                                                                                                                                                                                                     <  2:113796/248‑1 (MQ=255)                    ccTGGGCATCTACTCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGaa                                                                                                                                                                                                             <  2:357433/251‑1 (MQ=255)                    ccTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGaa                                                                                                                                                                                                             <  1:327253/251‑1 (MQ=255)                      tGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATg                                                                                                                                                                                                           >  1:230320/1‑251 (MQ=255)                                         cATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTGTAATTGCTTCCAGACCTTCGGAAGGCGACGGGATCTTCACCACAATTCCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGATTGCCCGCATAACCAATCGCtt                                                                                                                                                                                        <  2:134138/251‑1 (MQ=255)                                            aCGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCttttt                                                                                                                                                                                     >  1:310950/1‑251 (MQ=255)                                              ggggCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGc                                                                                                                                                                                   <  2:199366/248‑1 (MQ=255)                                                    aCGTATTTTGCCCCTGCCAGTGCGGATAATAAACCTTGTGCGGCGCTAGATACAGCAGTGCCAAGTGTAGTAATACCCTCTGTTATCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGTCGCTTCGCTTCTTCCACCA‑ACCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCAGCACCAATCGCTTTTTGCAGAcgc                                                                                                                                                                             <  2:263669/251‑1 (MQ=255)                                                                ccTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTcc                                                                                                                                                                 >  1:226536/1‑251 (MQ=255)                                                                   gCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCat                                                                                                                                                              <  1:384241/251‑1 (MQ=255)                                                                                        tgtgCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGc                                                                                                                                          <  1:366865/250‑1 (MQ=255)                                                                                             ggcgCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATTACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAAt                                                                                                                                    >  2:19948/1‑251 (MQ=255)                                                                                                                 cAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTATTTGTCAccc                                                                                                                <  1:158148/251‑1 (MQ=255)                                                                                                                      gtagtaATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGc                                                                                                            >  1:74300/1‑250 (MQ=255)                                                                                                                                ccTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAat                                                                                                  >  1:311206/1‑250 (MQ=255)                                                                                                                                            gTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGAc                                                                                     <  1:88357/251‑1 (MQ=255)                                                                                                                                               ttttAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGtt                                                                                  >  1:195006/1‑251 (MQ=255)                                                                                                                                                                               gggATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCGTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                                                  <  1:289316/251‑1 (MQ=255)                                                                                                                                                                                  aTTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATTACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                                               <  1:19948/251‑1 (MQ=255)                                                                                                                                                                                          cacaATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgat                                       <  2:310950/251‑1 (MQ=255)                                                                                                                                                                                            caATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt                                     <  1:89538/251‑1 (MQ=255)                                                                                                                                                                                              aTACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGttt                                   >  2:256421/1‑251 (MQ=255)                                                                                                                                                                                                ccccGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCACCCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTGTGCAGACGCGGCAGCACTTCCCATATGGACTCCGTGCTGGCAGCGATAATGCTCGGGTTAGTTGTGTCCCCGGCAATGGGGAATATGCGTGCCAGACGTGCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtttac                                 <  2:358652/250‑3 (MQ=255)                                                                                                                                                                                                 cccGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt                                <  2:351772/251‑1 (MQ=255)                                                                                                                                                                                                    ggAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCCCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt                             <  1:160441/251‑1 (MQ=255)                                                                                                                                                                                                     gAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCtttt                             >  1:342808/1‑250 (MQ=255)                                                                                                                                                                                                           cGTCGCGCATGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTGGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttc                      <  1:286968/251‑1 (MQ=255)                                                                                                                                                                                                            gTCGCGCAGGCGCTTCGCTTCTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttcc                     >  1:176214/1‑251 (MQ=255)                                                                                                                                                                                                                           cgcttcTTCCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTgcagc      >  2:101035/1‑251 (MQ=255)                                                                                                                                                                                                                               tcttcCACCATCCC‑TGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGc  >  1:109922/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                              GTACGAATGCCGTCTCCGCCCTGGGCATCTACGCGGTTAACATACGGAGCAACGTATTTTGCCCCTGCCAGTGCGGCTAATAACCCTTGTGCGGCGCTATATACAGCAGTGCCAAGTGTAGTAATACCCTCTTTTTTCAGTATTTTAATTGCTGCCAGACCTTCGGAAGTCACCGGGATTTTCACCACAATACCCGGAATAGCGTCGCGCAGGCGCTTCGCTTCTTCCACCATCCCCTGCGCGTCGCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC  >  NC_000913/4139275‑4139748

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