Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,235,262	G→A	6.1%	Q156* (CAA→TAA)	tdk ←	Thymidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,235,262	0	G	A	6.1%	196.1 / 3.0	66	Q156* (CAA→TAA)	tdk	Thymidine kinase
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base A (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 = 6.82e-01

AACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAATAACAC  >  USA300TCH1516_ALE/2235123‑2235399                                                                                                                                            |                                                                                                                                          aaCACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTg                                                                                                                                           <  1:507356/140‑1 (MQ=255)   cacTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTTGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCg                                                                                                                                         <  2:105625/140‑1 (MQ=255)    acTTGTAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTTGGATAAGCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGt                                                                                                                                        <  1:290313/140‑1 (MQ=255)              cctcctTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGAt                                                                                                                              <  1:245418/140‑1 (MQ=255)               ctcctTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATg                                                                                                                             <  1:733663/140‑1 (MQ=255)                    tattattatCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCTATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATgaagaa                                                                                                                        <  1:719575/140‑1 (MQ=255)                     attattatCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAAc                                                                                                                       <  1:593800/140‑1 (MQ=255)                         ttatCACTAGGTGCTACGAACTCAGACGCACAACAACGAGGCTAATAACATTCATTGGTACCTACGAATATGATCTGATCAACGATTTCAGATGGTTTAAGGGTTATTAATTGATGAGTATTATGAGATgaagagtgaca                                                                                                                   >  2:488977/1‑134 (MQ=255)                            tCACTAGGTGCTACGATATGGTGCGAGCGACAACGTTGCTCATAACTTTCATGGGCACCTGAGAAGATAATCGGATCATCGATTATAGCTGGTTTACCATTTATTAATCGTAGCGTACTAATAAAAGAAGACCCAaaaaa                                                                                                                >  2:70627/1‑139 (MQ=255)                               cTAGGTGCTACGATATGGTGCGCTCGACAACGTGGCTCATAAATTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAAAAGc                                                                                                             >  2:424862/1‑140 (MQ=255)                                      cTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACt                                                                                                      <  2:344639/140‑1 (MQ=255)                                          gATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTTGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCtt                                                                                                  <  1:756750/140‑1 (MQ=255)                                          gATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCtt                                                                                                  <  1:751599/140‑1 (MQ=255)                                            tatGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCATTGCGTTCTACTAGATGAAGCACCACAAACAGCGCAACCTGCCTGa                                                                                                >  2:304099/1‑140 (MQ=255)                                               ggtccccgcgACAACGTGGCTCATAACTTGCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGt                                                                                             <  1:458025/134‑1 (MQ=255)                                                gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAggt                                                                                            >  1:725426/1‑138 (MQ=255)                                                gTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGtt                                                                                            >  1:150211/1‑140 (MQ=255)                                                  gcgcgcGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTc                                                                                          >  2:335958/1‑140 (MQ=255)                                                   cgcgcgACAACGTGGCTTATAACTTTCATTGGCTCCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATGTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCg                                                                                         <  1:196788/140‑1 (MQ=255)                                                   cgcgcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCg                                                                                         <  1:171098/140‑1 (MQ=255)                                                   cgcgcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCg                                                                                         >  1:565239/1‑140 (MQ=255)                                                      gcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTa                                                                                      <  1:355509/140‑1 (MQ=255)                                                      gcgACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATATATTAATCGTTGCGATCTACTAGATGAACAACCCCAACCAGCACAAACCGCTACAAGTTTCGtac                                                                                      >  2:206108/1‑138 (MQ=255)                                                      gcgACAACGTGGCTCATAACTTTCAGTGGCACCAACAAAGATATTAAGGTATTCGTTTTGAGCTGGTTTACCATTTATTAAGCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTa                                                                                      <  1:371934/140‑1 (MQ=255)                                                      cggacgacGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTa                                                                                      <  2:266817/138‑1 (MQ=255)                                                        gACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACt                                                                                    >  1:615228/1‑140 (MQ=255)                                                           aaCGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTCGATGAAGAACCACAAAAAGCACAAACTGCTTGACGTTGCGTCACTCGt                                                                                 >  2:565523/1‑140 (MQ=255)                                                           aaCGTGGCTCATAACTTGCATTGGCACCTGCTAAGATAATCGGATCATCGATTGTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGt                                                                                 <  1:753900/140‑1 (MQ=255)                                                             cGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACACACTGCTTGAAGTTTCGTTACTTGTTc                                                                               >  1:79878/1‑140 (MQ=255)                                                                ggCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACt                                                                            <  2:739616/140‑1 (MQ=255)                                                                            tCATTGGCACCTACTAAGATAATCGGATCATCGATGTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCaa                                                                >  1:419739/1‑140 (MQ=255)                                                                               ttGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTACGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAAttt                                                             >  2:244755/1‑140 (MQ=255)                                                                               ttGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTACGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAAttt                                                             >  2:209534/1‑140 (MQ=255)                                                                                tGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAAAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAACTTa                                                            >  1:588490/1‑140 (MQ=255)                                                                                  gCACCTACTAAGATAATCGGATCATCTATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAgg                                                          <  2:66681/140‑1 (MQ=255)                                                                                   cACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGc                                                         >  1:356616/1‑140 (MQ=255)                                                                                   cACCTACTAAGATAATCGGATAATCGATTTTAGCTGGTTTACCATTTATAAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGc                                                         >  2:534992/1‑140 (MQ=255)                                                                                          taAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATAAATATAAGCATTGGtt                                                  >  1:226849/1‑140 (MQ=255)                                                                                          taAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGACATCAATTTAGGCATTGGtt                                                  >  2:596181/1‑140 (MQ=255)                                                                                             gATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGa                                               >  2:262396/1‑140 (MQ=255)                                                                                                aaccGGATCAGCGATTTTAGCTGGTTTACAATTTAGTAGTCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACg                                            <  1:110057/137‑1 (MQ=255)                                                                                                aaTCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACg                                            >  2:748725/1‑140 (MQ=255)                                                                                                aaTCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACg                                            >  1:525750/1‑140 (MQ=255)                                                                                                   cGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTATACTAGATGTAGACCCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGtt                                         >  1:377752/1‑140 (MQ=255)                                                                                                      atcatcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGc                                      >  2:477987/1‑140 (MQ=255)                                                                                                      atcatcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGc                                      >  1:35786/1‑140 (MQ=255)                                                                                                       tcatcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTCACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGca                                     >  2:219449/1‑139 (MQ=255)                                                                                                        catcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCACTTTAGGCATTGGTTCGAACGGTTCGccc                                    >  1:386679/1‑140 (MQ=255)                                                                                                        catcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGccc                                    <  2:338035/140‑1 (MQ=255)                                                                                                        catcGATTTTAGCTCTTTTACCATTTATTAATCGCTGCGTTCTACTAGATGAAGAACCACCAACAGCACAAACTGCTCGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGccc                                    <  2:454239/140‑1 (MQ=255)                                                                                                         atcGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGcccc                                   >  1:637123/1‑140 (MQ=255)                                                                                                                  aGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCGCTCAAATCCa                          >  2:122416/1‑140 (MQ=255)                                                                                                                   gCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACACCCATCAATTTAGGCATTGGTTCGAACGCTTCGCCCCTAAAATCCat                         >  2:578755/1‑140 (MQ=255)                                                                                                                      ggTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATc                      >  1:451105/1‑140 (MQ=255)                                                                                                                       gTTTAACATTTATTACTCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACGTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCt                     <  2:298693/140‑1 (MQ=255)                                                                                                                        ttcaCCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTATTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTa                    <  1:424862/137‑1 (MQ=255)                                                                                                                          tACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAAt                  >  1:242589/1‑140 (MQ=255)                                                                                                                          tACCATGTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAGCTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAAt                  <  2:27722/140‑1 (MQ=255)                                                                                                                            ccATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATcc                >  2:416221/1‑140 (MQ=255)                                                                                                                                ttattaATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGca            >  2:155521/1‑140 (MQ=255)                                                                                                                                 tattaATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGcaa           >  2:794987/1‑140 (MQ=255)                                                                                                                                    aaaTCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGcaacaa        <  1:534992/139‑1 (MQ=255)                                                                                                                                     aaTCGTTACGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAAt       <  1:209534/140‑1 (MQ=255)                                                                                                                                     aaTCGTTACGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAAt       <  1:244755/140‑1 (MQ=255)                                                                                                                                       tCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAATaa     >  1:774376/1‑140 (MQ=255)                                                                                                                                       tCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAATaa     >  1:560177/1‑140 (MQ=255)                                                                                                                                          ttGCGTTCTACTAGATGAAGAACCACAAACAGCACAACCTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCTATATCTAATCCAGCAACAATAacac  <  1:396162/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          AACACTTATAATTCCTCCTTATTATTATCACTAGGTGCTACGATATGGTGCGCGCGACAACGTGGCTCATAACTTTCATTGGCACCTACTAAGATAATCGGATCATCGATTTTAGCTGGTTTACCATTTATTAATCGTTGCGTTCTACTAGATGAAGAACCACAAACAGCACAAACTGCTTGAAGTTTCGTTACTTGTTCACTGACAGCCATCAATTTAGGCATTGGTTCGAACGGTTCGCCCCTAAAATCCATATCTAATCCAGCAACAATAACAC  >  USA300TCH1516_ALE/2235123‑2235399

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