Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,532,200	G→A	6.9%	R223C (CGT→TGT)	USA300HOU_RS07470 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,532,200	0	G	A	6.9%	176.0 / 3.6	62	R223C (CGT→TGT)	USA300HOU_RS07470	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/18);  new base A (2/2);  total (42/20)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

GAATAATCTACCATTAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGTTTCATCGCTCTTGAACTCGCAT  >  USA300TCH1516_ALE/1532060‑1532340                                                                                                                                             |                                                                                                                                             gAATAATCTACCATTAGAAAAATAAGCACAAAAAAACTACCAACCATAGAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACg                                                                                                                                              <  2:80153/141‑1 (MQ=255) aaaaaatCTAACAATAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACg                                                                                                                                              <  2:79171/137‑1 (MQ=255)        cTACCATTAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGc                                                                                                                                       >  1:89550/1‑141 (MQ=255)        cTACCATTAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGc                                                                                                                                       >  1:296027/1‑141 (MQ=255)          aCCATTAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGctc                                                                                                                                     <  2:552069/141‑1 (MQ=255)                    aaaTAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCCAAATAATACGTACAACCTCTGCTACAGCg                                                                                                                           >  2:180332/1‑141 (MQ=255)                    aaaTAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCg                                                                                                                           <  1:319203/141‑1 (MQ=255)                         aGCACAAAAAAACTAACAACCACACAAGAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAATATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGc                                                                                                                      <  2:103573/141‑1 (MQ=255)                          gCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTAAAGCGACTGCa                                                                                                                     >  1:704005/1‑141 (MQ=255)                           cacaAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAg                                                                                                                    >  1:219424/1‑141 (MQ=255)                               aaaaaaCTAGCAACCCAACAAAAATTGCATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTtgc                                                                                                                <  1:207613/141‑1 (MQ=255)                               aaaaaaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTtgc                                                                                                                <  2:158452/141‑1 (MQ=255)                                 aaaaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTtgctg                                                                                                              >  1:355408/1‑141 (MQ=255)                                  aaaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACATACAAGCTCTGCTACAGCGACTGCAGTtgctgc                                                                                                             >  1:697251/1‑141 (MQ=255)                                  aaaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACATACAAGCTCTGCTACAGCGACTGCAGTtgctgc                                                                                                             >  1:308842/1‑141 (MQ=255)                                  aaaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTAAAAGAGCCACTAAAATTATACGTACAAGCTCTGCTACAGCGACTGCAGTtgctac                                                                                                             >  2:349929/1‑139 (MQ=255)                                  aaaCTAGCAACCACACAAAAATGTGATTAGCGAGTTAATAAGTGTCTAATGTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAAAGCGAAATCACAACACGTCCTTCAAGCTCTGCTACAGCGACTGCAGTtgctac                                                                                                             >  2:741192/1‑139 (MQ=255)                                   aaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTCATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCa                                                                                                            >  2:456499/1‑141 (MQ=255)                                   aaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCa                                                                                                            >  1:194838/1‑141 (MQ=255)                                   aaCTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCa                                                                                                            >  2:634917/1‑141 (MQ=255)                                     gtAGCAACCACACGAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATGTATAAGATTAATCACTTGAACGGGCAATCAAAATAATACGTGCAAGCTCTGCTACTGCGACTGCAGTTGCTGCAAc                                                                                                          >  2:607412/2‑141 (MQ=255)                                     cTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTATTCTATAAGGTTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAAc                                                                                                          >  2:607421/1‑141 (MQ=255)                                     cTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTCCTACAGCGACTGCAGTTGCTGCAAc                                                                                                          >  2:100217/1‑141 (MQ=255)                                       aGCAACCACACAAACATGTGATTAGCTAGTTAATAAGTGTCTACTTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAACTAATACGTACAAGCTCTGCTACAACAACTGAAGTTACTGCAACaa                                                                                                        >  2:124931/1‑140 (MQ=255)                                       aGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCACAATCAAACTAATACGTACAAGCACTGCTACAGCGACTCACATTGCTGCAACAt                                                                                                        >  2:407441/1‑141 (MQ=255)                                          aaCCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAg                                                                                                     <  1:14704/141‑1 (MQ=255)                                             cacacaAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAAAGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGGGACTGCAGTTGCTGCAACATAAGTct                                                                                                  >  2:530918/1‑140 (MQ=255)                                                acaAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtg                                                                                               >  1:390743/1‑141 (MQ=255)                                                acaAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACATACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtg                                                                                               <  2:697251/141‑1 (MQ=255)                                                acaAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACATACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtg                                                                                               <  2:308842/141‑1 (MQ=255)                                                 caAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCGAGAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtga                                                                                              >  2:726148/1‑140 (MQ=255)                                                 caAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGGTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtgc                                                                                              <  2:467714/141‑1 (MQ=255)                                                    aaaTGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTTCTACAGCGACTGCAGTTGCTGCAACATAAGTCATtgctgc                                                                                           <  1:553122/141‑1 (MQ=255)                                                         tgATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAg                                                                                      >  1:57711/1‑141 (MQ=255)                                                              aGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGataat                                                                                 <  1:680209/141‑1 (MQ=255)                                                                  aGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACtt                                                                             <  2:604941/141‑1 (MQ=255)                                                                   gTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACttt                                                                            >  1:317053/1‑141 (MQ=255)                                                                        aaaGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTAcgcg                                                                       <  2:152904/140‑1 (MQ=255)                                                                           gtgtCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCGCCAGATAATACTTTCCGCGCAt                                                                    >  2:721064/1‑141 (MQ=255)                                                                             gtCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGc                                                                  >  2:740025/1‑141 (MQ=255)                                                                                 aaTTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGt                                                              <  1:663949/141‑1 (MQ=255)                                                                                  aTTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTTCAAGCTCTGCTACAGCGACTGCAGTTGCTGCACCATAAGTCATTGCTGCTGCAGATAATACTTTCAGCGCATGCTTCTa                                                             >  2:155442/1‑141 (MQ=255)                                                                                    ttAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCAGGCTTGTAtt                                                           >  2:41289/1‑141 (MQ=255)                                                                                       aGTTAATTGTTAATCTATAAGATTAATCACGTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCtt                                                        <  2:489550/141‑1 (MQ=255)                                                                                        gTTAATTGATAATCTATCAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCttt                                                       <  2:710954/141‑1 (MQ=255)                                                                                         ttAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATCATACGTACAAGCTGTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTACTGCTGCAGATGATACTTTACGCGCATACTTGTATTCtttt                                                      >  2:18598/1‑141 (MQ=255)                                                                                            aTTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTATATTCTTTTTCa                                                   >  2:602907/1‑141 (MQ=255)                                                                                              tGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCAtt                                                 <  1:634917/141‑1 (MQ=255)                                                                                               gTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCAttt                                                >  1:628740/1‑141 (MQ=255)                                                                                                     cTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATg                                          <  2:718960/141‑1 (MQ=255)                                                                                                         aaGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCa                                      >  1:270948/1‑141 (MQ=255)                                                                                                               aaTCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCg                                >  1:394585/1‑141 (MQ=255)                                                                                                               aaTCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCg                                >  1:389275/1‑141 (MQ=255)                                                                                                               aaTCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCg                                >  2:608159/1‑141 (MQ=255)                                                                                                                aTCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGt                               >  2:700570/1‑141 (MQ=255)                                                                                                                 tCACTTGAAACCGCACTCAAAATAATACGTACAAGCTCTGCCACAGCGACTGCAGTTGCTGCAACATAAATCATTGCTGCTGCAGATAATACTTTGCGCGCATGCTTGTATTCTTTTCCATTTACAATGTTCAATGCCGTa                              >  2:750853/1‑141 (MQ=255)                                                                                                                    cTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTGGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAAta                           >  2:472583/1‑140 (MQ=255)                                                                                                                       taaCGCGCAATCAAAATAATACGTACAAGCTCCGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCAGTTACAATGTTCAATGCCGTAATTTGt                        <  1:183288/140‑1 (MQ=255)                                                                                                                         aCGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGttt                      >  1:578679/1‑141 (MQ=255)                                                                                                                         aCGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGttt                      >  1:233918/1‑141 (MQ=255)                                                                                                                          cgcgCAATCAAAATAATCCGTACAAGCTCGGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCACGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGTTTc                     <  1:126514/141‑1 (MQ=255)                                                                                                                                aTCAAAATAATACGTAAAAGCGCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCAGGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTACATTGtttgctcccc               >  2:41599/1‑134 (MQ=255)                                                                                                                                 tCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGTTTCATCGctc              >  2:67921/1‑141 (MQ=255)                                                                                                                                         ataCGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGGTTCATCGCTCTTGAACTc      >  2:146981/1‑141 (MQ=255)                                                                                                                                             gTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGTTTCATCCATCTTGAACTCGCAt  >  2:95140/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                             GAATAATCTACCATTAGAAAAATAAGCACAAAAAAACTAGCAACCACACAAAAATGTGATTAGCTAGTTAATAAGTGTCTAATTTAAGTTAATTGTTAATCTATAAGATTAATCACTTGAACGCGCAATCAAAATAATACGTACAAGCTCTGCTACAGCGACTGCAGTTGCTGCAACATAAGTCATTGCTGCTGCAGATAATACTTTACGCGCATGCTTGTATTCTTTTTCATTTACAATGTTCAATGCCGTAATTTGTTTCATCGCTCTTGAACTCGCAT  >  USA300TCH1516_ALE/1532060‑1532340

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