Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,744,362	G→A	H104Y (CAC→TAC)	apt ←	Adenine phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,744,362	0	G	A	96.8%	91.5 / NA	31	H104Y (CAC→TAC)	apt	Adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/16);  minor base C (0/1);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTTCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAC  >  USA300TCH1516_ALE/1744217‑1744499                                                                                                                                                  |                                                                                                                                          tatTTCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACa                                                                                                                                         >  1:1080920/1‑148 (MQ=255)   tttCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAtt                                                                                                                                       >  1:730553/1‑148 (MQ=255)   tttCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAtt                                                                                                                                       >  1:41627/1‑148 (MQ=255)              ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAtt                                                                                                                            <  1:262263/148‑1 (MQ=255)              ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGCACATTGTTAAAACAtt                                                                                                                            <  1:32809/148‑1 (MQ=255)              ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTGTAATTGCATCTTTGTACATTGTTAAAACAtt                                                                                                                            <  1:497499/148‑1 (MQ=255)              ttATAAATGCAATACCTACTACGATACCGCCTAATCTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAtt                                                                                                                            <  1:393229/148‑1 (MQ=255)                            ccTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTa                                                                                                              <  1:313355/148‑1 (MQ=255)                               actacGATACCGCCTAATTTTTCCACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGt                                                                                                           <  1:214964/148‑1 (MQ=255)                                    gATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATAc                                                                                                      >  1:688478/1‑148 (MQ=255)                                       tccGCCCAAGTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCCCTACCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCa                                                                                                   <  1:974413/147‑1 (MQ=255)                                           ccTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCCCCAGTAGCTAATAAAACATCTGTAATTAACACACGTTGACCTGGTTTACTTGCATCTTTGTACATTGTTAAAACATTTGTCCCATATTCTAGGTCATACTCATAAc                                                                                               >  1:863736/1‑148 (MQ=255)                                                    tCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAAGATTTGTACCATATTCTAGGTGCTACTCATCACGCATGACtt                                                                                      >  1:55549/1‑148 (MQ=255)                                                                          cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACCCGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                                <  1:52685/148‑1 (MQ=255)                                                                          cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCttt                                                                <  1:723855/148‑1 (MQ=255)                                                                                        gTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaa                                                  <  1:430885/148‑1 (MQ=255)                                                                                         tAGCTAATAAATCATCTGTAATTAACACGCGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCaaa                                                 >  1:643256/1‑148 (MQ=255)                                                                                         tAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGACTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCcaa                                                 >  1:839435/1‑148 (MQ=255)                                                                                                  aaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTCCCCCGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATccc                                        <  1:1106222/148‑1 (MQ=255)                                                                                                   atcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                       <  1:117926/148‑1 (MQ=255)                                                                                                   agcatcTTGAATTAACACACGTTGCCCTGGTTTAATTGCATCTTTGTACATTGTCAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                       <  1:1176917/146‑1 (MQ=255)                                                                                                                  cacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCTCGATTGAATAAGCTACa                        >  1:1150222/1‑148 (MQ=255)                                                                                                                        ttGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAg                  >  1:1024923/1‑148 (MQ=255)                                                                                                                        ttGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAg                  <  1:730534/148‑1 (MQ=255)                                                                                                                        ttGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAg                  <  1:436623/148‑1 (MQ=255)                                                                                                                           aCCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCa               <  1:11861/148‑1 (MQ=255)                                                                                                                             cTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAt             <  1:85825/148‑1 (MQ=255)                                                                                                                                    aaTTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTACCAGGTGCAAACCCAATCCCCATTTAATAAGCTACAGGACAACCAATGGTAAAg      >  1:45336/1‑148 (MQ=255)                                                                                                                                       tGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGccc   >  1:1097022/1‑147 (MQ=255)                                                                                                                                        gCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGCGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACGGGACAGCCAATGATACAGCCAc  >  1:121358/1‑148 (MQ=255)                                                                                                                                        gCATCTTTGTACATTGTTAAAACATTTGGACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAc  >  1:653786/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                          TATTTCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAC  >  USA300TCH1516_ALE/1744217‑1744499

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