Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,024,507	C→A	10.1%	A37E (GCA→GAA)	USA300HOU_RS10210 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,024,507	0	C	A	10.1%	185.7 / 3.5	68	A37E (GCA→GAA)	USA300HOU_RS10210	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base A (2/6);  total (28/40)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATGCAAATTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGCGAAAGAAACA  >  USA300TCH1516_ALE/2024367‑2024646                                                                                                                                             |                                                                                                                                            aGATGCAAATTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGc                                                                                                                                             >  1:610196/1‑141 (MQ=255)     gCAAATTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGAaaaa                                                                                                                                         >  2:57074/1‑141 (MQ=255)       aaaTTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCaaaaaa                                                                                                                                       >  1:539052/1‑141 (MQ=255)       aaaTTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTGAGCTGATAGATCTAAACGTCAAGCTATAGATCAATCAGTTATAGATaaaaaaaa                                                                                                                                       >  2:445804/1‑141 (MQ=255)        aaTTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAAt                                                                                                                                      <  2:232375/141‑1 (MQ=255)                      gTAGTTTTTTATGGAAAATAAATTTGTGCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCACTCCGTTCAAGATGCAAAAAATGGTAACCGCACTCg                                                                                                                        <  1:404875/141‑1 (MQ=255)                       tAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTAAATCAGTTAAAGATGAAAAAAATGGTAACCGCACTAGt                                                                                                                       >  2:376844/1‑141 (MQ=255)                        tgTTTTTTATGGATAATAAAGTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTa                                                                                                                      <  2:222817/140‑1 (MQ=255)                        aGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATCGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGAAAAAAAAAAAAAAAATGGTAACCGCACTCGTa                                                                                                                      <  2:522074/141‑1 (MQ=255)                         gTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGGTAAAGATGCAAAAAATGGTCACCGCACTCGTaa                                                                                                                     >  1:430991/1‑141 (MQ=255)                          ttttttATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAg                                                                                                                    <  2:667263/141‑1 (MQ=255)                          ttttttATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAg                                                                                                                    <  2:187485/141‑1 (MQ=255)                          ttttttATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAg                                                                                                                    <  1:163068/141‑1 (MQ=255)                          ttatttATGGAAAACAAATTTGTTTCTGGTATTCTAATGGGTGCAGTAATGGGTGGTGCAATAAGTTTATCTGATAAATCCACACGTCAATCTTTAGTACACTCAGTCAAAGATGCAAAAAATGGTAACCGCACTCGTAAg                                                                                                                    <  2:336575/138‑1 (MQ=255)                           tttttATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGc                                                                                                                   >  1:340745/1‑141 (MQ=255)                            ttttATGGAAAAAAAATTTGTTCCGGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGcc                                                                                                                  <  2:99195/141‑1 (MQ=255)                               tATGGAAAATAAAGTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTc                                                                                                               <  2:168079/141‑1 (MQ=255)                                    aaaaTATATTTGTTCCTGGTATTTTACTTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTCTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAg                                                                                                          <  2:447421/141‑1 (MQ=255)                                     aaataaatTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGt                                                                                                         <  2:89578/141‑1 (MQ=255)                                          aatTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGGAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAgta                                                                                                    >  1:149935/1‑139 (MQ=255)                                          aatTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGgcagca                                                                                                    >  1:150406/1‑141 (MQ=255)                                          aatTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCa                                                                                                    >  1:680060/1‑141 (MQ=255)                                           atTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGAGGCAAAAAATGGTAACCGCACTCGTTAGCCTTCTAAAGTCAGCaa                                                                                                   <  1:423245/141‑1 (MQ=255)                                           atTGGTTCCTGCAATTTTAATCGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCaa                                                                                                   <  1:204277/141‑1 (MQ=255)                                              tGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGAt                                                                                                <  1:278228/141‑1 (MQ=255)                                                   cTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAg                                                                                           <  2:495751/141‑1 (MQ=255)                                                    aggTATTTTATTTGGGGCCGGAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCCTTAGTTCAGTAAGAAAAAAAAAAAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGa                                                                                          <  2:519187/140‑1 (MQ=255)                                                     ggTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTACGCCTTCTAAAGTCCGCAAGATTAAAGAc                                                                                         <  1:505602/141‑1 (MQ=255)                                                      gTATTTTAATTGGCGGCGTAATTGGCGGGGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAAAAAAAAAAAAAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACg                                                                                        <  2:693833/141‑1 (MQ=255)                                                             aaTTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGtttt                                                                                 <  2:750509/141‑1 (MQ=255)                                                                     ccGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATTAGTAAAAAAAAAAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGa                                                                         <  2:564958/141‑1 (MQ=255)                                                                      cGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGaa                                                                        <  1:57074/141‑1 (MQ=255)                                                                        tAATTGGGGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAg                                                                      <  2:208840/141‑1 (MQ=255)                                                                         aagtggtggTGCCATTAGGTCAGCTGATAAATCTACACGTTAAGCTTTCGTTCAATCAGGTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTTTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGa                                                                     <  2:53456/138‑1 (MQ=255)                                                                                 gtgCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGttgttg                                                             >  2:188980/1‑141 (MQ=255)                                                                                   gCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTgaa                                                           <  1:692955/141‑1 (MQ=255)                                                                                   gCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTGGTTgaa                                                           >  1:588999/1‑141 (MQ=255)                                                                                     aaTTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTgaaga                                                         <  2:241339/141‑1 (MQ=255)                                                                                        tAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAAt                                                      <  2:595026/141‑1 (MQ=255)                                                                                        tAGTTGACCTGATAAATCTAGACGTGGAGCTTCTGTTCAAAAAGTAAAAGATGAAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAAt                                                      <  2:740572/141‑1 (MQ=255)                                                                                            ttAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATtcgt                                                  >  1:32132/1‑141 (MQ=255)                                                                                               gCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTcgtcgt                                               <  2:233978/141‑1 (MQ=255)                                                                                                  gATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAGTGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGtaat                                            <  1:559728/141‑1 (MQ=255)                                                                                                  gATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGtaat                                            >  1:39047/1‑141 (MQ=255)                                                                                                    tAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTCACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGCCGAAGTTTTATACTGGACAGCTGTTGTTCAAGCAACTCCTCGtcatca                                          >  2:772953/1‑139 (MQ=255)                                                                                                    tAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGtaataa                                          >  2:772956/1‑141 (MQ=255)                                                                                                        tCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTACCCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACAAAGTTTTATACTGGAAAGATGTTGTTAAAGAATTTCGTAGTAATAATCCt                                      >  2:310167/1‑141 (MQ=255)                                                                                                        tCTACACGTCAAGCTTTAGTACAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAACCTTCTAAAGTCAGCAAGATCAAAAACAAACATTCATACTGGAAACCTATTGATAAAGAAATTAGACGCAATAATcat                                      >  2:213701/1‑139 (MQ=255)                                                                                                        tCTACACGCCAAGCTTTAGTTCGATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCt                                      >  2:564775/1‑141 (MQ=255)                                                                                                            cacCTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACCCCTAACCCTCCTGTTGGTAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAAt                                  <  1:367751/141‑1 (MQ=255)                                                                                                            caagTCAACCGTTAGTTCAATGAGTGAAACAGGGAAAATATGGGAACCGCCCGCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGGTGTTGAAGAAATTCGTCTTAATACTCCTGAAt                                  <  2:367550/138‑1 (MQ=255)                                                                                                                tCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAATTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGa                              <  1:23029/141‑1 (MQ=255)                                                                                                                tCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGa                              <  2:340745/141‑1 (MQ=255)                                                                                                                tCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGa                              <  1:667399/141‑1 (MQ=255)                                                                                                                  aaGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAAc                            >  1:93244/1‑141 (MQ=255)                                                                                                                  aaGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGACGTAATAATCCTGAATTAGAAc                            >  2:106269/1‑141 (MQ=255)                                                                                                                        tAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCAt                      >  1:206199/1‑141 (MQ=255)                                                                                                                        tAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCAt                      >  1:673655/1‑141 (MQ=255)                                                                                                                        tAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCAt                      >  2:293787/1‑141 (MQ=255)                                                                                                                             gaaTCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAg                 <  2:215716/140‑1 (MQ=255)                                                                                                                             cAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAg                 >  1:625149/1‑141 (MQ=255)                                                                                                                              aaacAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAAGTCGGCGTAATAATCCTTAATTAGAACGTTCATTAAAgg                <  2:221587/138‑1 (MQ=255)                                                                                                                              aaTCAGTAAAAAAAGAAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAgg                <  2:599154/141‑1 (MQ=255)                                                                                                                                tCAGTAAAAGATGCAAAAAATGGTAACCGCACGCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGAt              <  1:711679/141‑1 (MQ=255)                                                                                                                                 cAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATg             >  2:602202/1‑141 (MQ=255)                                                                                                                                      aaaGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATCATCCTGAATTAGAACGTTCATTAAAGGATGCgaaa        <  1:600302/141‑1 (MQ=255)                                                                                                                                      aaaGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGCgaaa        >  1:610571/1‑141 (MQ=255)                                                                                                                                       aaGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGCgaaag       >  1:200078/1‑141 (MQ=255)                                                                                                                                         gATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGGgaaagaa     >  2:496909/1‑141 (MQ=255)                                                                                                                                            gCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGCGAAAGAAACa  >  1:376229/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                            AGATGCAAATTTGAAAGAGGGGTAGTTTTTTATGGAAAATAAATTTGTTCCTGGTATTTTAATTGGTGCCGTAATTGGTGGTGCAATTAGTTTAGCTGATAAATCTACACGTCAAGCTTTAGTTCAATCAGTTAAAGATGCAAAAAATGGTAACCGCACTCGTAAGCCTTCTAAAGTCAGCAAGATTAAAGACGAAGTTTTATACTGGAAAGATGTTGTTGAAGAAATTCGTCGTAATAATCCTGAATTAGAACGTTCATTAAAGGATGCGAAAGAAACA  >  USA300TCH1516_ALE/2024367‑2024646

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