Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	591,934	Δ1 bp	5.3%	coding (2150/3624 nt)	rpoC →	DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	591,934	0	A	.	5.3%	258.2 / 7.7	76	coding (2150/3624 nt)	rpoC	DNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (44/28);  new base . (2/2);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCTGGTAAGATTA  >  USA300TCH1516_ALE/591798‑592071                                                                                                                                         |                                                                                                                                          gATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGCCAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGt                                                                                                                                       >  1:185315/1‑141 (MQ=255) gATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGt                                                                                                                                       >  1:237478/1‑141 (MQ=255) gATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGt                                                                                                                                       >  1:169912/1‑141 (MQ=255) gATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGt                                                                                                                                       >  1:173614/1‑141 (MQ=255) gATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGt                                                                                                                                       >  1:18490/1‑141 (MQ=255)   tAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGa                                                                                                                                     >  2:144050/1‑141 (MQ=255)     aGCAACAACTACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAAt                                                                                                                                   >  2:215419/1‑141 (MQ=255)     aGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAAt                                                                                                                                   >  2:443839/1‑141 (MQ=255)      gCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAAtt                                                                                                                                  >  1:42378/1‑141 (MQ=255)      gCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAAtt                                                                                                                                  >  1:292373/1‑141 (MQ=255)       caacaaATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTg                                                                                                                                 >  2:367461/1‑141 (MQ=255)          caaATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATg                                                                                                                              <  2:155965/141‑1 (MQ=255)          caaATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATg                                                                                                                              <  2:439315/141‑1 (MQ=255)          caaATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATg                                                                                                                              <  2:226236/141‑1 (MQ=255)                 ttGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGCAGGAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCat                                                                                                                       <  1:436028/141‑2 (MQ=255)                 ttGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCAc                                                                                                                       <  2:121375/141‑1 (MQ=255)                    atgaGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTg                                                                                                                    >  2:318946/1‑141 (MQ=255)                    atgaGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTg                                                                                                                    >  1:316998/1‑141 (MQ=255)                        gCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATaa                                                                                                                >  1:335600/1‑141 (MQ=255)                            gAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGCAATTTGGACAGATGCAAAAGATCAAATTAAAGGTGAATTGATGCAATCACTTGATAAAACt                                                                                                            >  1:92611/1‑141 (MQ=255)                            gAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACt                                                                                                            <  1:6573/141‑1 (MQ=255)                             aaaaaTTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGCGAATTGATGCAATCACTTGATAAAACTa                                                                                                           >  2:31691/1‑141 (MQ=255)                               aaaTTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAAc                                                                                                         >  2:240406/1‑141 (MQ=255)                                  ttAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCa                                                                                                      <  1:131979/141‑1 (MQ=255)                                  ttAGTCGAAAGAATTACAAAACAATTCAACCGTGGTTTAATCACAGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCa                                                                                                      <  1:357663/141‑1 (MQ=255)                                         aCAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCa                                                                                               <  1:440699/141‑1 (MQ=255)                                               ttACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATTTTCAtgatga                                                                                         >  1:144815/1‑141 (MQ=255)                                               ttACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTC‑AGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAg                                                                                        >  2:68397/1‑141 (MQ=255)                                               ttACAAAACAATTCAACCGCGGTATAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTC‑AGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAg                                                                                        >  2:460461/1‑141 (MQ=255)                                                   aaaaCAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGa                                                                                     <  2:259003/141‑1 (MQ=255)                                                    aaaCAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGAt                                                                                    <  2:42879/141‑1 (MQ=255)                                                        aaTTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAg                                                                                <  2:100362/141‑1 (MQ=255)                                                              aCCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGGTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCGTCATGAGGAGTGAGTCAGGTGccc                                                                          >  2:160063/1‑141 (MQ=255)                                                                  tGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTgg                                                                      <  2:214933/141‑1 (MQ=255)                                                                      ttAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAAc                                                                  >  2:367310/1‑141 (MQ=255)                                                                      ttAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAAc                                                                  >  2:257921/1‑141 (MQ=255)                                                                      ttAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAAc                                                                  >  1:103879/1‑141 (MQ=255)                                                                      ttAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAAc                                                                  >  1:29753/1‑141 (MQ=255)                                                                      ttAAGCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTTAAGGTGAATTGATGCAATTACTTGATAAAACTAACCCAATTTTCATGATGAGTGATTCAGGTGCCCGTGGTacc                                                                  >  2:253056/1‑139 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:446257/141‑1 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:185315/141‑1 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:169912/141‑1 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:73563/141‑1 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:117563/141‑1 (MQ=255)                                                                       tAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTGCATGATGAGTGATTCAGGTGCCCGTGGTAACg                                                                 <  2:18490/141‑1 (MQ=255)                                                                         aTCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCa                                                               <  2:453216/141‑1 (MQ=255)                                                                         aTCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCa                                                               <  1:81671/141‑1 (MQ=255)                                                                            aCTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGTCCGTGGTAACGCATCt                                                            >  1:80711/1‑141 (MQ=255)                                                                                  gaagaaAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACttt                                                      <  1:167498/141‑1 (MQ=255)                                                                                      aaagaTATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTacac                                                  >  1:442828/1‑141 (MQ=255)                                                                                        agaTATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAg                                                <  2:229398/141‑1 (MQ=255)                                                                                         gaTATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGt                                               >  1:57154/1‑141 (MQ=255)                                                                                         gaTATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGt                                               >  1:277011/1‑141 (MQ=255)                                                                                              aaTGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGGGGTAACGCATCTAGCTTTACAGAGTTAGCa                                          >  1:15551/1‑141 (MQ=255)                                                                                              aaTGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGGGATTCAGGTGCCCGTGGTCACGCATCTAACTTTACACAGTTAGCa                                          >  1:330985/1‑141 (MQ=255)                                                                                                       gttgAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGt                                 >  2:300085/1‑141 (MQ=255)                                                                                                        ttgAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTg                                <  2:316998/141‑1 (MQ=255)                                                                                                        ttgAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTg                                <  1:318946/141‑1 (MQ=255)                                                                                                          gAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGa                              >  1:157560/1‑141 (MQ=255)                                                                                                          gAAATTTGGACAGATGCAAAAGATCAAATTC‑AGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGAt                             <  1:460461/141‑1 (MQ=255)                                                                                                          gAAATTTGGACAGATGCAAAAGATCAAATTC‑AGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGAt                             <  1:68397/141‑1 (MQ=255)                                                                                                           aaaTTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGAt                             >  1:15178/1‑141 (MQ=255)                                                                                                                 ggACAGCTGCAACAGATCAAATTCAAGGTGAATTGATGCAATCACATGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATgg                       >  2:112664/1‑141 (MQ=255)                                                                                                                 ggACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATgg                       >  1:227829/1‑141 (MQ=255)                                                                                                                 ggACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATgg                       >  2:245346/1‑141 (MQ=255)                                                                                                                   aCAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGcc                     >  2:415084/1‑141 (MQ=255)                                                                                                                   aCAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGcc                     >  1:166623/1‑141 (MQ=255)                                                                                                                          cAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCAt              >  2:354212/1‑141 (MQ=255)                                                                                                                          cAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGAGTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCAt              >  2:434872/1‑141 (MQ=255)                                                                                                                            aaaGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCt            >  1:211012/1‑141 (MQ=255)                                                                                                                            aaaGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCt            >  1:181035/1‑141 (MQ=255)                                                                                                                            aaaGAGCAAATGCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGAATCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCt            <  2:369646/141‑1 (MQ=255)                                                                                                                                  cAAATTCAAGGTGAATTGATGCAAGCACTTGATAAAACTAACCCAATCTTAATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCAGCACCATCTGGTAAg      >  2:193629/1‑141 (MQ=255)                                                                                                                                     aTTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCTGGTAAGAtt   <  2:292373/141‑1 (MQ=255)                                                                                                                                      ttCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCTGGTAAGAttt  >  1:27003/1‑140 (MQ=255)                                                                                                                                      ttCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCTGGTAAGATTa  <  2:355206/141‑1 (MQ=255)                                                                                                                                         |                                                                                                                                          GATAAGCAACAAATACTTGATGAGCATGAAAAATTAGTCGACAGAATTACAAAACAATTCAACCGTGGTTTAATCACTGAAGAAGAAAGATATAATGCAGTTGTTGAAATTTGGACAGATGCAAAAGATCAAATTCAAGGTGAATTGATGCAATCACTTGATAAAACTAACCCAATCTTCATGATGAGTGATTCAGGTGCCCGTGGTAACGCATCTAACTTTACACAGTTAGCAGGTATGCGTGGATTGATGGCCGCACCATCTGGTAAGATTA  >  USA300TCH1516_ALE/591798‑592071

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑