Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,839,384	T→A	intergenic (‑120/‑140)	ascG ← / → ascF	asc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,839,384	0	T	A	98.0%	156.4 / NA	50	intergenic (‑120/‑140)	ascG/ascF	asc operon transcriptional repressor, prpBC operon repressor/cellobiose/arbutin/salicin‑specific PTS enzymes, IIB and IC components
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (28/21);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGAT  >  NC_000913/2839251‑2839507                                                                                                                                      |                                                                                                                            aGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        >  1:58984/1‑139 (MQ=255)        tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        >  1:237869/1‑132 (MQ=255)        tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        <  2:237869/132‑1 (MQ=255)        tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                 >  1:110444/1‑139 (MQ=255)        tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                 >  2:1352287/1‑139 (MQ=255)        tcgtcATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGc                                                                                                                 >  2:1653556/1‑139 (MQ=255)                 ttttCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGt                                                                                                        <  1:1352287/139‑1 (MQ=255)                           tttACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTAtt                                                                                              >  2:332642/1‑139 (MQ=255)                             tACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCa                                                                                            >  2:1594058/1‑139 (MQ=255)                                    aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActctc                                                                                                                     <  1:958738/107‑1 (MQ=255)                                    aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActctc                                                                                                                     >  2:958738/1‑107 (MQ=255)                                    aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        <  1:413024/104‑1 (MQ=255)                                    aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        >  2:413024/1‑104 (MQ=255)                                    aaTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGAtttt                                                                                     >  2:393911/1‑139 (MQ=255)                                     aTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAt                                                                                                             <  1:1681516/114‑1 (MQ=255)                                     aTCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTAt                                                                                                             >  2:1681516/1‑114 (MQ=255)                                           gCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAc                                                                              <  2:110444/139‑1 (MQ=255)                                            cAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAAcc                                                                             <  1:332642/139‑1 (MQ=255)                                              gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                      >  1:1486917/1‑96 (MQ=255)                                              gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                      <  2:1486917/96‑1 (MQ=255)                                              gtgtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                           <  1:1594058/139‑1 (MQ=255)                                                gtACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGtt                                                                         >  1:1627947/1‑139 (MQ=255)                                                    tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                                                                     >  2:108344/1‑139 (MQ=255)                                                    tGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCtt                                                                     >  1:825181/1‑139 (MQ=255)                                                         cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                           <  1:1420171/128‑1 (MQ=255)                                                         cGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCgg                                                                           >  2:1420171/1‑128 (MQ=255)                                                           tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        >  1:1613444/1‑81 (MQ=255)                                                           tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCAct                                                                                                                        <  2:1613444/81‑1 (MQ=255)                                                           tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGAc                                                                                                     <  1:1102944/100‑1 (MQ=255)                                                           tttGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGAc                                                                                                     >  2:1102944/1‑100 (MQ=255)                                                                     gTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAGCCGGTTTCTTAATTCCGTTTCAGCATc                                                    <  1:393911/139‑1 (MQ=255)                                                                        cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                          >  2:961235/1‑66 (MQ=255)                                                                        cTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCa                                                                                                                          <  1:961235/66‑1 (MQ=255)                                                                         ttGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                      >  1:495567/1‑69 (MQ=255)                                                                         ttGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCActct                                                                                                                      <  2:495567/69‑1 (MQ=255)                                                                                  cggcggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACTGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGt                                       <  1:1653556/139‑1 (MQ=255)                                                                                     cggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCAc                                    <  2:825181/139‑1 (MQ=255)                                                                                     cggGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCAc                                    >  1:86166/1‑139 (MQ=255)                                                                                       ggAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGt                                  <  2:58984/139‑1 (MQ=255)                                                                                                        aCCGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTAc                 >  1:502320/1‑139 (MQ=255)                                                                                                          cGGTTTCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACAt               >  2:1263159/1‑139 (MQ=255)                                                                                                             tttCACAAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTCGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTa            >  2:1193690/1‑139 (MQ=255)                                                                                                                 acaAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAGCCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTAcc          <  1:392748/137‑1 (MQ=255)                                                                                                                 acaAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTAcc          >  2:392748/1‑137 (MQ=255)                                                                                                                  caAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGt                                       >  2:1164173/1‑107 (MQ=255)                                                                                                                  caAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGt                                       <  1:1164173/107‑1 (MQ=255)                                                                                                                  caAATATAAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTa       <  2:86166/139‑1 (MQ=255)                                                                                                                       ataAAAAATGAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGAt  >  1:971298/1‑139 (MQ=255)                                                                                                                                gAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCAt                                                     <  1:1262589/79‑1 (MQ=255)                                                                                                                                gAACAAATCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCAt                                                     >  2:1262589/1‑79 (MQ=255)                                                                                                                                      |                                                                                                                            AGCATCGTCGTCATCATTTTCATCCCTTTACACGCAATCAACGCAGTGTACTGCACCGTTTGCCGATTGTCCTTGCACAATCGGCGGGAAAAATATTCAGGTGACCGGTTTCACAAATATAAAAAATGAACAATTCACTCTCTTGCTTATTTAGTGACAACTATTCATGATTTTGTGAAACCGGTTTCTTAATTCCGTTTCAGCATCGGCATTTTTCCGTCACGTCGACTGATAACAACTACATCTACCCTACTGAT  >  NC_000913/2839251‑2839507

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