Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	111,516	C→A	37.7%	Q196K (CAA→AAA)	USA300HOU_RS00520 →	putative lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	111,516	0	C	A	37.7%	63.0 / 56.6	68	Q196K (CAA→AAA)	USA300HOU_RS00520	putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base A (13/13);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111376‑111651                                                                                                                                             |                                                                                                                                        cTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTc                                                                                                                                         <  1:295196/141‑1 (MQ=255)  tGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCa                                                                                                                                        <  2:185421/141‑1 (MQ=255)    aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGTTGCCGAGTTATTCAGCTAAat                                                                                                                                      >  2:43844/1‑141 (MQ=255)    aGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGGTATTCGGCTAAat                                                                                                                                      >  1:489664/1‑141 (MQ=255)     ggAAGTTAAAGAAGAAAAAATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAata                                                                                                                                     <  2:387382/141‑1 (MQ=255)          ttAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAAt                                                                                                                                <  1:366283/141‑1 (MQ=255)          ttAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAAAATAGCAAGTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGAAGAGTTATGCGGCTAAATATCAAt                                                                                                                                >  2:257732/1‑141 (MQ=18)            aaaGAAGAAAACATAAAAGAAGAAATCGAAAACGTTAAGTTTTTTGCGCACTCTAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                              <  1:305861/141‑1 (MQ=18)            aaaGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTa                                                                                                                              <  1:249116/141‑1 (MQ=255)             aaaaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTaa                                                                                                                             <  2:91232/138‑1 (MQ=255)              agaagaAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAAc                                                                                                                            >  1:483141/1‑141 (MQ=255)               gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCATAGGTGCCGAGTTATTCAGCTCAATATCAATTAACt                                                                                                                           >  2:5694/1‑141 (MQ=255)               gaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                           >  2:69423/1‑141 (MQ=255)               aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                           <  1:270968/140‑1 (MQ=255)               aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                           <  1:108941/140‑1 (MQ=255)               aaagaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATGTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACt                                                                                                                           <  2:191095/140‑1 (MQ=255)                  gaaAACATAAAAAAAGAAATTGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACtcaa                                                                                                                        <  2:464333/141‑4 (MQ=255)                  gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTCAAGATTGATTTAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                        <  2:301658/141‑1 (MQ=255)                  gaaAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAat                                                                                                                        >  1:45872/1‑141 (MQ=255)                  gaaAACATAAAAAAAGAAATCGAAAACCTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATAAAGTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAat                                                                                                                        <  1:314474/141‑1 (MQ=17)                     aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                     >  2:123979/1‑141 (MQ=255)                     aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgat                                                                                                                     >  2:401533/1‑141 (MQ=255)                     aaCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTAATTCGGCTAAATATCAATTAACTAatgat                                                                                                                     >  1:304861/1‑141 (MQ=255)                     aaCATAAAAAAAGAAATCGAAAACTTGAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAatgat                                                                                                                     >  2:471721/1‑141 (MQ=255)                      aCATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAatgatg                                                                                                                    <  1:406341/141‑1 (MQ=255)                           aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTACTGCATTAATGAATTATAAAGCTGCAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                               <  1:342025/141‑1 (MQ=21)                           aaaaaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGATAAATATCAATTAACTAATGATGATtat                                                                                                               >  2:158623/1‑141 (MQ=255)                           aaaaaaGAAATCGAAAACTTTAAGTTTGGTGCGCAATATAGCAATTTTAAAGATTTAATGAAGTATAAAGATGGATATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATtat                                                                                                               <  2:464786/141‑1 (MQ=255)                              aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                            >  2:205508/1‑141 (MQ=255)                              aaaGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAAt                                                                                                            >  1:168603/1‑141 (MQ=255)                                      cGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACa                                                                                                    >  1:378729/1‑141 (MQ=255)                                       gAAAACTGTAAGTTCTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATAACAATTAACTCATGATGATTCTAATGTAAAACaa                                                                                                   >  2:24711/1‑141 (MQ=17)                                          aaCTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTa                                                                                                >  2:474265/1‑141 (MQ=255)                                          aaCTTTAAGTTGTTTGCGAAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATGAACTAATGATGATTATAATGTAAAACAATTa                                                                                                <  2:346435/141‑1 (MQ=255)                                             tttAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgc                                                                                             <  1:12931/141‑2 (MQ=255)                                              ttAAGTTGTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgca                                                                                            <  2:40874/141‑3 (MQ=255)                                               tAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACgcaa                                                                                           <  1:161224/141‑4 (MQ=255)                                                aaGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCaaa                                                                                          <  1:476830/141‑1 (MQ=255)                                                  gTTTTTTGAGCAATATATTAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAg                                                                                        <  2:407500/141‑1 (MQ=255)                                                   ttttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGa                                                                                       <  1:45380/141‑1 (MQ=255)                                                    tttttGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGat                                                                                      >  1:388307/1‑141 (MQ=255)                                                        tgggcAAGATAGCAATTTTAAGGATTTAATGAATTAGAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGa                                                                                  <  2:468380/138‑1 (MQ=255)                                                        tGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGa                                                                                  >  2:273222/1‑141 (MQ=255)                                                           gcAATATAGCAATTTTAAATATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAgatatgatat                                                                               <  1:132947/141‑1 (MQ=14)                                                           gcAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCGGCTAAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGCAAAAGATATGatat                                                                               >  1:213674/1‑138 (MQ=14)                                                                  agtaaTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaaca                                                                        <  2:375290/138‑1 (MQ=255)                                                                   gCAATTTTAAATATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGGTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACcaacaa                                                                       <  1:416692/141‑1 (MQ=255)                                                                      aTTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATa                                                                    >  1:489610/1‑141 (MQ=255)                                                                        tttAAAGATTTAATGAATTATAAAGATGGAGATATATCATAAAATACAGACATACCAATCTATTCAACTCAATATCAATTAACAAATGAACATTATAATATAAAACACTCAACTACCACATATAACATACCAACCAATaaa                                                                  >  2:480414/1‑141 (MQ=255)                                                                         ttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATaaaa                                                                 >  1:254202/1‑140 (MQ=255)                                                                         ttAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGATCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTAAGTAAATTATATTAAATACCAAAAAATaaaa                                                                 >  2:324862/1‑140 (MQ=18)                                                                          tAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAACATCGGCGTGTAAGTTAGGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAgc                                                                >  1:235415/1‑141 (MQ=18)                                                                                    aTGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtg                                                      <  1:59007/141‑1 (MQ=255)                                                                                       aagtataAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCtgttg                                                   <  1:479993/138‑1 (MQ=255)                                                                                                  tGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGATGTAATTGAAAGGTAc                                        >  2:53361/1‑141 (MQ=255)                                                                                                    ggaatatatCATATAATGCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAg                                      <  2:284489/138‑1 (MQ=255)                                                                                                     cgaTATAGCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACAAACGATGATTATAATGTAAAACAATTACGTAAAAGACAGGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAgg                                     <  1:416102/140‑1 (MQ=17)                                                                                                      aatatatGATATAATACAGAGTTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACGAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAggg                                    <  2:290388/140‑1 (MQ=255)                                                                                                        tatatCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGaa                                  >  1:77595/1‑141 (MQ=255)                                                                                                               tataATCCAGAGGTGACGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATATAGCGCCGAAGCTGGTGTTGAACGGTACAGGGAATTTaaaa                           >  1:16710/1‑141 (MQ=255)                                                                                                                   aTCCAGAGGTGCCGAGTTATTAAGCTCAATAACAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATAACAACAAATAAAGCGCGGAAGCTGTTGTTGAAAGGTACAGGAAATTAAAAAGGtt                       >  2:97898/1‑141 (MQ=14)                                                                                                                       agagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatc                   <  2:49866/141‑1 (MQ=255)                                                                                                                        gagGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGATGTTGTTGAAAGGTACAGGGAATTTAAAAGGTtcatca                  >  2:97648/1‑141 (MQ=255)                                                                                                                         agGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGAGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAACAGGTTCATCAg                 >  2:341514/1‑141 (MQ=255)                                                                                                                            tGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACTAAACAAAAACGAGCCTACGCTGTTGTTGAAAGATACAGGGAATTACAAAGGTCCATCAGTTg              >  2:105667/1‑141 (MQ=17)                                                                                                                             gCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAAAAAAGACTCGTAGCAGATGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTgg             >  2:299941/1‑141 (MQ=14)                                                                                                                                gAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAACGCTTCATCAGTTGGata          >  2:246732/1‑141 (MQ=255)                                                                                                                                 aGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGatat         >  2:445560/1‑141 (MQ=255)                                                                                                                                        cAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATaaaaaaa  >  2:10869/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                        CTGAGGAAGTTAAAGAAGAAAACATAAAAAAAGAAATCGAAAACTTTAAGTTTTTTGCGCAATATAGCAATTTTAAAGATTTAATGAATTATAAAGATGGAGATATATCATATAATCCAGAGGTGCCGAGTTATTCAGCTCAATATCAATTAACTAATGATGATTATAATGTAAAACAATTACGTAAAAGATATGACATACCAACAAATAAAGCGCCGAAGCTGTTGTTGAAAGGTACAGGGAATTTAAAAGGTTCATCAGTTGGATATAAAAAAA  >  USA300TCH1516_ALE/111376‑111651

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