BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,497,391	C→G	14.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	14.1%	770.8 / 96.3	261	N147K (AAC→AAG)	USA300HOU_RS12865	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (112/112); new base G (12/25); total (124/137)
Fisher's exact test for biased strand distribution p-value = 5.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 264 total aligned reads displayed.

AAATTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGTT  >  USA300TCH1516_ALE/2497301‑2497481

 aaTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACt                                                                                           >  1:2604424/1‑91 (MQ=255)

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

    ttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTAc                                                                                        <  1:291169/91‑1 (MQ=255)

    ttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTAc                                                                                        <  2:1292902/91‑1 (MQ=255)

     tACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACg                                                                                       <  1:3646351/91‑1 (MQ=255)

     tACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACg                                                                                       <  1:1921788/91‑1 (MQ=255)

      aCCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGa                                                                                      <  1:2837973/91‑1 (MQ=255)

      aCCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGa                                                                                      <  2:1118636/91‑1 (MQ=255)

       ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAt                                                                                     >  1:565011/1‑91 (MQ=255)

       ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAt                                                                                     >  1:2469595/1‑91 (MQ=255)

        cTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGAtt                                                                                    <  2:724185/91‑1 (MQ=255)

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

          aaGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGATTat                                                                                  <  1:354807/91‑2 (MQ=255)

          aaGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGATTat                                                                                  <  1:3112812/91‑2 (MQ=255)

          aaGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACGATTat                                                                                  <  2:1007031/91‑2 (MQ=255)

           aGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACt                                                                                 <  1:2026249/91‑1 (MQ=255)

            gCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACtt                                                                                <  2:199416/91‑1 (MQ=255)

             cGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTta                                                                               <  1:1523864/91‑1 (MQ=255)

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

               tGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaat                                                                             >  2:2175197/1‑91 (MQ=255)

               tGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaat                                                                             >  1:2813770/1‑91 (MQ=255)

                   atTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATcaaa                                                                         >  1:3128786/1‑91 (MQ=255)

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

                   atTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATaaaa                                                                         <  1:2256788/91‑1 (MQ=255)

                     tAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGt                                                                       >  2:2282369/1‑91 (MQ=255)

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

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

                       aTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAg                                                                     >  2:2059930/1‑91 (MQ=255)

                       aTTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAATAg                                                                     >  1:3058186/1‑91 (MQ=255)

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

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

                        tttAAAATTTCACAGCAGATTTGCATGAGAAATTTTAACTATTTTTGCCGACCATTATGATGGTAACTTAGGATTAGTTAATAAAGGTAGa                                                                    <  2:285809/91‑1 (MQ=255)

                         ttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGaa                                                                   <  2:621929/91‑1 (MQ=255)

                         ttAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGaa                                                                   <  2:2186073/91‑1 (MQ=255)

                          tAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAAt                                                                  <  1:3253350/91‑1 (MQ=255)

                            aaaTTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTc                                                                >  2:1967625/1‑91 (MQ=255)

                             aaTTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCg                                                               <  2:941980/91‑1 (MQ=255)

                              aTTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACGTACGATTACTTAATAAAAGTAGAATTCGt                                                              <  2:662036/91‑1 (MQ=255)

                               tttCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTa                                                             <  2:413836/91‑1 (MQ=255)

                               tttCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTa                                                             <  1:59811/91‑1 (MQ=255)

                                ttCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTAt                                                            <  2:489548/91‑1 (MQ=255)

                                 tCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTAtt                                                           >  1:2015128/1‑91 (MQ=255)

                                 tCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTAtt                                                           >  2:276764/1‑91 (MQ=255)

                                 tCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTAtt                                                           >  2:2720342/1‑91 (MQ=255)

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

                                     agcagATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTgaga                                                       >  1:3781261/1‑91 (MQ=255)

                                      gcagATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaa                                                      <  1:1750913/91‑1 (MQ=255)

                                       cagATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaaa                                                     <  1:1435900/91‑1 (MQ=255)

                                            ttGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACt                                                >  2:1346389/1‑91 (MQ=255)

                                            ttGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACt                                                >  2:3357126/1‑91 (MQ=255)

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

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

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

                                                 ttAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTaaaaa                                           <  1:1017521/91‑1 (MQ=255)

                                                  tAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAAAATTCGTATTGAGAAAGAACTaaaaaa                                          <  2:2262208/91‑1 (MQ=255)

                                                   aGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAAt                                         >  2:2607217/1‑91 (MQ=255)

                                                     aaaTTTTAACTACTATTGCCGACCATTATAATCGTAACTTAAGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATaa                                       <  2:2377189/91‑1 (MQ=255)

                                                      aaTTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAAc                                      >  2:2701500/1‑91 (MQ=255)

                                                         tttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACAtt                                   >  1:2984315/1‑91 (MQ=255)

                                                         tttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACAtt                                   >  2:1809103/1‑91 (MQ=255)

                                                         tttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACAtt                                   <  2:1863745/91‑1 (MQ=255)

                                                            aaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACt                                <  1:3088104/91‑1 (MQ=255)

                                                             actactATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                               >  1:1983780/1‑91 (MQ=255)

                                                             actactATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                               >  2:370837/1‑91 (MQ=255)

                                                             actactATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                               >  2:44431/1‑91 (MQ=255)

                                                             actactATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTa                               >  2:1691250/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

                                                                 ctaTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaacaa                           >  1:1292094/1‑91 (MQ=255)

                                                                      gCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAaacaac                      >  2:1235340/1‑91 (MQ=255)

                                                                      gCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAaacaac                      >  2:3682928/1‑91 (MQ=255)

                                                                      gCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAaacaac                      >  2:3497239/1‑91 (MQ=255)

                                                                         gACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACttt                   >  2:2341811/1‑91 (MQ=255)

                                                                          aCCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACtttt                  >  2:2932681/1‑91 (MQ=255)

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

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

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

                                                                           ccATTATAATCATAACTTACGATCACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACGAAAAAATAACATTACTAACAAACAACTTTTc                 <  2:2963237/91‑1 (MQ=255)

                                                                             aTTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCaa               >  1:2594769/1‑91 (MQ=255)

                                                                              aagtaATCGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCaaa              >  1:720621/4‑91 (MQ=255)

                                                                               tataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAt             <  1:805681/91‑1 (MQ=255)

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

                                                                                 taATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTa           <  2:3466314/91‑1 (MQ=255)

                                                                                  aaTCGTAAGTTACGATTACTTAATAACAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTaa          >  1:1048365/1‑91 (MQ=255)

                                                                                  aaTCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTaa          >  2:1927829/1‑91 (MQ=255)

                                                                                  aaTCGTAACTTACGACTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTaa          >  1:3307871/1‑91 (MQ=255)

                                                                                   aTCGTAAGTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAAt         >  2:1980348/1‑91 (MQ=255)

                                                                                   aTCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGagaa                                                      >  1:2627141/1‑46 (MQ=255)

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

                                                                                     cGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATgg       <  2:3073014/91‑1 (MQ=255)

                                                                                      gTAAGTTACGATTACTTAATAAAAGCAGAATTCGTATTGAGAAAGAACCAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGa      >  2:90972/1‑91 (MQ=255)

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

                                                                                         aCTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGt   >  1:1615519/1‑91 (MQ=255)

                                                                                         aCTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGt   >  1:2167743/1‑91 (MQ=255)

AAATTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGTT  >  USA300TCH1516_ALE/2497301‑2497481

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