Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,804	T→G	A357A (GCT→GCG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,804	0	T	G	93.5%	70.8 / ‑4.3	32	A357A (GCT→GCG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (11/18);  total (13/19)
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.92e-01

GAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCAT  >  NC_000913/4558‑5031                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      aaccgacaCCGTCCCACGTTCCCTGAACGAGGGTCAGGGGTCCCCCAAAGCGCATAAGGAGACGTTATCCAACGCGAGGGAGGGTAGCCAGCCGCACAGCGGTCCGCGATG‑GGAAGATTTGGTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTGAACCTCGGAGCCGCGCGCTGCCGGAGCGTATc                                                                                                                                                                                                                                   <  2:188350/244‑1 (MQ=255)    cGATCCCGTTCCACGTTTCCTGCCCGACGGTCAGTGGTCACCCAAAGCGACGCAGGCGACGTTATCAAACGCGATGGCCTTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTg                                                                                                                                                                                                                                <  1:47416/251‑1 (MQ=255)      aTACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTgc‑g                                                                                                                                                                                                                             >  1:152739/1‑251 (MQ=255)              ccaagTTTCCGGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGGCGTTATCCAACGAGATGGTCGTGAGTCAGTCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGTTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCGACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GATGCGTg                                                                                                                                                                                                                      <  1:343431/246‑1 (MQ=255)                             gacgGTCAGTGATCACCCAAAACGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTACGCCGCAAAATCTGGAAACTGAAAGAGCTGGGATATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGCACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCGGCGTGACCAGGTGAATCCAgg                                                                                                                                                                                                      >  2:370674/1‑251 (MQ=255)                              acagcgcgagggaccccAAAACGCCTCGGGCGACTTTATACACCGCTAAGGACGTGAGTCAGCTAAACAACTGGCCGCG‑TGGGAAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCTGCACGCTGCCGTAGCTTATCGTGT‑GCTGCGTGACCAGTTGAATCCAGGc                                                                                                                                                                                                     <  2:200407/237‑1 (MQ=255)                                          cccccAAAGCGACTCAGGCCACGTTATCCAACGCGTTGGCCGTGATTCGGCCGAACACCTGGCCGGT‑TGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGGGAATATGGCTTg                                                                                                                                                                                         <  2:222125/249‑1 (MQ=255)                                                 aGCGACTCAGGCGGCGTTATCCAACGCGAGGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGGGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                  <  1:101714/251‑1 (MQ=255)                                                              aCGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATcc                                                                                                                                                                    <  1:379567/251‑1 (MQ=255)                                                               cGTTATCCAACGCGATGGACGTGAGTCAGCCGAACCACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGAGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGAGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTGCGCACC‑GCGCCTCCg                                                                                                                                                                   >  2:186860/1‑251 (MQ=255)                                                                         cgcgATGGCCGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAAttt                                                                                                                                                          >  1:9081/1‑250 (MQ=255)                                                                                    tGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGGAATTTAAagagggagtg                                                                                                                                              >  1:197684/1‑245 (MQ=255)                                                                                               aacaacGGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtct                                                                                                                                   <  1:141914/251‑1 (MQ=255)                                                                                                 caacTGGCAGCG‑TGTGGAAGCGTTGTTCCGCCGCAAAAACTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATACAGGCGCATATGGAGTGTTCCTACGCACCTCCGC‑CCCGGCGAAATTTTAAGAGAGCGGGTAAGCGGTtcgcg                                                                                                                                 >  2:35345/1‑248 (MQ=255)                                                                                                         ccggtgggggAAGAGTGATTCCGCCGCAAAATCTGGCAACGGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACGCCTCGGAGCCGCACGCTGCCGTAGCGTAGCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAc                                                                                                                          <  2:118377/243‑1 (MQ=255)                                                                                                                tggcagagTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGTCGCACGCGGCCGAAGAGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCATCGGCACC‑GCGCATCCGGCGAAATTTAAAAAGGGGGTGGTATCGATTCCCGGTGGAACGtaggtt                                                                                                                   >  2:85658/5‑246 (MQ=255)                                                                                                                gggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                   <  1:216401/250‑1 (MQ=255)                                                                                                                 ggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATc                                                                                                                  >  1:262021/1‑251 (MQ=255)                                                                                                                                      aaTCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTTAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATTCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCCAAAAGAGCTGGCAGAAc                                                                                             <  1:184688/251‑1 (MQ=255)                                                                                                                                          tGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGc                                                                                         >  1:166399/1‑251 (MQ=255)                                                                                                                                                  gAAAGAGCTGGGTTATGCAGACGTGGGTGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGAttt                                                                                   <  1:63037/249‑1 (MQ=255)                                                                                                                                                    aaGAGATGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAACAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGGGTAACGTTGGATCTGCCAAAAGAGCTGGAAGAACGTGCTGATTTAccc                                                                               >  1:104402/1‑251 (MQ=255)                                                                                                                                                                atgcagcagtgtctggtgacaCCATGCAACATAAAATGCGTGAGTTAACAGAACTGGGCTACACATCGAAGGCGCACGATTCCGCAGCATATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGCATATGGCTTGCTCCTCGGCACC‑GCGCAGCAGGCGATATTTACAGAGAGCGTGGAAACGCTTCCCGGTGTATCGTTGGATGTGCCAAACGACCTGGCAGAACGTGCTGATTTACCCTTGCtttacaat                                                                   >  2:32161/21‑246 (MQ=255)                                                                                                                                                                         gggtgtatgaAACCACGCAACAGACAATGCGTGAGCTCAAAGACCTGGCATACACCACAGAGCCGCCCGCTGCTGTAGCGTATCGTGC‑GCTGCGTGACAAGTTGAAATTAGGCGCATATGGCGTGTTCCTCGTCCCC‑GCTCATCCGGCGAAATATAAAGAGAGCGTGTAAGCGAGTCTCGGTGAAACGTTGGAGCTGCCAAAAGAGCTGGCAGAACGTGCTGATGTTCCCTTTCTGTCACATAATCCGccc                                                          <  2:295107/245‑1 (MQ=255)                                                                                                                                                                               atgaAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCccgccggt                                                    >  1:331417/1‑249 (MQ=255)                                                                                                                                                                                              aGACAATGCGTGAGATAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgcgt                                     >  1:303320/1‑251 (MQ=255)                                                                                                                                                                                                aCAATGTGTGATTTAAAATAACAGGGCTACACCTCGGAGCCTCGCTCTGCCTTAGTGTGTCGTTC‑GCTGCGTGACTAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACC‑TCGCATCCGGCTAAATGTAAAGAGAGCGTGGAATCGATTCTCGGTGAAACTTTGGATCTGCAAAAAGAGCTGGCAGAACGTGCTTTTTTCCCCTTGTTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTa                                    <  1:215049/250‑1 (MQ=255)                                                                                                                                                                                                  aaTGCGTGAGGAAAAAGAACTGGGCTAAACCGCGGAGCGGCACGCTGCCGTAGCGTATCGGGG‑GATGGGGGACGAGTTGAATCCAGGCGAATAGGGCGTGTTCCTGGGCACG‑GGGCAAACGGAGAAATTTAAAGAGAGCGGGGAAGCGATTCGCGGGGAAACGTGGAATATGCAAAAAGAGAGGGCAGAACGGGCTGATTTCCCCTTCCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAAt                                 <  1:68793/251‑1 (MQ=255)                                                                                                                                                                                                              aaaaGAACTGGGCTACGCCTCGGAGCCGCACGCGGCCGTAGCGTTTCGTGC‑GCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGAAGAACGTGGTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatc                     <  1:34778/251‑1 (MQ=255)                                                                                                                                                                                                                  aaaCTGGGCTACTCTTGGGAGCCGCAGGCTGCGGTAGCGTTTC‑TGCCGCTGCGTCACCAGTTGACTCCAGGCGAATAGGGCTTGTTCCTCGGCCCC‑GCGCAGCCGGAAAAATTAAAAGAAAGTGTGGAAGCGATTCTCGGTGAAACGTTGGATCTCCCAAAAGAGCTGGCAGAACGTGCTGATTTCCCCTTTCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAAtcatca                 <  1:307314/250‑1 (MQ=255)                                                                                                                                                                                                                       gggCTACACCCCGGATCAGCACGGTGCCGTAGCGTACCGTGC‑GCTGCGTAACCAGTTGAAACCAGGCGAATATGGCTTGTTCCCCGGCACC‑GCGCTTCCGGCGAAATTAAAAGATAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGATCTGGCAGACCGGGCTAATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATTATGAATCATCAGTaa             <  1:338652/250‑1 (MQ=255)                                                                                                                                                                                                                                 ccGGAGCCGCACGCTGCTGTGGCGTATCGTGC‑GTTGCTGGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCAt  <  1:370674/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      GAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCG‑TGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGC‑GCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACC‑GCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCAT  >  NC_000913/4558‑5031

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