Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	463,014	C→A	E47* (GAA→TAA)	USA300HOU_RS02190 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	463,014	0	C	A	100.0%	129.8 / NA	36	E47* (GAA→TAA)	USA300HOU_RS02190	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/19);  total (17/19)

GTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCATATCTATC  >  USA300TCH1516_ALE/462885‑463159                                                                                                                                  |                                                                                                                                                  gTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATg                                                                                                                              >  1:221923/1‑151 (MQ=255)      gCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATgtagta                                                                                                                         >  1:283842/1‑151 (MQ=255)      gCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTATGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATgtagta                                                                                                                         >  1:284037/1‑151 (MQ=255)           gttgATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGt                                                                                                                    >  1:310918/1‑151 (MQ=255)                 aaTGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAAATTGCATCATCATGTCGTACGAGTCCAttt                                                                                                              >  2:295638/1‑151 (MQ=255)                         tgattgaATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATa                                                                                                      >  1:104802/1‑151 (MQ=255)                           attgaATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATaaa                                                                                                    <  1:248415/151‑1 (MQ=255)                                  gTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGc                                                                                             >  1:192904/1‑151 (MQ=255)                                            aTCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGa                                                                                   >  1:294089/1‑151 (MQ=255)                                             tCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGAt                                                                                  <  1:182581/151‑1 (MQ=255)                                                gCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTAc                                                                               <  2:112943/151‑1 (MQ=255)                                                 cGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTAcc                                                                              <  2:259697/151‑1 (MQ=255)                                                      tATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATc                                                                         >  1:186039/1‑151 (MQ=255)                                                              tAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACg                                                                 <  2:51787/151‑1 (MQ=255)                                                                    gTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGtt                                                           <  1:228369/151‑1 (MQ=255)                                                                     tCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTg                                                          <  2:192104/151‑1 (MQ=255)                                                                      cTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGTATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGc                                                         >  2:179240/1‑151 (MQ=255)                                                                      cTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGc                                                         >  1:295805/1‑151 (MQ=255)                                                                          tGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTaaa                                                     >  2:106890/1‑151 (MQ=255)                                                                           gTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTaaaa                                                    >  1:187094/1‑151 (MQ=255)                                                                            tAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAac                                                   <  1:117921/151‑1 (MQ=255)                                                                               cTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACAcc                                                <  1:1540/151‑1 (MQ=255)                                                                                   gTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTAcc                                            <  1:139249/151‑1 (MQ=255)                                                                                        cATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGc                                       <  1:191192/151‑1 (MQ=255)                                                                                            cTTATCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTa                                   >  1:85302/1‑151 (MQ=255)                                                                                                tCTAATAACACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGcta                               <  1:112955/151‑1 (MQ=255)                                                                                                  taataaCACATCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGctact                             <  2:26103/151‑1 (MQ=255)                                                                                                         acaTCAATATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCt                      <  2:192904/151‑1 (MQ=255)                                                                                                              aaTATAACTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTa                 >  2:237117/1‑151 (MQ=255)                                                                                                                    aCTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATATGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTc           >  2:44639/1‑151 (MQ=255)                                                                                                                     cTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCa          <  1:148289/151‑1 (MQ=255)                                                                                                                     cTTAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCa          <  1:123129/151‑1 (MQ=255)                                                                                                                      ttAATTCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCat         <  2:22344/151‑1 (MQ=255)                                                                                                                          ttCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCAtatct     <  2:46291/151‑1 (MQ=255)                                                                                                                          ttCATGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCAtatct     >  2:87028/1‑151 (MQ=255)                                                                                                                             aTGTTATAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAATACCTTTTTAAATTTCAtatctatc  <  2:142500/151‑1 (MQ=255)                                                                                                                                  |                                                                                                                                                  GTTTAGCTTTGTTGATAATGTCTTTGATTGAATGTAGTCCTTTATCGGCGAAGTATGATCTTAAGTTGTCTTTTGTAGCTTGGTCAGCATTCTTATCTAATAACACATCAATATAACTTAATTCATGTTCTAAGAAGTTTGCATCATCATGTAGTACGAGTCCATTTTGAGAATAAACTTTCGCATCTGCTTGATTACCATATCCAACAACGCCAGTTGCTAAAACACCTACCATTGCCGTAGCTACTAAAACCTTTTTAAATTTCATATCTATC  >  USA300TCH1516_ALE/462885‑463159

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