Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,307,299	(G)5→6	coding (501/930 nt)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,307,294	1	.	G	75.6%	101.0 / 30.7	41	coding (506/930 nt)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand):  ref base . (1/9);  new base G (25/6);  total (26/15)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

GTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑A‑GGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTC  >  NC_000913/4307155‑4307424                                                                                                                                                |                                                                                                                                    gTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑a                                                                                                                                      <  2:1367927/140‑1 (MQ=255)    ttCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑T‑‑A‑ggg                                                                                                                                  >  1:218963/1‑140 (MQ=255)       cGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑A‑gggggg                                                                                                                               >  2:468223/1‑139 (MQ=255)                   aaGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGt                                                                                                                   <  1:99966/140‑1 (MQ=255)                        cgcgCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGcc                                                                                                              >  1:41840/1‑140 (MQ=255)                          cgcAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTg                                                                                                            >  2:49116/1‑140 (MQ=255)                               tGGGTAATTTCGGGGTTGTTGTTCGTCCCAGCGTTTTAGCGCCATTCCAGAGCAATTGGTTCCCGGGCACCGAGGATTGCCACACGCGTGGTGTTGGGTTATATCTCC‑CG‑G‑GGGGGATAT‑GCCCCAGTCCGCCTGCcaca                                                                                                       <  2:840107/140‑1 (MQ=255)                                         cGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCAc                                                                                             >  2:1455450/1‑140 (MQ=255)                                           gggCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCAccc                                                                                           >  1:996772/1‑140 (MQ=255)                                               tgttgtTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑CG‑G‑GGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTcc                                                                                       <  1:388838/140‑1 (MQ=255)                                                        tACCAGCTTCTTAGCGCCATTCCAGAGCAATTGTTTCCCAGGCACCAAGGATTGCCACCCGCGCAGTGTGGGGTCATATCCCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAc                                                                              <  2:451813/140‑1 (MQ=255)                                                        tACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCAc                                                                              <  1:49116/140‑1 (MQ=255)                                                         cccAGCGTCTTAGCCCCATCCCAGAGAATTTGGTTCCCAGCCCCCCAGGATTGCCACCCGCGCAGTGTGGGGTCATATCCCC‑C‑‑A‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                            <  2:522699/139‑1 (MQ=255)                                                         aCCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCCCC‑C‑‑AGGGGGGATAT‑GACCCAGTTC‑CCTGCCACACCGTGTGCACCCGTCCACGGCGCACCg                                                                            <  2:1058584/140‑1 (MQ=255)                                                              cGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                        >  1:1021737/1‑140 (MQ=255)                                                              cGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                        >  1:369260/1‑140 (MQ=255)                                                              cGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCa                                                                        >  1:1395489/1‑140 (MQ=255)                                                              cGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑A‑GGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCAt                                                                       >  1:522699/1‑140 (MQ=255)                                                                tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATGTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCAACCGTGCACGGCGCACCGTTGATc                                                                      >  2:1029566/1‑140 (MQ=255)                                                                tCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTG‑‑CCCGTCCACGGCGCACCGTTCATCca                                                                    >  2:338627/1‑140 (MQ=255)                                                                 cTTAGCGCCATTCCAGAGAAATTGTTTCCCAGGCACCAAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcc                                                                     <  2:569508/140‑1 (MQ=255)                                                                  ttAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATcct                                                                    >  1:1254616/1‑139 (MQ=255)                                                                        ccATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGc                                                              >  2:362998/1‑140 (MQ=255)                                                                          atcccAGAGCAATTGGTTCCCGGCCCCCGAGGATGGCCCCACGCGCAGTGTGGGGT‑ATTTCCCC‑CGGG‑GGGGGATAT‑GACCCGGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGa                                                            <  2:362889/137‑1 (MQ=255)                                                                           ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGaa                                                           >  1:582264/1‑140 (MQ=255)                                                                           ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGaa                                                           >  1:744412/1‑140 (MQ=255)                                                                           ttCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGaa                                                           >  1:362889/1‑140 (MQ=255)                                                                              cggagCATTTGGTCTCCAGGCCCCCAGGATTGCCACACGCGCGGTGTGGGGTAATACCCCCGG‑‑G‑GGGGGATAG‑CACCCAGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                        <  2:1068277/138‑1 (MQ=255)                                                                              cAGAGCATTTGGTTCCCAGGCACCCAGGATTGCCACACGCGCAGTGTGGGGTAATATCCCCGG‑‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                        <  2:1395489/140‑1 (MQ=255)                                                                              cAGAGCAATTGTTTCCCAGGCCCCCAGGATTGCCACCCGCGCAGTGTTGGGTAATATCCCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                        <  2:1021737/140‑1 (MQ=255)                                                                              aagagCATTTGTTTCCCAGCCCCCCAGGATGCCCCCACGCGCGGTGTGGGGTCATCCCCCC‑CG‑G‑GGGGGATAT‑GACCCGGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAAccc                                                        <  2:1150127/139‑1 (MQ=255)                                                                                gagCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCa                                                      >  1:402245/1‑140 (MQ=255)                                                                                     aTTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATcccc                                                 >  2:1107174/1‑140 (MQ=255)                                                                                      ttGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCg                                                >  2:1177341/1‑140 (MQ=255)                                                                                        ggTTCCCAGGCCCCCAGGATTCCCACCCGCGCAGTGTGGGGTCATATCTCC‑CG‑G‑GGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTa                                              <  1:468223/140‑1 (MQ=255)                                                                                           ttCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCg                                           >  2:950740/1‑140 (MQ=255)                                                                                                  cccccAGGATTGCCACACGCGCGGTGTTGGGTCATATCCCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAgg                                    <  2:109640/138‑1 (MQ=255)                                                                                                   aCCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTACATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgc                                   >  1:785993/1‑140 (MQ=255)                                                                                                       aGGATTGCCCCCCGCGCGGTGTTGGGTCATACCCCC‑CG‑G‑GGGGGATAT‑GACCCGGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgcc                               <  2:402245/140‑1 (MQ=255)                                                                                                       aGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGgccgcc                               >  1:109640/1‑140 (MQ=255)                                                                                                        ggATTCCCACCCGCGCAGTGTTGGGTTATCTCCCCGG‑‑G‑GGGGGTTAT‑GCCCATGTCCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCa                              <  2:132617/140‑1 (MQ=255)                                                                                                            tGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAAc                          >  2:1336434/1‑140 (MQ=255)                                                                                                               cacaCGCGCAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAg                       >  2:819985/1‑140 (MQ=255)                                                                                                                 cccgcgcgGTGTTGGGTCATTTCCCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAgtt                     <  2:509760/138‑1 (MQ=255)                                                                                                                      gcAGTGTTGGGTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGc                <  1:1177341/140‑1 (MQ=255)                                                                                                                            ttGGGTCATATCCCC‑CG‑G‑GGGGGATAT‑GCCCCTTTCCCCCTGCCACACCGTGTG‑‑CCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg        <  1:338627/140‑1 (MQ=255)                                                                                                                            tggggTCATATCTCCGG‑‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAgg          <  1:777522/138‑1 (MQ=255)                                                                                                                              gggTCATCCCTCC‑A‑‑A‑GGGGGAAATGGACCAGGTTGCCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg        <  1:1029566/140‑1 (MQ=255)                                                                                                                              gggTCATACCCCC‑CG‑G‑GGGGGATAT‑GCCCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCg        <  2:820780/140‑1 (MQ=255)                                                                                                                               ggTCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCgg       <  1:950740/140‑1 (MQ=255)                                                                                                                                 tCATATCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGtt     <  1:1336434/140‑1 (MQ=255)                                                                                                                                    tatCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc  >  1:820780/1‑140 (MQ=255)                                                                                                                                    tatCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc  >  1:1156501/1‑140 (MQ=255)                                                                                                                                    tatCTCC‑C‑‑AGGGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTc  >  1:910374/1‑140 (MQ=255)                                                                                                                                                |                                                                                                                                    GTTTTCCGCATGGACGAAAAGATCGCGCAATGGGTAATTTCGGGGCTGTTGTTCGTACCAGCGTCTTAGCGCCATTCCAGAGCAATTGGTTTCCAGGCACCCAGGATTGCCACACGCGCAGTGTTGGGTCATATCTCC‑C‑‑A‑GGGGGATAT‑GACCCAGTTCGCCTGCCACACCGTGTGCACCCGTCCACGGCGCACCGTTCATCCACACTGCGAACCCCATCCCCGTACCGAGATAGGCCGCCAGAACCAGTTGTTGCGTAAGGCGGTTTTC  >  NC_000913/4307155‑4307424

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