Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,309,712	C→T	W423* (TGG→TGA)	alsA ←	D‑allose ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,309,712	0	C	T	97.0%	104.4 / ‑6.2	33	W423* (TGG→TGA)	alsA	D‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/17);  minor base A (0/1);  total (15/18)
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.54e-01

CATGCCATAATCTCCTCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGA  >  NC_000913/4309467‑4309953                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  cATGCCATAAGCGCCGCTTCTCTCATGTCATCGCGATTCTTCAGGAGTTGCGTCTGTCTTCCTTCGCTGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGTTAGTTCAGATTACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCAGCGATTTGTAAATTTCTGCTTTCGCGCCAACGTCTATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCAttt                                                                                                                                                                                                                                              <  2:273210/251‑1 (MQ=255)           tctcCTCTTCGCTCATGTAATCGCGATTCGTCTGGATTTGCGTCAGTCGTCCTTCGCAGCACACGGCTATGCGGTCGCAGACGGTGATAATTTCAGGTCGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGACAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCCACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCCTTTGGAGATCAgg                                                                                                                                                                                                                                    >  2:224366/1‑251 (MQ=255)              cctcTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGAc                                                                                                                                                                                                                                  <  1:326940/250‑1 (MQ=255)                       tCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTctgctg                                                                                                                                                                                                                         >  2:70741/1‑250 (MQ=255)                        cATGTCAACGCGATTCGTCAGGATTTGCGTCAGTCGTCATTCGCAGAACACGGCGAGGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTGTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTAGCGCCAACGTCGATGCCGCGGGTAGGTGCATCGAAAATAATCACTGCCGGGCAACAGCACAGACATTTGGAGATCAGGACTTTctgctg                                                                                                                                                                                                                         <  2:129214/249‑1 (MQ=255)                                cgcgATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTGTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGACCAGGACTTTCTGCTGATTTCCCCCg                                                                                                                                                                                                               <  2:173203/251‑1 (MQ=255)                                        gTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTc                                                                                                                                                                                                       >  2:140518/1‑251 (MQ=255)                                        gTCAGGATTTGCGTCAGTCGTCCTTCGAAGAACACGGCGATGCGGTAGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTAGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGAGGGTAGGTTCATCGAAAATAATCAATTCCAGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGCTTTCCCCCGGAGAGTTc                                                                                                                                                                                                       >  2:223144/1‑251 (MQ=255)                                          cAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCgg                                                                                                                                                                                                     <  1:272216/251‑1 (MQ=255)                                            ggATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTg                                                                                                                                                                                                   <  1:316400/251‑1 (MQ=255)                                                   cGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATTACTTTTCCGGCGTCCGCCAGTTGGCGCAGCACTTTGTAAATTTCGTCTTGCGCGCCAACGTCGATGCCGCGGGTATGTTCAGCGAAAATAATCACTTCCGGGCAACAGCGCATTCATGTGGAGTTCAGTACTGTCTACTGATTTCCCCCGGATAGTTCGTTGATATTCt                                                                                                                                                                                            <  2:209333/251‑1 (MQ=255)                                                      cAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTgg                                                                                                                                                                                          <  1:70741/250‑1 (MQ=255)                                                          cgtcCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTAc                                                                                                                                                                                     >  2:56186/1‑251 (MQ=255)                                                          cgtcCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCAGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTAc                                                                                                                                                                                     >  2:114426/1‑251 (MQ=255)                                                                                  cGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCcagcag                                                                                                                                                             <  1:223642/251‑1 (MQ=255)                                                                                        cAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATTACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCTCTTCCGGGCAACAGCACAGTCATTTGGAGATCATGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTcgcg                                                                                                                                                       <  1:62438/251‑1 (MQ=255)                                                                                               tGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGAttt                                                                                                                                                >  1:279949/1‑251 (MQ=255)                                                                                                                   gATGACCCCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTAATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTc                                                                                                                            >  2:101631/1‑251 (MQ=255)                                                                                                                               aGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTTATTTTCAGCGGTACGTTGCTCGTCAACTTCATGaaa                                                                                                                <  1:156546/251‑1 (MQ=255)                                                                                                                                    gACTTTTCCGTCGTCCGCCAGGTGGCGGATCACTTTGTAAAGTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGTAGAGTTCGGTGTTATGCTGGTTTACCGAATGTCATTTCAGCGTCAGCAGTTCGTGTTGATTTTGATCGGTACGTTGCTCGTCAACTTCATGAAACAAGc                                                                                                           <  1:200013/251‑1 (MQ=255)                                                                                                                                     ccTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGc                                                                                                                                 <  1:87712/227‑1 (MQ=255)                                                                                                                                     aCTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGc                                                                                                                                 >  2:87712/1‑228 (MQ=255)                                                                                                                                        tttCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGTAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCa                                                                                                        <  2:311008/250‑1 (MQ=255)                                                                                                                                            cgtcgtCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGCTTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATcgcg                                                                                                   >  1:276400/1‑251 (MQ=255)                                                                                                                                                           ggCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                     <  1:34301/250‑1 (MQ=255)                                                                                                                                                                                 ggatttCGCGCCAACGGCGATGCCGCGGGTAGGTTCATCGAAAATACTCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAAGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGAAATTTCAGCGCCAGCAGTTAGCGTTGATTTTCAGCGGTACTATTCTCGTCAACTTCATGAAACAAGCACATCGCGCCTTTATAACAGCCGTCGTTCAGACTGCGGCTGATCg                                                              >  2:995/4‑251 (MQ=255)                                                                                                                                                                                  gCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGc                                                             >  1:165901/1‑251 (MQ=255)                                                                                                                                                                                                             gTAGGTTCATCGAAACTAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACAtt                                                                                                                                                                          >  2:19226/1‑115 (MQ=255)                                                                                                                                                                                                             gTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACAtt                                                                                                                                                                          <  1:19226/115‑1 (MQ=255)                                                                                                                                                                                                                   tCATCGAAAATAATCACTTCCGGGCAACAGCACAGTCAGTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaaa                             <  1:28248/250‑1 (MQ=255)                                                                                                                                                                                                                       cGAAAATAATCACTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGCCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATCTTCTGGTTTACCGAATGAAATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                        >  2:324301/1‑251 (MQ=255)                                                                                                                                                                                                                                   cTTCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCa                                                                                                                                             <  2:94074/122‑1 (MQ=255)                                                                                                                                                                                                                                    ttCCGGGCAACAGCACAGTCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCggcggc           >  1:230559/1‑251 (MQ=255)                                                                                                                                                                                                                                             acagcacagtCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGa  <  2:279949/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  CATGCCATAATCTCCTCTTCGCTCATGTCATCGCGATTCGTCAGGATTTGCGTCAGTCGTCCTTCGCAGAACACGGCGATGCGGTCGCAGACGGTGATAATTTCAGGTAGTTCAGATGACACCATCAGGATGACTTTTCCGTCGTCCGCCAGTTGGCGCATCACTTTGTAAATTTCGGCTTTCGCGCCAACGTCGATGCCGCGGGTAGGTTCATCGAAAATAATCACTTCCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGA  >  NC_000913/4309467‑4309953

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