Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	1,002,200	T→G	M173R (ATG→AGG)	mecA →	Adapter protein MecA 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,002,200	0	T	G	97.3%	237.3 / ‑6.0	74	M173R (ATG→AGG)	mecA	Adapter protein MecA 1
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base G (28/44);  total (28/46)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTA  >  USA300TCH1516_ALE/1002062‑1002332                                                                                                                                           |                                                                                                                                     aaGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAATTTGAAGATTTTTTAAATAggg                                                                                                                                    >  2:262422/1‑141 (MQ=255)  aGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGt                                                                                                                                   <  1:944622/141‑1 (MQ=255)  aGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGt                                                                                                                                   <  1:740628/141‑1 (MQ=255)   ggtttaaaaCGACAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGtt                                                                                                                                  <  2:801320/135‑1 (MQ=255)   ggTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGtt                                                                                                                                  <  1:97924/141‑1 (MQ=255)   ggTGCTAAACGTCCAATGTCTTCAGCATGTAAAAATACAAGAACAATCATTGTGAAATGTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGTTTTGTTATATAGGGtt                                                                                                                                  <  1:586399/141‑1 (MQ=255)    gtgcTAAAGGGCAAAGGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTg                                                                                                                                 <  2:862965/138‑1 (MQ=255)      tCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGAt                                                                                                                               <  1:580400/141‑1 (MQ=255)        tAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACGACAGAGTTTGAAGATTTGTTATATAGGGTTGATgg                                                                                                                             <  2:825054/141‑1 (MQ=255)        aaaaCGTCAAAAGACTTCATCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATgg                                                                                                                             <  2:620973/140‑1 (MQ=255)           aCGTCAAAAGTCTTCAGCACGGAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTAc                                                                                                                          <  1:204362/141‑1 (MQ=255)                 aaaGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatta                                                                                                                    <  2:495922/141‑1 (MQ=255)                 aaaGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatta                                                                                                                    <  2:295070/141‑1 (MQ=255)                       tgcAGCACGTAAAAATACAAGAACAAGCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGCTATATAGGGTTGATGGTACTTCTTCTTATGc                                                                                                              <  1:726074/139‑1 (MQ=255)                           gCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGta                                                                                                          >  1:467760/1‑141 (MQ=255)                           gCAAGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGta                                                                                                          <  1:214986/141‑1 (MQ=255)                               gTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATAtt                                                                                                      <  2:99090/141‑1 (MQ=255)                                  aaaaCACAAGAATAATCATTGAGAAATTTAACGATGTAGAAGATGTTATTAATTACGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGACGGTACTTATTATTATGCTGTATATTTTg                                                                                                   <  1:748779/141‑1 (MQ=255)                                   aaaTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGa                                                                                                  <  2:666074/141‑1 (MQ=255)                                              cAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTg                                                                                       >  1:939819/1‑141 (MQ=255)                                               aagcATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTgt                                                                                      >  2:411397/4‑140 (MQ=255)                                                  cATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTCTCTCGGGTAGATGGTACTAATTATTATGCTGTATATTTTGATGGACATGTTGATCa                                                                                   >  1:871053/1‑141 (MQ=255)                                                    ttGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaag                                                                                 >  1:846778/1‑141 (MQ=255)                                                     tgtgAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaaga                                                                                <  2:118606/141‑1 (MQ=255)                                                     tgtgAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaaga                                                                                >  1:877548/1‑141 (MQ=255)                                                          aaTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGACTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCa                                                                           <  2:269245/141‑1 (MQ=255)                                                           aTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGAACTCATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCAt                                                                          >  1:321442/1‑141 (MQ=255)                                                             ttAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTa                                                                        >  2:658147/1‑141 (MQ=255)                                                             ttAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTa                                                                        >  2:77021/1‑141 (MQ=255)                                                              tAACGATTTAGAAGCTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTaa                                                                       <  1:252279/141‑1 (MQ=255)                                                               accGATTTCGAAGATGTTATTCATTCTGAATATCCTAGCAATCCAATCACCACCGAGCTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAAt                                                                      <  1:447235/139‑1 (MQ=255)                                                               aaCGATTTAGAAAATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAAt                                                                      <  2:180774/141‑1 (MQ=255)                                                                aCGATTTAGAAGATGTTAGTAATTATGCATATCATAGCAAGTCATTAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATg                                                                     <  2:53578/141‑1 (MQ=255)                                                                   aTTTAGAAGATGTTATTAATTATGCATAGCATATCAATCCAATAACTACAGAGCTAGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGTTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATa                                                                  <  2:544124/141‑1 (MQ=255)                                                                   aTTTAGAAGACGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTGGTTATATAGGGTTGATGGTACTTATTAGTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATa                                                                  <  2:365236/141‑1 (MQ=255)                                                                         aagaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACa                                                            <  2:67612/141‑1 (MQ=255)                                                                               ttattaATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAAt                                                      >  2:854736/1‑141 (MQ=255)                                                                                 attaattaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTg                                                    >  2:41027/1‑141 (MQ=255)                                                                                   taattaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCt                                                  <  1:87862/141‑1 (MQ=255)                                                                                   taattaCTCATATCCAAGCAATTCTATAACTTCATAGTTTGAAGATTAATTAAAAAGCGTTGATGGTTCATGTTCGTCAGAGGTTGTATTGGCTTGACGTGTTGATCAAGGAATCATTAATGATAGTTACAGTCCATTGCt                                                  <  2:388763/141‑1 (MQ=255)                                                                                    aattaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCtt                                                 >  1:71264/1‑141 (MQ=255)                                                                                        aTGCATATCATAGCAATCCAATAACTACAGGGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAAt                                             <  2:471978/141‑1 (MQ=255)                                                                                         tGCATATCATAGCAATCCAATAACTAGAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAAtt                                            <  2:39124/141‑1 (MQ=255)                                                                                           cATATCAAAGCAATCCAATAACTACAGATTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGCCATTAATGATAGTTACAGTCAATTGCTTGAATTTg                                          <  2:876674/141‑1 (MQ=255)                                                                                               tCATAGCAATGCAATAACTACAGATTTTGAAGATTTGTTATATAGGGTTGATTGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATCATAGTTACAGTCAATTGCTTGAATTTGCTTa                                      <  2:386354/141‑1 (MQ=255)                                                                                               tCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATGAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTa                                      >  1:546368/1‑141 (MQ=255)                                                                                               acATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGACAAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTa                                      <  1:434406/140‑1 (MQ=255)                                                                                                 aTAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATc                                    >  2:296808/1‑141 (MQ=255)                                                                                                 aTAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATc                                    >  1:350613/1‑141 (MQ=255)                                                                                                 aTAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATc                                    >  1:894328/1‑141 (MQ=255)                                                                                                 aTAGCAATCCAATAACTAAGGAGAAGGAAGCTTAGTTATATAGGGTTGATGGTACTTATTATTATGCTGTAGATTTTGATCGTCATGTTGATCAAGAAGCCATTAATGTTAGTTACAGTAAATTGCTTGAATTTGCtttca                                    >  2:271812/1‑138 (MQ=255)                                                                                                 aTAGCAATCCAATAAATACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCGATTGCTTGAATTTGCTTATc                                    >  1:731837/1‑141 (MQ=255)                                                                                                         ccAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATACAccagac                            >  2:337767/1‑141 (MQ=255)                                                                                                           aaTAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagacag                          >  1:847716/1‑141 (MQ=255)                                                                                                            atcaCTACAGAGTTTCAAGCTTTGTTATATATGGTTGATGGTACTTATCATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAAcagacaga                         <  1:62625/138‑1 (MQ=255)                                                                                                             aaaCTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaa                        <  2:35752/140‑1 (MQ=255)                                                                                                               aCTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaaca                      <  2:939819/141‑1 (MQ=255)                                                                                                                cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     <  1:439871/141‑1 (MQ=255)                                                                                                                cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     >  1:833187/1‑141 (MQ=255)                                                                                                                cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     <  1:259544/141‑1 (MQ=255)                                                                                                                cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     <  1:796328/141‑1 (MQ=255)                                                                                                                cTACAGAGTTGGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     <  1:252987/141‑1 (MQ=255)                                                                                                                cTACAGAGGTTGAAGATTTGTTAGATAGGGTTGATGGTACTTATTATTACGCTTGATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                     <  1:131517/141‑1 (MQ=255)                                                                                                                    agagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGt                 >  1:747769/1‑141 (MQ=255)                                                                                                                       gTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGATTGAATTTGCTTATCCAACAGACAGAACAGAAGttta              <  2:871053/141‑1 (MQ=255)                                                                                                                        ttttAAGAGTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTACATTTTGATAGTCATGTTGATCAAGAGGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttat             <  2:99716/141‑1 (MQ=255)                                                                                                                         ttGAAGATTTTTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttatt            <  1:957152/141‑1 (MQ=255)                                                                                                                          tGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttattt           >  1:639720/1‑141 (MQ=255)                                                                                                                            aaGAATTGTTAAATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaa         <  1:797367/141‑1 (MQ=255)                                                                                                                             aGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaaa        >  1:372459/1‑141 (MQ=255)                                                                                                                              gATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAGTTTGCTTATCCAACAGACAGCCAAGATGGTTATTTAAAt       >  2:152790/1‑141 (MQ=255)                                                                                                                                 ttGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGAc    <  2:583363/141‑1 (MQ=255)                                                                                                                                   gTTATATAGGGTTGATGGTAGTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGGCATTAATGATAGTTACAGTGAATTGCTTGAATTTGCTTATCCAACAGACAGGACAGAAGTTTATTTAAATGACTa  >  1:562690/1‑141 (MQ=255)                                                                                                                                   gTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAGGTCATTAATGATAGTTACAGTCAATTGCTCGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACtt  >  2:838538/1‑140 (MQ=255)                                                                                                                                   gTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTa  >  2:254334/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                     AAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTA  >  USA300TCH1516_ALE/1002062‑1002332

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: ‑