BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,789,551	C→T	9.2%	intergenic (‑32/+138)	gapA2 ← / ← coaE	Glyceraldehyde‑3‑phosphate dehydrogenase 2/Dephospho‑CoA kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,789,551	0	C	T	9.2%	302.0 / 16.0	98	intergenic (‑32/+138)	gapA2/coaE	Glyceraldehyde‑3‑phosphate dehydrogenase 2/Dephospho‑CoA kinase
Reads supporting (aligned to +/- strand): ref base C (38/51); new base T (3/6); total (41/57)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 101 total aligned reads displayed.

ATACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCATT  >  USA300TCH1516_ALE/1789461‑1789636

aTACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATAc                                                                                       >  1:3058779/1‑91 (MQ=255)

aTACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATAc                                                                                       >  1:3058969/1‑91 (MQ=255)

aTACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATAc                                                                                       >  1:4036067/1‑91 (MQ=255)

 tACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACt                                                                                      <  1:792648/91‑1 (MQ=255)

  aCGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACtt                                                                                     >  1:2054303/1‑91 (MQ=255)

  aCGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACtt                                                                                     >  1:2053628/1‑91 (MQ=255)

    gTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTa                                                                                   >  1:153754/1‑91 (MQ=255)

     tAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTAt                                                                                  >  1:1339527/1‑91 (MQ=255)

     tAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTAt                                                                                  <  1:260764/91‑1 (MQ=255)

      aaTACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATg                                                                                 <  1:2121133/91‑1 (MQ=255)

       aTACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATgg                                                                                <  1:4002515/91‑1 (MQ=255)

       aTACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATgg                                                                                >  1:2786353/1‑91 (MQ=255)

       aTACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATgg                                                                                >  1:2786833/1‑91 (MQ=255)

          ccATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAAt                                                                             >  1:641454/1‑91 (MQ=255)

               cttcCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGta                                                                        >  1:3614056/1‑91 (MQ=255)

               cttcCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGta                                                                        <  1:450280/91‑1 (MQ=255)

                 tcCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGtata                                                                      >  1:2724548/1‑91 (MQ=255)

                 tcCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGtata                                                                      >  1:2725047/1‑91 (MQ=255)

                  ccAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAg                                                                     <  1:3627474/91‑1 (MQ=255)

                  ccAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAg                                                                     <  1:3627336/91‑1 (MQ=255)

                    aaTTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCa                                                                   <  1:2823095/91‑1 (MQ=255)

                     aTTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCat                                                                  >  1:642701/1‑91 (MQ=255)

                       tCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCatat                                                                <  1:2821696/91‑1 (MQ=255)

                       tCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCatat                                                                <  1:2821670/91‑1 (MQ=255)

                           cccATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATtata                                                            >  1:3016849/1‑91 (MQ=255)

                           cccATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATtata                                                            <  1:2472927/91‑1 (MQ=255)

                            ccATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTatat                                                           <  1:2772902/91‑1 (MQ=255)

                            ccATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTatat                                                           >  1:3736941/1‑91 (MQ=255)

                            ccATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTatat                                                           <  1:3070537/91‑1 (MQ=255)

                             cATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATg                                                          <  1:3108857/91‑1 (MQ=255)

                             cATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATg                                                          <  1:3474710/91‑1 (MQ=255)

                              aTACCATTAATTGCAATATTCGTTGATATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGa                                                         <  1:2422619/91‑1 (MQ=255)

                              aTACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGa                                                         >  1:2749829/1‑91 (MQ=255)

                              aTACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGa                                                         <  1:2763981/91‑1 (MQ=255)

                              aTACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGa                                                         <  1:24120/91‑1 (MQ=255)

                                 ccATTATTTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAAt                                                      >  1:3226464/1‑91 (MQ=255)

                                 ccATTATTTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAAt                                                      >  1:3226386/1‑91 (MQ=255)

                                 ccATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAAt                                                      >  1:2882915/1‑91 (MQ=255)

                                        ttGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAg                                               <  1:1124985/91‑1 (MQ=255)

                                        ttGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAg                                               <  1:1462136/91‑1 (MQ=255)

                                        ttGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAg                                               <  1:1363424/91‑1 (MQ=255)

                                         tGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGc                                              <  1:2470668/91‑1 (MQ=255)

                                          gCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCg                                             <  1:3398550/91‑1 (MQ=255)

                                          gCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCg                                             <  1:1137590/91‑1 (MQ=255)

                                           cAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGt                                            <  1:3011271/91‑1 (MQ=255)

                                           cAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGt                                            >  1:1776836/1‑91 (MQ=255)

                                           cAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAAATAGTATAGCATATTATATGAAATTTGAAAGCGt                                            >  1:1776900/1‑91 (MQ=255)

                                             atatTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTCTGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGttt                                          <  1:3589015/91‑1 (MQ=255)

                                              tatTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGtttt                                         <  1:1890890/91‑1 (MQ=255)

                                              tatTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGtttt                                         <  1:2026087/91‑1 (MQ=255)

                                              tatTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGtttt                                         <  1:2027343/91‑1 (MQ=255)

                                                ttCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTcc                                       >  1:705456/1‑91 (MQ=255)

                                                 tCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCg                                      <  1:322688/91‑1 (MQ=255)

                                                   gTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGtt                                    <  1:1963654/91‑1 (MQ=255)

                                                    ttGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTg                                   >  1:3923855/1‑91 (MQ=255)

                                                    ttGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTg                                   >  1:3923761/1‑91 (MQ=255)

                                                     tGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGa                                  <  1:3585243/91‑1 (MQ=255)

                                                     tGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGa                                  >  1:543504/1‑91 (MQ=255)

                                                      gACATTAAAAGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGat                                 <  1:2668471/91‑1 (MQ=255)

                                                      gACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGat                                 >  1:2059630/1‑91 (MQ=255)

                                                         aTTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGattta                              >  1:869225/1‑91 (MQ=255)

                                                          ttAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttat                             <  1:1841178/91‑1 (MQ=255)

                                                          ttAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttat                             <  1:752960/91‑1 (MQ=255)

                                                          ttAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttat                             <  1:753586/91‑1 (MQ=255)

                                                            aaaaGCACCCCTTGATTTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAATTTTGAAAGCGTTTTCCGTTGatttagtt                           <  1:1458605/91‑1 (MQ=255)

                                                            aaaaGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttattt                           <  1:1548731/91‑1 (MQ=255)

                                                            aaaaGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttattt                           <  1:3105201/91‑1 (MQ=255)

                                                            aaaaGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGatttattt                           <  1:1723014/91‑1 (MQ=255)

                                                             aaaGCACCCCTTGATCTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGAtttattta                          <  1:1570494/91‑1 (MQ=255)

                                                             aaaGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGAtttattta                          <  1:46483/91‑1 (MQ=255)

                                                                 cACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCtt                      >  1:210564/1‑91 (MQ=255)

                                                                  aCCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTGTTTATCttt                     <  1:1898834/91‑1 (MQ=255)

                                                                  aCCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCttt                     >  1:1684692/1‑91 (MQ=255)

                                                                    cccTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCttttt                   >  1:1344481/1‑91 (MQ=255)

                                                                    cccTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCttttt                   >  1:1342982/1‑91 (MQ=255)

                                                                    cccTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCttttt                   >  1:1342925/1‑91 (MQ=255)

                                                                    cccTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCttttt                   >  1:1897050/1‑91 (MQ=255)

                                                                      cTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTta                 <  1:90700/91‑1 (MQ=255)

                                                                       ttGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTtat                <  1:459946/91‑1 (MQ=255)

                                                                         gATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATaa              >  1:1836070/1‑91 (MQ=255)

                                                                            cTTTATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATaaaaa           >  1:2344464/1‑91 (MQ=255)

                                                                            cTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATaaaaa           >  1:330250/1‑91 (MQ=255)

                                                                             tttATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATaaaaaa          >  1:673022/1‑91 (MQ=255)

                                                                              ttATGTGTTATATTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAAt         <  1:3017021/91‑1 (MQ=255)

                                                                              ttATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAAt         <  1:2876854/91‑1 (MQ=255)

                                                                              ttATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAAt         <  1:1023923/91‑1 (MQ=255)

                                                                              ttATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAAt         <  1:2876715/91‑1 (MQ=255)

                                                                               tATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATA        >  1:1695584‑M1/1‑91 (MQ=255)

                                                                                aTGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAG       >  1:1839642‑M1/1‑91 (MQ=255)

                                                                                aTGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAG       <  1:3161842‑M1/91‑1 (MQ=255)

                                                                                 tgtgTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGT      >  1:885324‑M1/1‑91 (MQ=255)

                                                                                 tgtgTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGT      <  1:152193‑M1/91‑1 (MQ=255)

                                                                                  gtgtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTC     <  1:1823781‑M1/91‑1 (MQ=255)

                                                                                  gtgtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTC     <  1:2032709‑M1/91‑1 (MQ=255)

                                                                                  gtgtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTC     <  1:2033887‑M1/91‑1 (MQ=255)

                                                                                    gtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCAT   <  1:1796779‑M1/91‑1 (MQ=255)

                                                                                    gtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCAT   <  1:3118388‑M1/91‑1 (MQ=255)

                                                                                    gtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCAT   >  1:2011073‑M1/1‑91 (MQ=255)

                                                                                    gtTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGCTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCAT   >  1:2010953‑M1/1‑91 (MQ=255)

                                                                                     ttATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCATT  >  1:2968953‑M1/1‑91 (MQ=255)

ATACGTAATACCATTCTTCCAATTCTACCCATACCATTAATTGCAATATTCGTTGACATTAAAAGCACCCCTTGATCTTTATGTGTTATACTTTATGGAATTAGTATAGCATATTATATGAAATTTGAAAGCGTTTTCCGTTGATTTATTTATCTTTTTTATAAAAAATAGTCATT  >  USA300TCH1516_ALE/1789461‑1789636

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