Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,400,393	C→T	6.1%	P73L (CCA→CTA)	dinB1 →	DNA polymerase IV 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,400,393	0	C	T	6.1%	188.7 / 3.7	66	P73L (CCA→CTA)	dinB1	DNA polymerase IV 1
Reads supporting (aligned to +/- strand):  ref base C (35/26);  new base T (2/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TTGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAT  >  USA300TCH1516_ALE/1400255‑1400532                                                                                                                                           |                                                                                                                                            ttGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACcac                                                                                                                                           >  1:199186/1‑141 (MQ=255)    aaaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACATAAATCACAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATa                                                                                                                                        >  2:617167/1‑141 (MQ=255)    aaaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGAGTCAGGGTTCTGTAATATTGGCTGCGCCACCTAAATTAAAAGAATTAGGCATCAAGTCCGGGTCGCGATTGTTTGAAATACCACATa                                                                                                                                        <  2:435815/141‑1 (MQ=255)     aaaGGGGCTAGATGCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAg                                                                                                                                       >  2:340682/1‑141 (MQ=255)     aaaGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAg                                                                                                                                       >  1:3467/1‑141 (MQ=255)         gggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  1:496080/1‑141 (MQ=255)         gggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  1:632195/1‑141 (MQ=255)         gggCTAAATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAATAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAAt                                                                                                                                   >  2:592349/1‑141 (MQ=255)          ggCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATg                                                                                                                                  <  2:57118/141‑1 (MQ=255)            cTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGat                                                                                                                                <  1:945079/141‑1 (MQ=255)            cTAGATCCATTAGAAACAAAGCCAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGACTTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGat                                                                                                                                <  1:265882/141‑1 (MQ=255)             tAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACAAAGAAATGata                                                                                                                               >  2:741916/1‑141 (MQ=255)                 tCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTa                                                                                                                           <  2:760803/141‑1 (MQ=255)                  ccATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTAc                                                                                                                          <  1:9043/141‑1 (MQ=255)                      tAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGTGATTGTTTGAAATACCACATAGAAATGATATTTACATTa                                                                                                                      <  1:346189/141‑1 (MQ=255)                        gAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATc                                                                                                                    >  1:700839/1‑141 (MQ=255)                        gAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGAAACTGTAATATTGGCTGCGACACCTAAATTAAAAGAAGTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATc                                                                                                                    >  2:417530/1‑141 (MQ=255)                        gAAACAAAGCTAGCTGTTGTTGACAAGACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATc                                                                                                                    >  1:457110/1‑141 (MQ=255)                             aaaGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATcc                                                                                                               >  2:836642/1‑141 (MQ=255)                              aaGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGGTTGTTTGAAATACCACATAGAAATGATGTTTACATGATAAATCCa                                                                                                              >  2:15671/1‑141 (MQ=255)                               aacttgctGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGAGATTGTTTGAAATACAACATAGAAATGATATTAACATTATCAATCCaa                                                                                                             >  2:112480/6‑141 (MQ=255)                                   agctGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTAt                                                                                                         >  2:127405/1‑141 (MQ=255)                                       gttgttGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGt                                                                                                     >  1:776514/1‑141 (MQ=255)                                            tgCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAata                                                                                                <  2:690857/141‑1 (MQ=255)                                               aGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCt                                                                                             <  2:640468/141‑1 (MQ=255)                                                   aCTAATCGTCAGGGTTCTGTAATATTGGCTGTGACACCTAGATTAAAAGAATTAGGCATCAAGACAAGGTCGCGATTGTTTGAAGTACCACAGAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAAt                                                                                         >  1:694399/1‑141 (MQ=255)                                                   aCTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATAACTTAAt                                                                                         >  2:88846/1‑141 (MQ=255)                                                    cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGTGTCGCTATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        <  2:315881/141‑1 (MQ=255)                                                    cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        <  2:363800/141‑1 (MQ=255)                                                    cTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATg                                                                                        >  2:922680/1‑141 (MQ=255)                                                       aGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGttt                                                                                     >  1:69587/1‑141 (MQ=255)                                                        gCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTc                                                                                    >  2:869572/1‑141 (MQ=255)                                                            cAGGGTTCTGTAATATTGGCTGCAACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGtt                                                                                >  2:674811/1‑141 (MQ=255)                                                               ggTTCTGTACAATTGGCTGCGACAACTAAATTAAACGCTGGATACATACAGTCAGTGCCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCt                                                                             >  1:690286/1‑141 (MQ=255)                                                                  tCTGTAATATTGGCTGCGAACGCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTAtt                                                                          >  1:226838/1‑141 (MQ=255)                                                                        agatggacTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAg                                                                    <  1:440556/134‑1 (MQ=255)                                                                          atgtgCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTATATTATCAACCCTAGTATGCGTAATTATCTTAATGTTTCAGTTGCTATTTCTAAGat                                                                  <  2:444141/137‑1 (MQ=255)                                                                              gCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTACCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgca                                                              >  1:297519/1‑141 (MQ=255)                                                                               cTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcat                                                             >  2:861732/1‑141 (MQ=255)                                                                               cTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcat                                                             >  1:533477/1‑141 (MQ=255)                                                                                tgggACACCGAAATTAAAAGAATTAGGGATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcatt                                                            <  1:382235/138‑1 (MQ=255)                                                                                tGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcatt                                                            <  1:245004/141‑1 (MQ=255)                                                                                  cGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcattgc                                                          <  2:700839/141‑1 (MQ=255)                                                                                  cGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACTACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGattgcattgc                                                          <  2:457110/141‑1 (MQ=255)                                                                                    acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCCTTGCGt                                                        >  2:461839/1‑141 (MQ=255)                                                                                    acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGt                                                        >  2:285826/1‑141 (MQ=255)                                                                                    acacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGt                                                        >  2:644428/1‑141 (MQ=255)                                                                                     cacCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGtt                                                       <  1:741916/141‑1 (MQ=255)                                                                                          aaaTTAAAAGAATTAGGCATCAAGACAGGTTCGCGATTGTTTGAAATACCACGTAGAAAGGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  2:382858/141‑1 (MQ=255)                                                                                          aaaTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTACATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  1:414003/141‑1 (MQ=255)                                                                                          aaaTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATAtt                                                  <  1:236526/141‑1 (MQ=255)                                                                                            aTTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTcc                                                <  1:127405/141‑1 (MQ=255)                                                                                                  aGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTAAGTTGCTATTTCTAAGATTGCATTGCGTTATTTTCAACCTGa                                          >  2:367562/1‑141 (MQ=255)                                                                                                   gAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGaa                                         >  1:743421/1‑141 (MQ=255)                                                                                                    aaTTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTACACCTGAAg                                        >  1:677805/1‑141 (MQ=255)                                                                                                        acacATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAATTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGAttt                                    <  1:1920/138‑1 (MQ=255)                                                                                                        aGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGAttt                                    >  1:39576/1‑141 (MQ=255)                                                                                                            aTCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTacac                                >  2:789464/1‑141 (MQ=255)                                                                                                             tCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACAcc                               <  2:406068/141‑1 (MQ=255)                                                                                                                  aCAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAatat                          >  1:153098/1‑141 (MQ=255)                                                                                                                   cAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGCAGATTTACACCAAtata                         >  2:630257/1‑141 (MQ=255)                                                                                                                           cgATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGAc                 <  1:303035/141‑1 (MQ=255)                                                                                                                            gATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACg                <  1:340682/141‑1 (MQ=255)                                                                                                                            gATTGTTTAAAATACCACATAGAAATGATATTTACATGATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACg                <  1:153428/141‑1 (MQ=255)                                                                                                                                 tttGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAAttt           >  1:108688/1‑141 (MQ=255)                                                                                                                                    gAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTGATATTCCACCTGAAGATTTACACCAATATGGTATTGAGGAAtttttt        >  1:21422/1‑141 (MQ=255)                                                                                                                                     aaaTACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTGCTAAGATTGCAGTGGGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTa       <  2:690286/141‑1 (MQ=255)                                                                                                                                          ccACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAt  >  2:285383/1‑141 (MQ=255)                                                                                                                                           |                                                                                                                                            TTGAAAAGGGGCTAGATCCATTAGAAACAAAGCTAGCTGTTGTTGCAGATACTAAGCGTCAGGGTTCTGTAATATTGGCTGCGACACCTAAATTAAAAGAATTAGGCATCAAGACAGGGTCGCGATTGTTTGAAATACCACATAGAAATGATATTTACATTATCAATCCAAGTATGCGTAAATATCTTAATGTTTCAGTTGCTATTTCTAAGATTGCATTGCGTTATATTCCACCTGAAGATTTACACCAATATAGTATTGACGAATTTTTTATGGAT  >  USA300TCH1516_ALE/1400255‑1400532

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