Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	729,046	G→T	5.6%	M305I (ATG→ATT)	USA300HOU_RS03535 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	729,046	0	G	T	5.6%	269.8 / 2.4	87	M305I (ATG→ATT)	USA300HOU_RS03535	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/35);  new base T (3/2);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCTTTAATAGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGATA  >  USA300TCH1516_ALE/728912‑729182                                                                                                                                       |                                                                                                                                         ccTTTAATAGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTa                                                                                                                                    <  2:80717/141‑1 (MQ=255)    ttAATAGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCa                                                                                                                                 >  1:606181/1‑141 (MQ=255)     tAATAGCAACAGTGTTATCTGTTATTATGGTGGTGGTTTTGTTATACACTATTTTAGGAACCATAAACAGATTTTCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATTATGGTAGCAg                                                                                                                                >  1:1024539/1‑141 (MQ=255)     tAATAGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATTATGGTAGCAg                                                                                                                                >  1:428281/1‑141 (MQ=255)         aGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGtta                                                                                                                            >  2:595471/1‑141 (MQ=255)            aaCAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAAAGATGGTAGCAGGttattt                                                                                                                         >  2:807061/1‑141 (MQ=255)            aaCAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGAACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGttattt                                                                                                                         <  2:35845/141‑1 (MQ=255)               aGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATAGTATTCTAGGACTAATGTCTTCATGTGCAGCACGTTTTACAGAACTATACAGTAAAAATTATCATAAATTTAATATTAAAATTATGGTAGCAAgttagttact                                                                                                                      >  2:254000/1‑139 (MQ=255)                gtgtTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATttatta                                                                                                                     >  1:880670/1‑141 (MQ=255)                gtgtTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTAGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGTAGGTTATttatta                                                                                                                     >  2:962374/1‑141 (MQ=255)                  gtTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCCTACAGTAAAAATTATCATATCTTTATTCTTATGAACAGGTTAAAAGGGTATTTATTAAg                                                                                                                   >  1:28801/1‑141 (MQ=255)                   ttATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGt                                                                                                                  >  1:923680/1‑141 (MQ=255)                      tCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACGGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTc                                                                                                               >  1:773609/1‑141 (MQ=255)                      tCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACGGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTc                                                                                                               >  1:569535/1‑141 (MQ=255)                       cTGTTATTATGCTGGCGGTTATGTATAATACTATTCTGGGACTAATGTATTCATTTGGCACAGGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCg                                                                                                              <  1:1018395/141‑1 (MQ=255)                          ttattaTGTTGGCGGTTTTGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAg                                                                                                           <  2:73452/141‑1 (MQ=255)                              tatGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGAtt                                                                                                       >  1:590579/1‑141 (MQ=255)                                 gCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAATTTTAGTAGGATTTGc                                                                                                    >  2:998289/1‑141 (MQ=255)                                 gCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTAGTCATTTACGGCACGTTTTACAACCACATCGGTTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGc                                                                                                    <  1:915780/141‑1 (MQ=255)                                  cTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCt                                                                                                   <  2:274006/141‑1 (MQ=255)                                  cTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGCTTATTTAGTAAGTTTCGTAAGATTTGCt                                                                                                   <  2:178475/141‑1 (MQ=255)                                   tGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGCACCATACAGTAAAAATTCTCATATCTTTATTAGTATAATGACGGGAGCAGGTTATTTATTAAGTTTCGTAGCATTTGCTg                                                                                                  >  2:416764/1‑141 (MQ=255)                                   tGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTg                                                                                                  >  2:298708/1‑141 (MQ=255)                                     gctgTTATGTAGAATTATAGTCTATGACTAATATATTCATTTTTAGTACGTTTTACAGAACCATACAGTAAAAATTATCTTATCTTTATTATTATAATGATGCTAGGAGGTTATTTATTAAGTTTCGTAGGATTTGCTGaa                                                                                                <  2:536321/138‑1 (MQ=255)                                       ggTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGaatt                                                                                              <  2:185433/141‑1 (MQ=255)                                                    aCTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTta                                                                                 <  1:591377/141‑1 (MQ=255)                                                         tCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATAc                                                                            <  2:606181/141‑1 (MQ=255)                                                          cTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATTATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACa                                                                           <  2:428281/141‑1 (MQ=255)                                                          cTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACa                                                                           <  2:771827/141‑1 (MQ=255)                                                          cTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACa                                                                           >  2:45413/1‑141 (MQ=255)                                                          cTAGGACTAATGTATTCATTTGCAGCACGGTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTTAATTAATTAATAAGTtatatata                                                                           >  2:759099/1‑139 (MQ=255)                                                          cTAGGACTAATGTACTCATTTGCAGCACGTTTTACAGAATCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACa                                                                           <  1:191619/141‑1 (MQ=255)                                                          cTAGAACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATTATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACa                                                                           <  2:1024539/141‑1 (MQ=255)                                                            aGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAAt                                                                         >  1:824186/1‑141 (MQ=255)                                                            aGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAAt                                                                         >  2:875280/1‑141 (MQ=255)                                                              gagtAAGGTATTCATTTGCATGACGTTTTACAGAACGATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTa                                                                       <  2:749823/138‑1 (MQ=255)                                                              gACTAATGTATTCATTTGCAGTACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTa                                                                       <  2:645423/141‑1 (MQ=255)                                                               aCTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTAt                                                                      <  1:800279/141‑1 (MQ=255)                                                                  aaTGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATggg                                                                   <  1:183246/141‑1 (MQ=255)                                                                  aaTGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATggg                                                                   <  1:264476/141‑1 (MQ=255)                                                                   aTGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTATGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGa                                                                  <  1:82927/141‑1 (MQ=255)                                                                   aTGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGa                                                                  <  1:936604/141‑1 (MQ=255)                                                                       atttatttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATg                                                              <  2:826807/141‑1 (MQ=255)                                                                       aTTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATg                                                              <  2:79632/141‑1 (MQ=255)                                                                        ttCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGt                                                             >  1:569727/1‑141 (MQ=255)                                                                            tttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGc                                                         >  1:506025/1‑141 (MQ=255)                                                                            tttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGc                                                         >  1:951877/1‑141 (MQ=255)                                                                            tttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGc                                                         <  1:304056/141‑1 (MQ=255)                                                                            tgtGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGc                                                         <  1:64024/139‑1 (MQ=255)                                                                             ttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCt                                                        <  2:923680/141‑1 (MQ=255)                                                                             ttGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCt                                                        >  2:651758/1‑141 (MQ=255)                                                                             ttGCAGCACGTTTTACAGAACCATACAGTAAAAAGTATCATATCTTTGTTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCt                                                        <  2:924025/141‑1 (MQ=255)                                                                              tGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCtt                                                       >  2:977897/1‑141 (MQ=255)                                                                                cagcaCGTTTTACAGAACCATACCGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttat                                                     <  1:1012656/141‑1 (MQ=255)                                                                                cagcaCGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttat                                                     >  1:900191/1‑141 (MQ=255)                                                                                cagcaCGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttat                                                     <  1:363061/141‑1 (MQ=255)                                                                                    aCGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttatttat                                                 >  1:580212/1‑141 (MQ=255)                                                                                    aCGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttatttat                                                 >  1:582681/1‑141 (MQ=255)                                                                                    aCGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTAATATAATGATGGTAGCAGATTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCttatttat                                                 >  1:865690/1‑141 (MQ=255)                                                                                         ttACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTgtat                                            >  2:518812/1‑140 (MQ=255)                                                                                         ttACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGCAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGCCTTATGTATTgtac                                            >  2:153907/1‑140 (MQ=255)                                                                                         ttACAGAACCATACAGTAAAAAGTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTgtag                                            <  2:438057/141‑1 (MQ=255)                                                                                           aCAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTgtagta                                          >  2:381380/1‑141 (MQ=255)                                                                                                aCCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTACTTACTAGGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGt                                     >  1:334576/1‑141 (MQ=255)                                                                                                aCCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGCTAGGTAGGCTTATTTATCGTAGTAGCTGt                                     >  1:434393/1‑141 (MQ=255)                                                                                                aCCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGt                                     >  1:558907/1‑141 (MQ=255)                                                                                                aCCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGt                                     >  1:331779/1‑141 (MQ=255)                                                                                                aCCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTAGTTATTAAGGTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATAGAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGt                                     >  1:806421/1‑141 (MQ=255)                                                                                                           aaaaTTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAat                          >  1:280753/1‑141 (MQ=255)                                                                                                             aaTTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAatat                        >  2:244802/1‑141 (MQ=255)                                                                                                             aaTTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAatat                        >  1:153201/1‑141 (MQ=255)                                                                                                             aaTTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAatat                        >  1:153157/1‑141 (MQ=255)                                                                                                             aaTTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGATGTAATTATTAAatat                        >  2:692479/1‑141 (MQ=255)                                                                                                                 aTCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCa                    <  1:91522/141‑1 (MQ=255)                                                                                                                  tCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTCTTGAAGTAGCTGTAATTATTAAATATTTCaa                   >  2:500352/1‑141 (MQ=255)                                                                                                                  tCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTTTAGTAGCTGTAATGGTTAAttttttgta                   >  1:715332/1‑138 (MQ=255)                                                                                                                      atttttATTATTATAATGATGGTAGCATGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGt               <  2:674522/138‑1 (MQ=255)                                                                                                                      atCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGt               <  1:580070/141‑1 (MQ=255)                                                                                                                        cTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTaa             <  1:1011310/141‑1 (MQ=255)                                                                                                                           tattattatAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTGATTATTAAATATTTCAAACGTaaaaa          >  1:765768/1‑141 (MQ=255)                                                                                                                           tattattatAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTAAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTaaaaa          >  1:149166/1‑141 (MQ=255)                                                                                                                              tattatAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACTTAAAATTGc       >  2:593177/1‑141 (MQ=255)                                                                                                                                 ttgaaTGATGATAGAATGTTATTTATTATGTTTCGTAGGATTTGCTGAATTAATTAATTAGTTATATACAATTATGGGATATGTAGTCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGa    <  2:58117/138‑1 (MQ=255)                                                                                                                                  ataATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGAt   >  2:864671/1‑141 (MQ=255)                                                                                                                                  ataATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGAt   <  2:433487/141‑1 (MQ=255)                                                                                                                                  ataATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGAt   >  1:68194/1‑141 (MQ=255)                                                                                                                                   taATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGATa  >  1:430917/1‑141 (MQ=255)                                                                                                                                       |                                                                                                                                         CCTTTAATAGCAACAGTGTTATCTGTTATTATGCTGGCGGTTATGTATAATACTATTCTAGGACTAATGTATTCATTTGCAGCACGTTTTACAGAACCATACAGTAAAAATTATCATATCTTTATTATTATAATGATGGTAGCAGGTTATTTATTAAGTTTCGTAGGATTTGCTGAATTAATTAATAAGTTATATACAATTATGGGATATGTAGGCTTATTTATTGTAGTAGCTGTAATTATTAAATATTTCAAACGTAAAAATGCGGATA  >  USA300TCH1516_ALE/728912‑729182

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