BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,497,391	C→G	9.1%	N147K (AAC→AAG)	USA300HOU_RS12865 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	2,497,391	0	C	G	9.1%	450.1 / 27.4	131	N147K (AAC→AAG)	USA300HOU_RS12865	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (54/65); new base G (7/5); total (61/70)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 131 total aligned reads displayed.

AATTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGT  >  USA300TCH1516_ALE/2497302‑2497480

aaTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACt                                                                                          <  1:2847680/91‑1 (MQ=255)

 aTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACtt                                                                                         <  1:61722/91‑1 (MQ=255)

 aTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACtt                                                                                         <  1:2410650/91‑1 (MQ=255)

 aTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACtt                                                                                         >  1:28497/1‑91 (MQ=255)

  tttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTa                                                                                        <  1:3471749/91‑1 (MQ=255)

  tttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTa                                                                                        <  1:1568398/91‑1 (MQ=255)

   ttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTAc                                                                                       >  1:1900074/1‑91 (MQ=255)

   ttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTAc                                                                                       <  1:3962038/91‑1 (MQ=255)

     aCCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGa                                                                                     <  1:2002810/91‑1 (MQ=255)

      ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGAt                                                                                    <  1:1882074/91‑1 (MQ=255)

      ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGAt                                                                                    <  1:3064384/91‑1 (MQ=255)

       cTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAtt                                                                                   <  1:1017890/91‑1 (MQ=255)

       cTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAtt                                                                                   >  1:2248866/1‑91 (MQ=255)

       cTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAtt                                                                                   >  1:256849/1‑91 (MQ=255)

        tAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGATTa                                                                                  <  1:1168151/91‑1 (MQ=255)

        tAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATATTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTa                                                                                  >  1:2396438/1‑91 (MQ=255)

              tGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaat                                                                            <  1:2969244/91‑1 (MQ=255)

                atatTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaataa                                                                          >  1:1514542/1‑91 (MQ=255)

                  atTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATaaaa                                                                        >  1:1762755/1‑91 (MQ=255)

                   ttAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAg                                                                       <  1:3833530/91‑1 (MQ=255)

                    tAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGt                                                                      <  1:1352753/91‑1 (MQ=255)

                    tAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGt                                                                      <  1:583804/91‑1 (MQ=255)

                    gtagttaaaatttctaatgcaaatctgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGt                                                                      >  1:41760‑M2/42‑91 (MQ=255)

                    gtagttaaaatttctaatgcaaatctgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGt                                                                      >  1:42959‑M2/42‑91 (MQ=255)

                     aaTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTa                                                                     <  1:293429/91‑1 (MQ=255)

                     aaTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTa                                                                     >  1:1272855/1‑91 (MQ=255)

                     aaTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTAGTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTa                                                                     >  1:3299865/1‑91 (MQ=255)

                      aTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAg                                                                    <  1:1071583/91‑1 (MQ=255)

                      aTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAg                                                                    <  1:1071584/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:2921555/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:1606520/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:16168/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:1416043/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:2182884/91‑1 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   >  1:1797392/1‑91 (MQ=255)

                       tttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                   <  1:2257651/91‑1 (MQ=255)

                        ttaaaatttctaatgcaaatctgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGaa                                                                  >  1:2586052‑M2/38‑91 (MQ=255)

                        TTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGaa                                                                  <  1:304647‑M1/91‑1 (MQ=255)

                         TAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAAt                                                                 >  1:503258‑M1/1‑91 (MQ=255)

                          aaaatttctaatgcaaatctgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAAtt                                                                >  1:781647‑M2/36‑91 (MQ=255)

                          aaaatttctaatgcaaatctgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAAtt                                                                >  1:1457249‑M2/36‑91 (MQ=255)

                           AAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTc                                                               <  1:3813371‑M1/91‑1 (MQ=255)

                           AAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTc                                                               <  1:3320577‑M1/91‑1 (MQ=255)

                             ATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGt                                                             >  1:3279108‑M1/1‑91 (MQ=255)

                              TTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTa                                                            >  1:2465156‑M1/1‑91 (MQ=255)

                               ttctaatgcaaaactgctgtgaaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCatat                                                           >  1:3208290‑M2/31‑91 (MQ=255)

                                TCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTAtt                                                          <  1:3878192‑M1/91‑1 (MQ=255)

                                 CACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGATCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTg                                                         <  1:3592690‑M1/91‑1 (MQ=255)

                                 CACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTg                                                         <  1:1319707‑M1/91‑1 (MQ=255)

                                  acaGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTga                                                        >  1:861533/1‑91 (MQ=255)

                                   cagcagATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTgag                                                       >  1:2808287/1‑91 (MQ=255)

                                       agATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaaag                                                   <  1:2528999/91‑1 (MQ=255)

                                       agATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaaag                                                   <  1:915061/91‑1 (MQ=255)

                                       agATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaaag                                                   >  1:3745172/1‑91 (MQ=255)

                                        gATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaaaga                                                  <  1:1323949/91‑1 (MQ=255)

                                           ttGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACt                                               <  1:4009424/91‑1 (MQ=255)

                                           ttGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACt                                               <  1:4011143/91‑1 (MQ=255)

                                             gCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTaa                                             >  1:3607653/1‑91 (MQ=255)

                                                  aGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAAt                                        <  1:3353311/91‑1 (MQ=255)

                                                    aaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATaa                                      >  1:35113‑M2/10‑91 (MQ=255)

                                                    aaattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATaa                                      >  1:35727‑M2/10‑91 (MQ=255)

                                                     aattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAAc                                     <  1:1433447‑M2/83‑1 (MQ=255)

                                                     aattttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAAc                                     <  1:3205204‑M2/83‑1 (MQ=255)

                                                        tttaaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACAtt                                  <  1:865286‑M2/86‑1 (MQ=255)

                                                          taaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTAc                                <  1:4132100‑M2/88‑1 (MQ=255)

                                                            attactaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                              <  1:1398446/89‑1 (MQ=255)

                                                            aCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                              >  1:3031935‑M2/2‑91 (MQ=255)

                                                            aCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                              >  1:1497515‑M2/2‑91 (MQ=255)

                                                             ctactaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaa                             >  1:3047691/1‑91 (MQ=255)

                                                             ctactaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaa                             >  1:390756/1‑91 (MQ=255)

                                                              tactaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaac                            <  1:2199608/91‑1 (MQ=255)

                                                               actaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTcaca                           >  1:1920895/1‑91 (MQ=255)

                                                               actaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaaca                           >  1:3479536/1‑91 (MQ=255)

                                                                ctaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaacaa                          <  1:108130/91‑1 (MQ=255)

                                                                 taTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaacaaa                         <  1:416625/91‑1 (MQ=255)

                                                                  aTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaacaaac                        <  1:1422439/91‑1 (MQ=255)

                                                                   ttGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaacaaaca                       <  1:2112313/91‑1 (MQ=255)

                                                                      ccGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACt                    <  1:3341772/91‑1 (MQ=255)

                                                                       cGACCATGATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACtt                   <  1:2077026/91‑1 (MQ=255)

                                                                        gACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACttt                  <  1:1725503/91‑1 (MQ=255)

                                                                         aCCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACtttt                 <  1:2976303/91‑1 (MQ=255)

                                                                          ccATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTc                >  1:2423513/1‑91 (MQ=255)

                                                                          ccATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTc                <  1:894438/91‑1 (MQ=255)

                                                                               ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt           >  1:835150/1‑91 (MQ=255)

                                                                               ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt           <  1:947129/91‑1 (MQ=255)

                                                                               ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt           <  1:947328/91‑1 (MQ=255)

                                                                               ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt           >  1:3533599/1‑91 (MQ=255)

                                                                               ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt           >  1:2234844/1‑91 (MQ=255)

                                                                                taATCGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTa          >  1:398845/1‑91 (MQ=255)

                                                                                 aaTCGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTaa         >  1:3043964/1‑91 (MQ=255)

                                                                                  aTCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAATGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAAt        >  1:3111185/1‑91 (MQ=255)

                                                                                  aTCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCGAATTAAt        >  1:2076963/1‑91 (MQ=255)

                                                                                  aTCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAAt        <  1:3524936/91‑1 (MQ=255)

                                                                                   tCGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTCCTAACAAACAACTTTTCAAATTAATg       >  1:819134/1‑91 (MQ=255)

                                                                                    cGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATgg      >  1:1713730/1‑91 (MQ=255)

                                                                                    cGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATgg      >  1:1713871/1‑91 (MQ=255)

                                                                                     gTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGa     >  1:1971010/1‑91 (MQ=255)

                                                                                     gTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGa     <  1:3694432/91‑1 (MQ=255)

                                                                                       aaCTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAg   >  1:1274410/1‑91 (MQ=255)

                                                                                        aCTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGt  <  1:1914949/91‑1 (MQ=255)

AATTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGT  >  USA300TCH1516_ALE/2497302‑2497480

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