Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,406	A→T	G262G (GGA→GGT)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,406	0	A	T	92.9%	75.8 / ‑6.9	28	G262G (GGA→GGT)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (10/16);  minor base G (2/0);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCT  >  NC_000913/78171‑78647                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  gTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATTGGGGCGGCAGCTGGGCTGGGAATaaaatcg                                                                                                                                                                                                                                    >  2:36460/1‑245 (MQ=255) gTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATa                                                                                                                                                                                                                                          <  1:36460/245‑1 (MQ=255)           gattgCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTctgg                                                                                                                                                                                                                          <  1:318562/251‑1 (MQ=255)                  tttcttcTTCCGTCTGTTGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCGGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTTCATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTggctggct                                                                                                                                                                                                                    <  2:223009/244‑1 (MQ=255)                                                  gCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCGACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTGTGTCAAACGTTATGGTAATCGGCGAat                                                                                                                                                                                     <  1:112407/249‑1 (MQ=255)                                                         gAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTTCTGATGGGGACGGCAGCTGGTCTGGGAATACCAGCAATGATTCTGGCTGGCTTATATGTTAAAAGTTATTGTAAgcggcggattaatggcct                                                                                                                                                                            >  2:284080/1‑239 (MQ=255)                                                               gCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAg                                                                                                                                                                       <  1:13597/250‑1 (MQ=255)                                                                                 ggctggcAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGGACTggt                                                                                                                                                    >  1:115882/1‑249 (MQ=255)                                                                                             aGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGgat                                                                                                                                         <  2:208412/250‑1 (MQ=255)                                                                                                   gtatggATGTTAGTTGTCGCCTCGACGTGAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTTCCAGACAAACTGGCGGGTGTCCTGATTGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATAATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTGACACCGGATTGAtttt                                                                                                                                  <  2:281833/247‑1 (MQ=255)                                                                                                      cGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGAt                                                                                                                                                                                                                               >  2:102644/1‑155 (MQ=255)                                                                                                      cGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGAt                                                                                                                                                                                                                               <  1:102644/155‑1 (MQ=255)                                                                                                         aTGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTAtttca                                                                                                                            <  1:89801/251‑3 (MQ=255)                                                                                                            ttATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGGAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGGGGACTGTTTTACACCGGATTGATTTTAttttat                                                                                                                           >  1:115506/1‑246 (MQ=255)                                                                                                             tATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAg                                                                                                                         <  2:165143/250‑1 (MQ=255)                                                                                                               tttGTCGCATCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGGCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGGGTACTGTTTTACAACGGGTTGTTTTTAttttctagccc                                                                                                                      >  2:79957/1‑243 (MQ=255)                                                                                                                  gTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGCTAGGATTGCCAGACAAACTGGCGGGCTTCCTGATGGGGGCTGCCGCTGGTCTGGTAATCCCTGCAACGATTCTGTGTGGGTACTATGTCTAATGGTATAGCAAGCGGGGCATGATGGTCATAGCAGTGGGGGCACGCGGACTGTTTTACATCggcgttattttattttctagtcgtct                                                                                                                   >  2:95875/1‑227 (MQ=255)                                                                                                                             gTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGTCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTAtttc                                                                                                                             <  2:237250/230‑2 (MQ=255)                                                                                                                             gTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGTCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTAtttc                                                                                                                             >  1:237250/1‑229 (MQ=255)                                                                                                                               tAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGGCGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGGCATAGCAGTGGCGGCAGGAGTAATGTTTTACACCGGGTTGGTTTTTTTTCATAGCCGTCTGGCGTTGATGACg                                                                                                      >  2:16191/1‑251 (MQ=255)                                                                                                                                           gCAACACCATGTACATTATTGATATGCCGGTGTGGATCAGTAGCGAGTGAGGATTGCCAGACAAACTGGCGGGTTTCCTGATTGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGTGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACtttt                                                                                           <  1:150703/250‑1 (MQ=255)                                                                                                                                               cacCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGGCGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGGGGCGGCAGGAGGAATGTTTTTAACCGGGTTGATTTTTTTTCATAACCGTCTGGGGttttagtcgcttccaattttttac                                                                                       >  2:117982/1‑228 (MQ=255)                                                                                                                                                       aCATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGGTGGTCTGGAAATACCAGCAATGGTTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTAAAAGAAGTGGCGGGAGGAGGACTGTTTTACACACGATTGATTTTTTGTCATAGCCGTCTGGCGTTGATTACGCAGCAACTTTTTAACGCTGGAttt                                                                              >  1:282091/1‑251 (MQ=255)                                                                                                                                                                            ggATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTAGCGGCATTGTTGCGGGTAt                                                          <  1:294332/250‑1 (MQ=255)                                                                                                                                                                                cAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTTTTTTATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGGATTgggg                                                     >  1:265982/1‑250 (MQ=255)                                                                                                                                                                                  gTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATg                                                   <  2:238949/251‑1 (MQ=255)                                                                                                                                                                                                          aaCTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCt                           <  1:2935/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggCAGCTGGTCTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTATTTCATAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCt  <  1:320353/251‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  GTCTGGTATTGATTGCATTTATGCTTCCGTCTGTGGCGCGGGTAGAACTGCCGTCGGAAAATGCTTTATCAATGCAAGGTGGCTGGCAGGATAGTAACGTACGGATGTTATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACTGGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAGTACTGTTTTACACCGGATTGATTTTCTTTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCT  >  NC_000913/78171‑78647

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