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	98.6%	233.9 / NA	69	M173R (ATG→AGG)	mecA	Adapter protein MecA 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (25/43);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGAAGCTCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTAAAA  >  USA300TCH1516_ALE/1002054‑1002340                                                                                                                                                   |                                                                                                                                             agaagCTCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAGTAACTACAGAGATTGAAGATTTGTTATATAggg                                                                                                                                            >  1:80616/1‑149 (MQ=255) agaagCTCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAggg                                                                                                                                            >  1:145599/1‑149 (MQ=255)     gCTAAAGGTTTTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACATTCATTTGGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGa                                                                                                                                        <  1:192225/149‑1 (MQ=255)       tCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGCTTTGTTATATAGGGTTGATg                                                                                                                                      >  1:251048/1‑149 (MQ=255)       tCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATg                                                                                                                                      >  1:898303/1‑149 (MQ=255)          aGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAAGACAAGCAAATACGCATTTAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTa                                                                                                                                   <  1:913327/149‑1 (MQ=255)             ttCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACtt                                                                                                                                <  1:204688/149‑1 (MQ=255)             ttCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACtt                                                                                                                                <  1:411225/149‑1 (MQ=255)              tCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACtta                                                                                                                               <  1:11777/149‑1 (MQ=255)                tAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatt                                                                                                                             <  1:133104/149‑1 (MQ=255)                   aCGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttattatt                                                                                                                          >  1:749861/1‑149 (MQ=255)                         aaaGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATACCTACAGGGTTAGAAGTTTTGTTATATAGGGTTGATGGTACTTATTATTATGCTg                                                                                                                    >  1:23131/1‑149 (MQ=255)                          aaGTCTTCAACAAGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGt                                                                                                                   <  2:320952/149‑1 (MQ=255)                           aGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGta                                                                                                                  >  2:8238/1‑149 (MQ=255)                             tCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGtata                                                                                                                >  1:360774/1‑149 (MQ=255)                               ttCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGCTGCTGGTACTTATTATTATGCTGTATAtt                                                                                                              <  1:731729/149‑1 (MQ=255)                                  aGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGGTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTg                                                                                                           >  2:214878/1‑149 (MQ=255)                                     aCGTAAAAATACAAGAACAATCATTGTGATATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGGTGGTACTTATTATTATGCTGTATATTTTGATa                                                                                                        >  1:481325/1‑149 (MQ=255)                                     aCGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATa                                                                                                        >  1:712525/1‑149 (MQ=255)                                     aCGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTAGGCATATCATAGCAAGCCAATAACTACAGAGTTTGAAGATTTGTGATATAGGGCTGATGGTACGTATTATTATGCGGTATATTTGGATa                                                                                                        >  1:420321/1‑149 (MQ=255)                                      cGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAg                                                                                                       <  1:356184/149‑1 (MQ=255)                                      cGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAg                                                                                                       <  2:56505/149‑1 (MQ=255)                                          aaaaTACAAGAACAATCATTGTGAAATTTAACGCTTTAGAAGATGTTATTAATTATGCATATAATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTAATTATTATTATGCTATATTTTTTtatatgctt                                                                                                   >  2:543566/1‑144 (MQ=255)                                                cAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGAt                                                                                             >  2:472464/1‑149 (MQ=255)                                                    aaCAATCATTGTGAAATTTAACGATTTAGAAGAACTTTTTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaag                                                                                         >  1:449500/1‑149 (MQ=255)                                                       aaTCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaagaag                                                                                      >  1:560008/1‑149 (MQ=255)                                                        ttCATTGTGAAATTTATCGGTTTAGAGGTTGTTATTTATTTTGCTTTTTATTGCAATCCTATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGt                                                                                     <  1:29822/148‑1 (MQ=255)                                                          taTTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCa                                                                                   <  2:626107/148‑1 (MQ=255)                                                          cATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCa                                                                                   <  2:156117/149‑1 (MQ=255)                                                           aTTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAAGTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGTTGTATATGTTGATAGTCATGTTGATCAAGAAGTCAt                                                                                  <  1:602604/149‑1 (MQ=255)                                                            ttGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCAtt                                                                                 <  1:73031/149‑1 (MQ=255)                                                            ttGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCAtt                                                                                 <  1:47615/149‑1 (MQ=255)                                                                 aaaTTTAGCGATCTATGAGATGTTATTAATTATGCATATCCTAGCAAGCCTATAATTATAGAGTTTGAATGTTTTTTATATGGGGTTGATGGTATTTTTTATTTTGCTGGATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGa                                                                            <  2:328305/149‑1 (MQ=255)                                                                 aaaTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGa                                                                            <  2:493535/149‑1 (MQ=255)                                                                   aTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGACGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATa                                                                          >  1:586414/1‑149 (MQ=255)                                                                   aTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATa                                                                          >  1:905921/1‑149 (MQ=255)                                                                       aaCGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTa                                                                      >  1:589433/1‑149 (MQ=255)                                                                        aCGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATATTCATGTTGATCAAGAAGTCATTAATGATAGTTAc                                                                     <  2:578974/149‑1 (MQ=255)                                                                        aCGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATTAACAAGGAATTAAAGATAATTAc                                                                     >  2:26528/1‑149 (MQ=255)                                                                        aCGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTAc                                                                     <  1:748776/149‑1 (MQ=255)                                                                          gATTGAGAAGATGTTAGTAATTATGCATATCATAGCAATCCAATAACCACAGATTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAg                                                                   <  2:272405/149‑1 (MQ=255)                                                                           aTTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAgg                                                                  >  1:577971/1‑148 (MQ=255)                                                                               agaagaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTTATTAATGATAGTTACAGTCAAt                                                              <  1:498421/149‑1 (MQ=255)                                                                                  agaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGc                                                           <  2:399073/149‑1 (MQ=255)                                                                                     tGTTATTAATTATTCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTg                                                        <  1:311081/149‑1 (MQ=255)                                                                                     tGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGTTTGATAGAACTTTGTTTCATGCGGTATATTTTGATAGTTATGTTGATCAATAAGTCATTAATGATAGTTACAGTCATTTGCTTg                                                        >  1:681810/1‑149 (MQ=255)                                                                                     tGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTg                                                        <  1:569394/149‑1 (MQ=255)                                                                                             attaTGCATATCATAGCAAGCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCt                                                <  1:247309/149‑1 (MQ=255)                                                                                              ttaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTCCTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCtt                                               >  1:258565/1‑149 (MQ=255)                                                                                                     tatCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAAGTGCTTGAATTTGCTTATCCAac                                        <  1:498042/149‑1 (MQ=255)                                                                                                          tAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagaca                                   <  1:919921/149‑1 (MQ=255)                                                                                                          tAGCAATCCAATAACTACAGAGTTTGAAGATGTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCAAacagaca                                   >  1:575646/1‑149 (MQ=255)                                                                                                           aGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagacag                                  >  2:839962/1‑149 (MQ=255)                                                                                                              aaTCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaac                               >  1:72026/1‑149 (MQ=255)                                                                                                                tCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaacag                             >  2:62366/1‑149 (MQ=255)                                                                                                                      aaCTACATAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGttt                       <  2:529129/149‑1 (MQ=255)                                                                                                                        cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttat                     <  1:544615/149‑1 (MQ=255)                                                                                                                        cTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttat                     <  1:839962/149‑1 (MQ=255)                                                                                                                         tACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttatt                    <  2:145599/149‑1 (MQ=255)                                                                                                                             gagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaaa                <  1:744030/149‑1 (MQ=255)                                                                                                                             gagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaaa                >  2:721408/1‑149 (MQ=255)                                                                                                                             gagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaaa                <  2:835655/149‑1 (MQ=255)                                                                                                                                 ttGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGAc            <  2:523679/149‑1 (MQ=255)                                                                                                                                     aGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATg        <  2:360774/149‑1 (MQ=255)                                                                                                                                     aGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATg        <  1:873295/149‑1 (MQ=255)                                                                                                                                     aGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATg        <  1:873171/149‑1 (MQ=255)                                                                                                                                     aGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAAGTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATg        <  2:852701/149‑1 (MQ=255)                                                                                                                                      gttttGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGc       <  2:227425/147‑1 (MQ=255)                                                                                                                                       aTTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCt      <  2:865312/149‑1 (MQ=255)                                                                                                                                         ttGTCACATAGGGTTGATGGCACTTATTAGTGTGGTGTGACTTTCGTTAGTCAGGTTGTTTAAGAAGTCATTAATGATCGTTTCATTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTaa    <  2:195005/149‑1 (MQ=255)                                                                                                                                           tttATATAGGGTTGATGTTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTaaaa  <  2:238576/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AGAAGCTCAAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTAAAA  >  USA300TCH1516_ALE/1002054‑1002340

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