Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,673,843	T→C	6.1%	K391K (AAA→AAG)	dnaG ←	DNA primase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,673,843	0	T	C	6.1%	189.0 / 3.4	66	K391K (AAA→AAG)	dnaG	DNA primase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base C (2/2);  total (33/33)
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.11e-01

TAATAATTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGCATCGTTGCCATACTTACC  >  USA300TCH1516_ALE/1673703‑1673971                                                                                                                                             |                                                                                                                                 taataaTTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCt                                                                                                                                  >  1:862648/1‑141 (MQ=255)  aattatTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTTACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCtt                                                                                                                                 <  1:6319/141‑1 (MQ=255)  aagaatTGGCTGGTGCTGGATTGAATTTTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCAGTTAAAGCCTTTTGTTGCAAAATCGATGATGTCATAAGCGAAATATCATGACTTAGTTCtt                                                                                                                                 <  1:309786/138‑1 (MQ=255)     aatTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCa                                                                                                                              <  2:855365/141‑1 (MQ=255)      atTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCaa                                                                                                                             <  1:867305/141‑1 (MQ=255)        tGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAAt                                                                                                                           >  2:87801/1‑141 (MQ=255)        tGGCTGGTGCTTGATTGAATTGGATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAAt                                                                                                                           <  1:513594/141‑1 (MQ=255)          gCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATaa                                                                                                                         >  1:352889/1‑141 (MQ=255)              gTGCTTGATTGAATTGTAGTTCGTTAGATAATTTCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGtt                                                                                                                     <  2:109356/141‑1 (MQ=255)                        gaATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAg                                                                                                           >  2:225614/1‑141 (MQ=255)                          aTTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTTTTTTTTTTTTTTTGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATACAGTTCATTTGaatgtt                                                                                                         >  2:377532/1‑137 (MQ=255)                          aTTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAgga                                                                                                         >  2:215784/1‑140 (MQ=255)                          aTTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGt                                                                                                         >  1:95193/1‑141 (MQ=255)                           ttGTATTTCGTTAGATAATAGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATAACTTGTTTCTTCCAAATAACGTTCATATGAAAgagc                                                                                                        >  1:518215/1‑138 (MQ=255)                             ttATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCAt                                                                                                      <  2:329023/140‑1 (MQ=255)                             ggagttCGGTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCAt                                                                                                      <  1:366920/137‑1 (MQ=255)                             gTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCAt                                                                                                      <  1:3614/141‑1 (MQ=255)                             gTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCAGATGAAAGGTCAt                                                                                                      >  2:830223/1‑141 (MQ=255)                                 ttCGTTGGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTAtg                                                                                                  <  2:396238/141‑1 (MQ=255)                                    tttAGATAACTGCTCAGGACTAAAATTAAAAAAGGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGc                                                                                               <  2:33096/140‑1 (MQ=255)                                    gTTAGCTAATTGCTGAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTTGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGc                                                                                               <  1:673485/141‑1 (MQ=255)                                    gTTAGCTAATTGCTCAGGCCTAACATTGTAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGc                                                                                               <  1:446011/141‑1 (MQ=255)                                    gTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGc                                                                                               <  1:765520/141‑1 (MQ=255)                                    gTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGc                                                                                               <  1:865526/141‑1 (MQ=255)                                       aGCTAATTGCTCAGGAAAAAAAAAAAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAAt                                                                                            <  2:498905/141‑1 (MQ=255)                                       aGCTAATTGCTCAAGACTAACAATGAAAAATGGTGCAACATCATTTAAAGCTTTTTGTGGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAGATAACGTTCATATGAAAGGTCATTATGTGCAAt                                                                                            <  2:622694/141‑1 (MQ=255)                                         cTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCCTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAAttt                                                                                          >  2:864471/1‑141 (MQ=255)                                         cTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCCTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAAttt                                                                                          >  2:494243/1‑141 (MQ=255)                                         cTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCAAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAAttt                                                                                          >  2:427164/1‑141 (MQ=255)                                           aaTTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCa                                                                                        <  2:244161/141‑1 (MQ=255)                                                cTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCttt                                                                                   <  2:352889/141‑1 (MQ=255)                                                 tCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTCGTTCGTGCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCtttt                                                                                  >  1:776379/1‑141 (MQ=255)                                                  cAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTa                                                                                 >  2:763741/1‑141 (MQ=255)                                                  cAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTa                                                                                 <  2:435835/141‑1 (MQ=255)                                                  cAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTa                                                                                 >  2:412962/1‑141 (MQ=255)                                                  cAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCAGGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGCCATTATGTGCAATTTCATCTTTTa                                                                                 >  2:335023/1‑141 (MQ=255)                                                   aGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTaa                                                                                <  1:284910/141‑1 (MQ=255)                                                   aGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTaa                                                                                <  2:420512/141‑1 (MQ=255)                                                     gACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAACTCGAGTATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAata                                                                              <  1:358942/141‑1 (MQ=255)                                                           cATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCa                                                                        >  2:737351/1‑141 (MQ=255)                                                           cATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCa                                                                        >  2:268348/1‑141 (MQ=255)                                                                     tGGTGCAACATCATTTAAAGCCTTTTGTTGCAATATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATTACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAAtg                                                              <  2:289628/141‑1 (MQ=255)                                                                     tGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAACGACATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTTACTTTATAAtg                                                              >  1:793438/1‑141 (MQ=255)                                                                         tcAACATCATTTAAAGCCTTTTGTCGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCa                                                          <  1:617076/140‑1 (MQ=255)                                                                          cAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCaa                                                         <  2:352188/141‑1 (MQ=255)                                                                            aCATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAAt                                                       >  2:440205/1‑141 (MQ=255)                                                                                 atTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACtt                                                  <  2:381347/141‑1 (MQ=255)                                                                                      aaGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGt                                             >  1:263447/1‑141 (MQ=255)                                                                                       aGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGGAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATAAAAACATTATAATGTGCACCTGCCTTTTTGTc                                            <  2:624434/141‑1 (MQ=255)                                                                                           ttttGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCCTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCAttt                                        <  1:864471/141‑1 (MQ=255)                                                                                           ttttGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCCTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCAttt                                        <  1:494243/141‑1 (MQ=255)                                                                                            tttGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCAtttt                                       >  2:682627/1‑141 (MQ=255)                                                                                            tttGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCAtttt                                       >  2:682654/1‑141 (MQ=255)                                                                                              tgttgCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTCTAATGCGCAAATGACTTTTTGTCAtttttt                                     <  1:225614/141‑1 (MQ=255)                                                                                              tgttgCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCAtttttt                                     <  2:830004/141‑1 (MQ=255)                                                                                               gttgCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTa                                    >  1:585710/1‑141 (MQ=255)                                                                                               gttgCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTa                                    >  1:90949/1‑141 (MQ=255)                                                                                                   cAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAAGTACTTTTTGTAATTTTTTACaaa                                >  1:422441/1‑141 (MQ=255)                                                                                                        tCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAg                           >  1:519519/1‑141 (MQ=255)                                                                                                        tCGATGATTTCATAAGCGAAAGATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCCACTTATAATGTGCATATGCCTGTTTGGCATGTTTTACtaatacgg                           >  2:647891/1‑136 (MQ=255)                                                                                                          gatgatTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTa                         >  2:84022/1‑141 (MQ=255)                                                                                                            tgatTTCAGAAGCGAAATATCATGTCTGAGTTGTTTCAAATAACGTTCATATGTAAGGTCATTATGTGCAATTTCATCTTTTACTATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTaaa                       <  2:170071/141‑1 (MQ=255)                                                                                                               tttCATAAGCGAAATCTCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGc                    >  1:700865/1‑141 (MQ=255)                                                                                                               tttCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGc                    <  1:66316/141‑1 (MQ=255)                                                                                                                           aaTATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGCATCGTTGCCATa        <  1:427164/141‑1 (MQ=255)                                                                                                                           aaTATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAACCAGTAAATGCATCGTTGCCATa        >  2:543924/1‑141 (MQ=255)                                                                                                                            atatCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGCGCATCTCCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAGAGGCAGTAAATGCATCTTTATCATAc       >  1:410285/1‑141 (MQ=255)                                                                                                                                cATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGCATCGTTGCAATAATTAc   >  2:306845/1‑141 (MQ=255)                                                                                                                                 aTGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGCATCGTTGCCATACTTAcc  >  2:757649/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                 TAATAATTGGCTGGTGCTTGATTGAATTGTATTTCGTTAGCTAATTGCTCAGGACTAACATTGAAAAATGGTGCAACATCATTTAAAGCCTTTTGTTGCAAAATCGATGATTTCATAAGCGAAATATCATGACTTAGTTCTTTCAAATAACGTTCATATGAAAGGTCATTATGTGCAATTTCATCTTTTAATATACTCACTTTATAATGTGCAAATGACTTTTTGTCATTTTTTACAAAAGCAGTAAATGCATCGTTGCCATACTTACC  >  USA300TCH1516_ALE/1673703‑1673971

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