Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BX571857	2,727,708	G→A	22.3%	T430I (ACC→ATC)	SAS2534 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BX571857	2,727,708	0	G	A	22.3%	75.5 / 15.9	31	T430I (ACC→ATC)	SAS2534	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/11);  new base A (3/4);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCAG  >  BX571857/2727614‑2727804                                                                                               |                                                                                                 aaTGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTg                                                                                               >  1:75653/1‑98 (MQ=255)  aTGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAGGGTTgg                                                                                              <  2:95186/98‑1 (MQ=255)  aTGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTgg                                                                                              <  1:108635/98‑1 (MQ=255)  aTGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTgg                                                                                              <  2:278418/98‑1 (MQ=255)   tGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTggg                                                                                             >  2:141227/1‑98 (MQ=255)   tGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTggg                                                                                             <  2:89811/98‑1 (MQ=255)      cATGAGCTAGTTTTGCGATAACATCATCCACTCTATCTGCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTa                                                                                          <  1:207522/98‑1 (MQ=255)      cATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTa                                                                                          >  2:295973/1‑98 (MQ=255)         gAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCa                                                                                       >  2:103868/1‑98 (MQ=255)         gAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGCTACCa                                                                                       >  2:4711/1‑98 (MQ=255)             tAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAAt                                                                                   <  1:72195/98‑1 (MQ=255)               gTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTc                                                                                 <  2:242993/98‑1 (MQ=255)                 tttGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTcac                                                                               <  2:147482/98‑1 (MQ=255)                  ttGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTcaca                                                                              <  2:142306/98‑1 (MQ=255)                           catcatCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCa                                                                     >  2:301146/1‑98 (MQ=255)                                          cTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCt                                                      >  1:103676/1‑98 (MQ=255)                                          cTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCt                                                      >  1:34239/1‑98 (MQ=255)                                              aTCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGc                                                  >  1:236393/1‑98 (MQ=255)                                                   tATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATACTTCACATTTAATCCAATCATACCTAGTCCTGAGCCCAcc                                             >  2:70437/1‑98 (MQ=255)                                                      cAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCt                                          >  2:58873/1‑98 (MQ=255)                                                         aaGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTcc                                       >  1:35135/1‑98 (MQ=255)                                                         aaGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTcc                                       <  1:156055/98‑1 (MQ=255)                                                          aGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCa                                      <  1:49361/98‑1 (MQ=255)                                                              aCCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAg                                  >  1:281431/1‑98 (MQ=255)                                                              aCCAAGTAGCCATTTTCCCCATCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAg                                  >  2:90174/1‑98 (MQ=255)                                                                          ttttCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCaaa                      <  1:143829/98‑1 (MQ=255)                                                                                  aTCTCGAATAAAGATTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGc              <  2:46204/98‑1 (MQ=255)                                                                                        aaTAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTg        >  2:188465/1‑98 (MQ=255)                                                                                        aaTAAAGATTGGGTTACCATAATTCACATTTAGTCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTg        <  1:153871/98‑1 (MQ=255)                                                                                          tAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCCTCGCTAGTCGAt      >  2:197458/1‑98 (MQ=255)                                                                                            aaGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGc    <  1:71113/98‑1 (MQ=255)                                                                                              ggTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCag  <  2:69667/98‑1 (MQ=255)                                                                                               |                                                                                                 AATGGCATGAGCTAGTTTTGCGATAACATCATCCACACTATCTTCATCGGTATCAAAAGGTACCAAGTAGCCATTTTCCCCATCTCGAATAAAGGTTGGGTTACCATAATTCACATTTAATCCAATCATACCTAGTCCTGAGCCTACCGCTTCCATTAGTGTTAACCCAAAACCTTCGCTAGTTGATGCAG  >  BX571857/2727614‑2727804

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