Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	848,832	T→G	23.3%	G227G (GGT→GGG)	SAS0756 →	thermonuclease precursor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	848,832	0	T	G	23.3%	62.5 / 24.4	30	G227G (GGT→GGG)	SAS0756	thermonuclease precursor
Reads supporting (aligned to +/- strand):  ref base T (7/16);  new base G (2/5);  total (9/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

CTAACAATACACATGAACAACTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATTTATTATCAATTCAGAATTAAA  >  BX571857/848738‑848924                                                                                               |                                                                                             cTAACAATACACATGAACAACTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCa                                                                                            <  1:264820/97‑1 (MQ=255)                aacaacTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTa                                                                             <  1:279015/97‑1 (MQ=255)                  caacTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTaaa                                                                           >  2:236793/1‑97 (MQ=255)                   aacTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTaaaa                                                                          <  1:76280/97‑1 (MQ=255)                     cTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGACCATTGTAAAAgt                                                                        >  1:94864/1‑97 (MQ=255)                      ttttAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAgtg                                                                       >  1:235397/1‑97 (MQ=255)                       tttAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAgtgt                                                                      <  1:215089/97‑1 (MQ=255)                        ttAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTc                                                                     >  1:321155/1‑97 (MQ=255)                           aGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAGTGTCAct                                                                  <  1:158050/97‑1 (MQ=255)                                aaGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAGTGTCACtgctgc                                                             <  2:178041/97‑1 (MQ=255)                                  gTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGGTTCAGGGCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTa                                                           <  2:332005/97‑1 (MQ=255)                                   tGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAg                                                          >  1:161527/1‑97 (MQ=255)                                    gAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGt                                                         <  2:258057/97‑1 (MQ=255)                                       gCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGc                                                      >  2:351135/1‑97 (MQ=255)                                            aGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACttt                                                 <  2:252694/97‑1 (MQ=255)                                                 aaaaaGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATaa                                            <  1:55774/97‑1 (MQ=255)                                                    aaGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAAttt                                         >  1:362993/1‑97 (MQ=255)                                                      gagaAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAAttttt                                       <  1:312151/97‑1 (MQ=255)                                                           aTTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATc                                  <  1:111682/97‑1 (MQ=255)                                                             tAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCAc                                <  1:339761/97‑1 (MQ=255)                                                             tAAATATTTGGAGCGAAGACAACGCTGATTCAGGGCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCAc                                <  1:230561/97‑1 (MQ=255)                                                               aaTATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGa                              <  2:282219/97‑1 (MQ=255)                                                                  atTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGatat                           >  1:230999/1‑97 (MQ=255)                                                                       gAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGAttt                      <  1:303437/97‑1 (MQ=255)                                                                         gCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATttat                    <  2:184330/97‑1 (MQ=255)                                                                                 aaCGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATTTATTATCAAtt            <  2:236034/97‑1 (MQ=255)                                                                                 aaCGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATTTATTATCAAtt            <  1:163335/97‑1 (MQ=255)                                                                                        aTTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCGCTTTTATAATTTTTGGATCACGATATGATTTATTATCAATTCAGAAtt     >  2:346608/1‑97 (MQ=255)                                                                                           cAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGCTCACGATATGATTTATTATCAATTCAGAATTaaa  <  2:319197/97‑1 (MQ=255)                                                                                           cAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATTTATTATCAATTCAGAATTaaa  <  2:319388/97‑1 (MQ=255)                                                                                               |                                                                                             CTAACAATACACATGAACAACTTTTAAGAAAAAGTGAAGCACAAGCGAAAAAAGAGAAATTAAATATTTGGAGCGAAGACAACGCTGATTCAGGTCAATAATGCTCATTGTAAAAGTGTCACTGCTGCTAGTGGCACTTTTATAATTTTTGGATCACGATATGATTTATTATCAATTCAGAATTAAA  >  BX571857/848738‑848924

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