Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000964	490,812	G→T	P78T (CCT→ACT)	ydzK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000964	490,812	0	G	T	95.0%	63.2 / NA	21	P78T (CCT→ACT)	ydzK	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (12/7);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GGCCCTAAGAGCCCTTGTTTTTTTTTTTTTTTTTGGAAACAAATGTATTTTCAGCATGTTTGAATTCATAAAATGATCGGTATATGTATTAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGGAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATCCACAAATAAAACAGAATAAAT  >  NC_000964/490564‑491035                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                ggCCCTAAGAGCCCTTGTTTTTTTTTTTTTTTTTGGAAACAAATGTATTTTCAGCATGTTTGAATTCATAAAATGATCGGTATATGTATTAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGGaa                                                                                                                                                                                                                               <  2:314354/251‑1 (MQ=255)                                   ggAAACAAATGTATTTTCAGCATGTTTGAATTCATAAAATGATCGGTATATGTATTAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATaaaa                                                                                                                                                                                              >  1:52551/1‑250 (MQ=255)                                                                   cATAAAAGGATCGGGTTATGTATTAGGAACTGATTGCGCCATACTCGAAACAGACAAAAGGACAAAGGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCATAACCGACCAAGAAGGTGCCGAATAAAAACATTCCGCCGATAAACCCGTGGAAATTtgtg                                                                                                                                                              <  2:271067/250‑1 (MQ=255)                                                                           gATCGGTATATGTATTAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGcc                                                                                                                                                     >  1:304442/1‑251 (MQ=255)                                                                               ggTATATGTAATAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGTTGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACGCTTCACATAATAGGTTATTTTTCTTTTTTCTTAAGCCCGAAGAATAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAACCATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTaa                                                                                                                                                 <  2:4014/251‑1 (MQ=255)                                                                                            agAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATATTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCATATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTACTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACTCAGAAGGTGCCGAATACAAACATGACGCCGGTAAACACGTGGAAGTTTGTTTGTCCGGCAGTAAACAACGTGAtt                                                                                                                                      >  2:36464/2‑249 (MQ=255)                                                                                                     ttCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGAt                                                                                                                           >  2:52721/1‑251 (MQ=255)                                                                                                                cGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACg                                                                                                                 >  2:290844/1‑250 (MQ=255)                                                                                                                                               ttCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCAttt                                                                                  <  1:313832/250‑1 (MQ=255)                                                                                                                                                            gCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAAAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGAttt                                                                     <  2:147113/250‑1 (MQ=255)                                                                                                                                                                        gCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAAc                                                         <  1:116417/250‑1 (MQ=255)                                                                                                                                                                           tGCTTATTCTTATAAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGACGCGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGa                                                      >  2:244064/1‑250 (MQ=255)                                                                                                                                                                                       tCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGACGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAACACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCg                                         >  1:334936/1‑251 (MQ=255)                                                                                                                                                                                             cTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAAc                                    >  2:277374/1‑250 (MQ=255)                                                                                                                                                                                                aCAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAAc                                                    >  2:187540/1‑231 (MQ=255)                                                                                                                                                                                                aCAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAAc                                                    <  1:187540/231‑1 (MQ=255)                                                                                                                                                                                                        attctCATAACAGGTAATTTTACTTTTTTCTTCCAACAGAAGGAGAAAGCAAGACCCGAAGCGACCAAGAAGATCCCGAATGAAAACATGCCGCGGACAAACACGCGGAAATTTCTGTGTCCGGCCGTCCACGCCGTCTTTAGGGAAATGATACTGAGCACGCTCCTGGCAATGCCCGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGCTAACAACAGGAACaa                         >  2:277558/2‑250 (MQ=255)                                                                                                                                                                                                             cATAATAGGTTATTTTGCTTTTTACTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATcc                      <  2:170158/248‑1 (MQ=255)                                                                                                                                                                                                                aataGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATCCACAAAt                >  1:76737/1‑251 (MQ=255)                                                                                                                                                                                                                   aGGTTATTTTTCTTTTTTCTTCAGCCCGAAGCAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGGGTCCGGCCGTAAACGCCGGGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATCCACAAATaa              >  1:1565/1‑250 (MQ=255)                                                                                                                                                                                                                               ttttttcttCAGCCCGAAGAAGAAAGTAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATCCACAAATAAAACAGaataaat  >  2:49453/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                GGCCCTAAGAGCCCTTGTTTTTTTTTTTTTTTTTGGAAACAAATGTATTTTCAGCATGTTTGAATTCATAAAATGATCGGTATATGTATTAGGAACTGATTTCGCCATACTCGAAACAGACAAAAGGACAAATGCATCTTGTTTCAAAGGAAGATGCATTTGTCCTTGCATGCTTATTCTTATCAAATCTTACAGAACTCTTCTCATAATAGGTTATTTTTCTTTTTTCTTCAGCCCGAAGAAGAAAGGAAAAGCAGAACCGACCAAGAAGGTGCCGAATAAAAACATGCCGCCGATAAACACGTGGAAATTTGTGTGTCCGGCCGTAAACGCCGTGATTAGGGAAATGATACTGAGCACGATCATGGCAATGCCGGGCGCGACTCTCATTTTTCCCTTGGATTTATGGGCAGCAACAGAACATAAAATGCCGATAACAATAGGAACAATCCACAAATAAAACAGAATAAAT  >  NC_000964/490564‑491035

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