BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	2,497,391	C→G	11.2%	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	11.2%	812.2 / 73.2	259	N147K (AAC→AAG)	USA300HOU_RS12865	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (110/119); new base G (7/22); total (118/141)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

Only 100 of 265 total aligned reads displayed.

AAATTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGTT  >  USA300TCH1516_ALE/2497301‑2497481

aaaTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAc                                                                                            <  2:933943/91‑1 (MQ=255)

aaaTTTACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAc                                                                                            >  1:2980938/1‑91 (MQ=255)

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

   tttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTa                                                                                         <  1:1133079/91‑1 (MQ=255)

   tttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTa                                                                                         >  1:2576908/1‑91 (MQ=255)

   tttACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTa                                                                                         >  1:2577271/1‑91 (MQ=255)

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

     tACCTAAGCGTGATATTGATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAAGTTACg                                                                                       <  2:3354340/91‑1 (MQ=255)

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

     tACCTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAAACGTAAGTTACg                                                                                       <  2:2650827/91‑1 (MQ=255)

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

       ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAt                                                                                     <  1:2317935/91‑1 (MQ=255)

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

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

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

       ccTAAGCGTGATATTAATTTAAAATTTCACAGCAGATATGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGAt                                                                                     >  1:3345454/1‑91 (MQ=255)

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

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

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

          aaGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTAc                                                                                  <  1:868161/91‑1 (MQ=255)

          aaGCGTGATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTAc                                                                                  <  1:331284/91‑1 (MQ=255)

                gATATTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaata                                                                            <  2:2222349/91‑1 (MQ=255)

                 atatTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTtaataa                                                                           <  1:1006085/91‑1 (MQ=255)

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

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

                  tatTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATaaa                                                                          <  2:2044132/91‑1 (MQ=255)

                  tatTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATaaa                                                                          >  2:3620003/1‑91 (MQ=255)

                  tatTAATTTAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATaaa                                                                          >  2:3722627/1‑91 (MQ=255)

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

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

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

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

                        tttAAAATATCACAGCAGATTTGCATTAGAAATTTTAACTACTATAGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGa                                                                    >  2:1735249/1‑91 (MQ=255)

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

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

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

                          tAAAATTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTACCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAAt                                                                  <  2:3804723/91‑1 (MQ=255)

                              aTTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGt                                                              <  1:2276012/91‑1 (MQ=255)

                              aTTTCACAGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGt                                                              <  1:229294/91‑1 (MQ=255)

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

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

                                  cacaGCAGATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTg                                                          <  1:3148647/91‑1 (MQ=255)

                                  cacaGCAGATTAGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTg                                                          >  2:3353/1‑91 (MQ=255)

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

                                     agcagATTTGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTg                                                          >  1:2923484/1‑88 (MQ=255)

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

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

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

                                           tttGCATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAAc                                                 >  2:1617968/1‑91 (MQ=255)

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

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

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

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

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

                                               cATTAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTaaa                                             >  2:2801437/1‑91 (MQ=255)

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

                                                 ttAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTaaaaa                                           >  2:2418341/1‑91 (MQ=255)

                                                  tAGAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTaaaaaa                                          <  2:1107645/91‑1 (MQ=255)

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

                                                    gAAATTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATa                                        >  1:2149521/1‑91 (MQ=255)

                                                     aaaTTTTAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATaa                                       <  1:2163036/91‑1 (MQ=255)

                                                      aaTTTTAACTACTAGTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACAAAAAAATAAc                                      <  2:1001399/91‑1 (MQ=255)

                                                       aTTTTAACTACTATTGCCGACCATGATAATCGTAACTTAAGATTAATTAATAAAAGTAGAATTCGTATTGAGAAAAAAATAAAAAATAACa                                     <  2:2821589/91‑1 (MQ=255)

                                                        ttttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACAt                                    >  1:2381821/1‑91 (MQ=255)

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

                                                          ttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTa                                  <  2:3223837/91‑1 (MQ=255)

                                                          ttAACTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTa                                  >  2:1236839/1‑91 (MQ=255)

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

                                                            aaCTACTATTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACt                                >  1:669785/1‑91 (MQ=255)

                                                              ctactcTTGCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaa                              >  2:2464987/1‑91 (MQ=255)

                                                                actaTTGCCGACCATTATAATCGTAACTTACGAATACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaaca                            <  1:3501641/91‑1 (MQ=255)

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

                                                                   aTTGCCGTCCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTaaccaac                         <  2:3508264/91‑1 (MQ=255)

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

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

                                                                      gCCGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTCGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAAcaaccaac                      >  2:2524698/1‑91 (MQ=255)

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

                                                                        cGACCATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACtt                    <  1:2625636/91‑1 (MQ=255)

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

                                                                            cATTATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCa                <  2:1504821/91‑1 (MQ=255)

                                                                              ttATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCaaa              <  2:489818/91‑1 (MQ=255)

                                                                              ttATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCaaa              <  1:3000285/91‑1 (MQ=255)

                                                                              ttATAATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCaaa              <  1:2270631/91‑1 (MQ=255)

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

                                                                                gtaATCGTAAGTTACGATTACTTACTAAAAGTAGCATTCGCATTGCGAAAGAACTAAAAAATAACACTACTAACATACACCTTTTCAAAtt            >  2:3673247/2‑91 (MQ=255)

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

                                                                                ataATCGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAAAACTAAAAAATAACATTACTAACAAACAACTTTTCAAAtt            <  2:1358308/91‑1 (MQ=255)

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

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

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

                                                                                     cGTAAGTTACGATTACTTAATAAACGTAGAATTCGTATTGCGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATgg       >  2:44493/1‑91 (MQ=255)

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

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

                                                                                     cGTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATgg       >  2:2975534/1‑91 (MQ=255)

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

                                                                                      gTAACTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGa      >  1:3782984/1‑91 (MQ=255)

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

                                                                                       tAAGTTACGATTACTTAATAAAAGTAGAATTCGTAATGAGAAAGAACTAAAAAACAACATTACTAACAAACAACTTTTCAAATAAATGGaa     >  2:1843319/1‑91 (MQ=255)

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

                                                                                         aCTTACGATTACTTAATAAAAGTAGAATTCGTATTGAGAAAGAACTAAAAAATAACATTACTAACAAACAACTTTTCAAATTAATGGAAGt   >  2:1994285/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` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`