BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,972,377	T→C	9.0%	S24S (AGT→AGC)	prsA →	Foldase protein PrsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,972,377	0	T	C	9.0%	642.5 / 42.2	221	S24S (AGT→AGC)	prsA	Foldase protein PrsA
Reads supporting (aligned to +/- strand): ref base T (112/89); new base C (18/2); total (130/91)
Fisher's exact test for biased strand distribution p-value = 3.28e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 226 total aligned reads displayed.

TTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCAA  >  USA300TCH1516_ALE/1972287‑1972467

   cAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGcc                                                                                         <  1:309483/91‑1 (MQ=255)

   cAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGcc                                                                                         >  1:2953404/1‑91 (MQ=255)

           tttCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGActc                                                                                 >  1:224880/1‑91 (MQ=255)

           tttCAACTATGAAGATGATAAACAAATTAATCGTCCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGActc                                                                                 <  1:564925/91‑1 (MQ=255)

             tCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTa                                                                               <  2:3725207/91‑1 (MQ=255)

               aaCTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTaaa                                                                             <  2:502146/91‑1 (MQ=255)

                  gaggAAGATGATAAACTAAGTAATCGTTCCGGGAACAGCTAGTGCTTTATTATTATTCTCTTTTGGCGCTAGTGCCACAGACTCTAAAGaa                                                                          <  1:2524050/88‑1 (MQ=255)

                    tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                        <  2:698007/91‑1 (MQ=255)

                    tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                        >  2:1472573/1‑91 (MQ=255)

                    tGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGaaaa                                                                        >  2:39774/1‑91 (MQ=255)

                     gAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAAt                                                                       <  1:1025356/91‑1 (MQ=255)

                     gAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAAt                                                                       >  1:1798958/1‑91 (MQ=255)

                      aaGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATa                                                                      <  1:3382864/91‑1 (MQ=255)

                        gatgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGTCTCTAAAGAAAATACa                                                                    >  2:2973872/1‑91 (MQ=255)

                        gatgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAATTaat                                                                    >  1:3570043/1‑89 (MQ=255)

                        gatgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACa                                                                    >  2:3208779/1‑91 (MQ=255)

                         atgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAt                                                                   >  1:1705664/1‑91 (MQ=255)

                          tgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAtt                                                                  <  2:3812567/91‑1 (MQ=255)

                          tgatAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACAtt                                                                  >  1:2611857/1‑91 (MQ=255)

                            atAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                                >  2:1607816/1‑91 (MQ=255)

                            atAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                                <  2:1578789/91‑1 (MQ=255)

                            atAAACAAATTAAACGTTCCGGTAACAGCTAGTGCTTTATTATGAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTaa                                                                <  1:3874374/91‑1 (MQ=255)

                             tAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAt                                                               <  1:2171811/91‑1 (MQ=255)

                              aaaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAtt                                                              >  1:1979776/1‑91 (MQ=255)

                              aaaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAtt                                                              <  2:1864036/91‑1 (MQ=255)

                              aaaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCGTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAtt                                                              <  2:2054366/91‑1 (MQ=255)

                               aaCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAAttt                                                             >  1:623801/1‑91 (MQ=255)

                                aCAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCTCTAGTGCCACAGACTCTAAAGAAAATACATTAATttc                                                            <  1:3259259/91‑1 (MQ=255)

                                 cAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttct                                                           >  1:3295123/1‑91 (MQ=255)

                                 cAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATttct                                                           <  1:986654/91‑1 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGTTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                        <  2:3417615/91‑1 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCAGAGACTCTAAAGAAAATACATTAATTtcttct                                                        <  2:2968172/91‑1 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                        >  1:3852645/1‑91 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                        >  1:2049435/1‑91 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                        <  2:2598512/91‑1 (MQ=255)

                                    aTTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTtcttct                                                        >  2:2219267/1‑91 (MQ=255)

                                     ttAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTa                                                       <  2:60950/91‑1 (MQ=255)

                                      tAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTaa                                                      <  2:1045019/91‑1 (MQ=255)

                                        aTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAg                                                    >  1:2389252/1‑91 (MQ=255)

                                        aTCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAg                                                    >  2:2544496/1‑91 (MQ=255)

                                         tCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGc                                                   >  1:1786550/1‑91 (MQ=255)

                                          cGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCt                                                  <  2:1517972/91‑1 (MQ=255)

                                            ttCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTgg                                                <  1:1441116/91‑1 (MQ=255)

                                              ccGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGgag                                              <  1:3210714/91‑1 (MQ=255)

                                               cGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGgag                                              <  1:2023370/90‑1 (MQ=255)

                                               cGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGgag                                              >  2:2023370/1‑90 (MQ=255)

                                                ggTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGAc                                            >  2:683995/1‑91 (MQ=255)

                                                 gTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACg                                           <  1:3036030/91‑1 (MQ=255)

                                                 gTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACg                                           >  2:1984856/1‑91 (MQ=255)

                                                  tAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGt                                          <  2:2811438/91‑1 (MQ=255)

                                                  tAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGt                                          <  2:2333885/91‑1 (MQ=255)

                                                    aCAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTaa                                        >  2:2649033/1‑91 (MQ=255)

                                                     cAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAAc                                       >  1:463549/1‑91 (MQ=255)

                                                       gCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAg                                     >  1:3444148/1‑91 (MQ=255)

                                                        cTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                                    <  1:2492730/91‑1 (MQ=255)

                                                        cTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGt                                    >  1:2710768/1‑91 (MQ=255)

                                                         tAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGtt                                   >  2:2849395/1‑91 (MQ=255)

                                                         tAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGtt                                   >  1:814640/1‑91 (MQ=255)

                                                          aGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTg                                  >  2:482453/1‑91 (MQ=255)

                                                          aGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTg                                  <  2:816641/91‑1 (MQ=255)

                                                          aGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTg                                  >  1:2708318/1‑91 (MQ=255)

                                                             gCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAg                               >  1:324113/1‑91 (MQ=255)

                                                             gCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAg                               <  1:3224969/91‑1 (MQ=255)

                                                              cTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGa                              <  1:2777865/91‑1 (MQ=255)

                                                              cTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGa                              >  2:2071288/1‑91 (MQ=255)

                                                               tttATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGAt                             <  1:3652592/91‑1 (MQ=255)

                                                                 tattattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATAc                           >  2:2415277/1‑91 (MQ=255)

                                                                   ttattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACaa                         >  1:1903573/1‑91 (MQ=255)

                                                                    tattaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAAt                        >  2:1258141/1‑91 (MQ=255)

                                                                      ttaGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGa                      <  1:1149451/91‑1 (MQ=255)

                                                                       taGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaa                     >  1:3510684/1‑91 (MQ=255)

                                                                        aGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaa                    >  1:2075910/1‑91 (MQ=255)

                                                                         ggCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaa                   <  1:1174829/91‑1 (MQ=255)

                                                                         ggCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaa                   >  1:1372931/1‑91 (MQ=255)

                                                                          gcgcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaa                  <  1:3494664/91‑1 (MQ=255)

                                                                           cgcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATGCAATGaaaaaa                 >  2:354957/1‑91 (MQ=255)

                                                                            gcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaaa                >  2:1311411/1‑91 (MQ=255)

                                                                            gcTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGaaaaaaa                >  2:3725352/1‑91 (MQ=255)

                                                                             cTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAAt               >  1:3148745/1‑91 (MQ=255)

                                                                              ttGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATc              >  1:3156588/1‑91 (MQ=255)

                                                                               tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg             >  2:1900047/1‑91 (MQ=255)

                                                                               tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg             <  1:797549/91‑1 (MQ=255)

                                                                               tgtgGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCg             <  1:1904749/91‑1 (MQ=255)

                                                                                gtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCgg            >  2:3046753/1‑91 (MQ=255)

                                                                                gtgGCACTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCgg            >  1:115216/1‑91 (MQ=255)

                                                                                  ggcacTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa          >  2:3136073/5‑91 (MQ=255)

                                                                                  ggcacTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa          >  1:2317862/5‑91 (MQ=255)

                                                                                  ggCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa          <  2:2099643/91‑1 (MQ=255)

                                                                                  ggCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa          >  1:2874209/1‑91 (MQ=255)

                                                                                  ggCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTa          >  1:2086716/1‑91 (MQ=255)

                                                                                   gcacTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTaa         >  1:1300925/4‑91 (MQ=255)

                                                                                    cacTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTaaa        >  1:72065/3‑91 (MQ=255)

                                                                                    cacTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTaaa        >  1:2436521/3‑91 (MQ=255)

                                                                                      cTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGa      <  1:1478695/91‑1 (MQ=255)

                                                                                      cTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGa      >  1:1769177/1‑91 (MQ=255)

                                                                                      cTAGCGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGa      >  1:2917478/1‑91 (MQ=255)

                                                                                         gTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCa   <  1:2681054/91‑1 (MQ=255)

                                                                                          tGCCACAGACTCTACAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCaa  >  2:3297724/1‑91 (MQ=255)

                                                                                          tGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCaa  >  2:273955/1‑91 (MQ=255)

                                                                                          tGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCaa  <  1:2968463/91‑1 (MQ=255)

TTACAAAGGAGTTTCAACTATGAAGATGATAAACAAATTAATCGTTCCGGTAACAGCTAGTGCTTTATTATTAGGCGCTTGTGGCGCTAGTGCCACAGACTCTAAAGAAAATACATTAATTTCTTCTAAAGCTGGAGACGTAACAGTTGCAGATACAATGAAAAAAATCGGTAAAGATCAA  >  USA300TCH1516_ALE/1972287‑1972467

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