Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	78,761	G→A	V381I (GTC→ATC)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	78,761	0	G	A	93.3%	74.3 / ‑0.8	30	V381I (GTC→ATC)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (13/15);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

TTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCT  >  NC_000913/78522‑79001                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ttcaTAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGAATTGTTGCGGGTATTGAGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCAACCACCTTATTTACTAACAGCATTTCTACCGGGGTAATTCTGGCTGGCGTTCTTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTAttt                                                                                                                                                                                                                                       >  1:282906/4‑251 (MQ=255)   caTAGCCGTCTGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTAttttt                                                                                                                                                                                                                                     <  1:142829/250‑1 (MQ=255)               gCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGt                                                                                                                                                                                                                         <  1:356453/251‑1 (MQ=255)                    gatgaCGCTGCAACTTGTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAg                                                                                                                                                                                                                    <  1:58815/251‑1 (MQ=255)                                   ttttAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTAAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGActtg                                                                                                                                                                                                     >  2:26875/1‑248 (MQ=255)                                                  tATCGGCATTGGTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGaaa                                                                                                                                                                                       <  2:157678/250‑1 (MQ=255)                                                   aTCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCATGGAAGAGCGGGGGCAGCTACCCCCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTAATGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGTTAGCGGAAAGc                                                                                                                                                                                     >  1:132648/1‑251 (MQ=255)                                                                      gtGGGATGCTATGGTTTCAGGAATTAATGCCTGGTAGAGCGGGGGCAGCTACCGCCTTAGTGACTAACAGTAGTTCTACCGGTGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCa                                                                                                                                                                  <  2:221536/250‑1 (MQ=255)                                                                                  ggTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCTGTTATTTCTTTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTTATGACGTTGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAAtt                                                                                                                                                      <  1:236994/251‑1 (MQ=255)                                                                                   gTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTc                                                                                                                                                     >  1:356424/1‑251 (MQ=255)                                                                                        aGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATaaa                                                                                                                                                <  1:221113/251‑1 (MQ=255)                                                                                         ggATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAAc                                                                                                                                               >  1:184730/1‑251 (MQ=255)                                                                                                aTGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCG‑CATTATTTTTAACCGCAAAGGTTAAATACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTATTTCATAAACGCAGgtt                                                                                                                                        >  1:165238/1‑251 (MQ=255)                                                                                                              gggggCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGAAGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATa                                                                                                                           >  1:314133/1‑250 (MQ=255)                                                                                                               agtacagcgcccaccTTATTTACTAGCAGTATTTCTACCGGGGGAATACGGTCTGGCGTTATTAAGGGAGCAATTGCACAAAGGTCGTGGCACGTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCGACAT‑ACTTGTAACCGCACAGTCTAAAGACGTTTGATGACGTAGACCATAACGGACAGCCCGGTCATTTTTCCCGGCAAGGGGCGTGATGAGTAAACGCAGGTTGTTTTGCTTccta                                                                                                                           >  2:106250/11‑246 (MQ=255)                                                                                                                        accaccTTATTTACTAACAGTATTTATACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAAt                                                                                                                  <  2:317083/249‑1 (MQ=255)                                                                                                                              atATTTACTAACAGTATTTATACCGGGGTAATTCTGGCTGTCGTTATTCATGGAGCATTTGAACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTCTTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGGCGTTTGATGACGTTGACGATAGCGGCAAGCCCGGTCATTTGACGAGGCCAGGGGATTAATTCATAAACGCAGGTTGTTTTGATTCATAAGCGGCACTGGcggcg                                                                                                           >  2:162582/2‑247 (MQ=255)                                                                                                                                           gTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                             <  1:256020/251‑1 (MQ=255)                                                                                                                                            tATTTCTACCGGGGTAATTATGGCTGGCGTTATTCAGGGAGCAATTGCAAAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGAAGTGGACGATAGCGGAAAGCCCGGTCATTTGAACGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCCTAAGCGGCAATGGCGTCGTCGTGCTGCAAGGt                                                                                             >  2:65162/1‑250 (MQ=255)                                                                                                                                                  tACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGTAGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGccc                                                                                       >  1:321740/1‑250 (MQ=255)                                                                                                                                                                 ggctggcGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGAc                                                                         <  1:276298/249‑1 (MQ=255)                                                                                                                                                                           aTTCAGGGAGCAATTGCACAAAGTTGGGGGCCATTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCACAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtcg                                                                                                        >  1:408735/1‑208 (MQ=255)                                                                                                                                                                           aTTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATGCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGcgtcgtcg                                                                                                        <  2:408735/208‑1 (MQ=255)                                                                                                                                                                                  gAGCAATTGCACAAAGTTGGGGGCCCTGTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTAGCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTggcg                                                      <  2:321740/251‑1 (MQ=255)                                                                                                                                                                                              aaaGTTGGGGGCACTTTGCTGTCTACTGGGTAATGGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATc                                          <  2:96658/251‑1 (MQ=255)                                                                                                                                                                                                          aCTTTTCTGTCTCCTGGGTAATTGCTGGTATTGTTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTTGACGATAGCGGAAAGCCCGGTCATTTTACCTGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCTGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATTCAGTGTCTGCGGAAGGCACCGATGGTAAAGCg                              <  1:374964/251‑1 (MQ=255)                                                                                                                                                                                                                    cTACTGGGTAATTGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTtc                    >  1:224781/1‑251 (MQ=255)                                                                                                                                                                                                                           gTAATTGCGGTTATTTCTGTTATCG‑CCTTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCg             <  2:459080/251‑1 (MQ=255)                                                                                                                                                                                                                                tGCGGTTATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCttt         <  1:254473/250‑1 (MQ=255)                                                                                                                                                                                                                                      tATTTCTGTTATCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACctct  >  1:124364/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  TTAATAGCCGTATGGCGTTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTCTGTTGTCG‑CATTATTTTTAACCGCAAAGGTTAAAGACGTTTGATGACGTGGACGATAGCGGAAAGCCCGGTCATTTGACCGGGCAAGGGGATTAATTCATAAACGCAGGTTGTTTTGCTTCATAAGCGGCAATGGCGTCGTCGTGCTGCAAGGTAAGCCCAATACTGTCCAGACCGTTCATCATGCAGTGGCGGCGGAAGGCATCGATGGTAAAGCGATAGGTTTTCTCTCCCGCTTTCACCTCT  >  NC_000913/78522‑79001

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