BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	690,106	C→T	5.2%	128 (TAG→TAA)	USA300HOU_RS03325 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	690,106	0	C	T	5.2%	934.1 / 18.3	274	128 (TAG→TAA)	USA300HOU_RS03325	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (135/123); new base T (12/4); total (147/127)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATATT  >  USA300TCH1516_ALE/690016‑690196

gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                            <  2:2345236/91‑1 (MQ=255)

gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                            <  1:3199897/91‑1 (MQ=255)

  tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTa                                                                                          >  1:236923/1‑91 (MQ=255)

     cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAttt                                                                                       <  2:463261/91‑1 (MQ=38)

     cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                       >  1:2536131/1‑91 (MQ=255)

     cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                                       >  2:3002815/1‑91 (MQ=255)

      gTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTa                                                                                      >  2:647314/1‑91 (MQ=255)

         aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                                   <  1:1896517/91‑1 (MQ=255)

          ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                                  <  2:3709418/91‑1 (MQ=255)

            aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                                <  2:1507527/91‑1 (MQ=255)

            aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                                <  1:2733254/91‑1 (MQ=255)

            aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                                <  2:1220672/91‑1 (MQ=255)

               aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                             <  2:2003757/91‑1 (MQ=255)

               aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                             <  1:738028/91‑1 (MQ=255)

                aaTTATCTCTTAGACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttt                                                                            >  1:3036462/1‑91 (MQ=255)

                 aTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttt                                                                           <  1:2188858/91‑1 (MQ=255)

                  ttATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttt                                                                          <  1:1644079/91‑1 (MQ=255)

                   tATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttttt                                                                         <  1:81230/91‑1 (MQ=255)

                   tATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttttt                                                                         <  1:3127338/91‑1 (MQ=255)

                     tctcTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTa                                                                       >  1:1016683/1‑91 (MQ=255)

                      ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAGTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                      <  2:1005228/91‑1 (MQ=255)

                      ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                      <  1:370470/91‑1 (MQ=255)

                      ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                      <  2:1386291/91‑1 (MQ=255)

                       tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                     <  1:966283/91‑1 (MQ=255)

                       tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                     <  1:3306417/91‑1 (MQ=255)

                       tctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATc                                                                     <  2:1814534/91‑1 (MQ=255)

                          tAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                                  <  1:2666876/91‑1 (MQ=255)

                          tAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                                  <  2:710926/91‑1 (MQ=255)

                          tAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                                  >  2:3493027/1‑91 (MQ=255)

                          tAAACCTGCATACGGTGATTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATg                                                                  <  1:3389271/91‑1 (MQ=255)

                              ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                              <  1:1479729/91‑1 (MQ=255)

                                tGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTc                                                            >  2:2799798/1‑91 (MQ=255)

                                 gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                           >  1:3279612/1‑91 (MQ=255)

                                  cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                          <  2:2620501/91‑1 (MQ=255)

                                  cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                          <  2:1871386/91‑1 (MQ=255)

                                   aTACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAc                                                         >  1:3012633/1‑91 (MQ=255)

                                   aTACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAc                                                         <  1:483730/91‑1 (MQ=255)

                                    tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATCGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                        >  1:3339696/1‑91 (MQ=255)

                                    tACGGTGCTCTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATCTTCTAct                                                        >  1:1971425/1‑91 (MQ=255)

                                      cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                      >  2:3916253/1‑91 (MQ=255)

                                      cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                      >  2:3875977/1‑91 (MQ=255)

                                       ggTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtct                                                     <  1:3549177/91‑1 (MQ=255)

                                       ggTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtct                                                     <  1:1417184/91‑1 (MQ=255)

                                        gTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtctt                                                    <  1:1295368/91‑1 (MQ=255)

                                        gTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtctt                                                    >  1:1381971/1‑91 (MQ=255)

                                         tGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTGTTTATCATGATATTCTACttcttc                                                   <  2:3361706/91‑1 (MQ=255)

                                            ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtcatc                                                >  1:38290/1‑91 (MQ=255)

                                            ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTtcatc                                                <  1:2623600/91‑1 (MQ=255)

                                              ttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTcatcat                                              >  2:2156898/1‑91 (MQ=255)

                                              ttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTcatcat                                              >  1:2053666/1‑91 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           >  2:2605469/1‑91 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           <  1:1751850/91‑1 (MQ=255)

                                                 aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                           >  2:3068015/1‑91 (MQ=255)

                                                  aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                          >  1:1738054/1‑91 (MQ=255)

                                                  aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                          <  1:2910849/91‑1 (MQ=255)

                                                     cATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTaa                                       >  2:2761165/1‑91 (MQ=255)

                                                     cATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTaa                                       >  2:3222993/1‑91 (MQ=255)

                                                      aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                      >  1:3866910/1‑91 (MQ=255)

                                                       tGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAtt                                     >  2:60818/1‑91 (MQ=255)

                                                       tGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAtt                                     >  1:922353/1‑91 (MQ=255)

                                                        gTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAttt                                    >  2:2504044/1‑91 (MQ=255)

                                                         tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                                   <  1:2939536/91‑1 (MQ=255)

                                                          cATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTaa                                  >  1:3499199/1‑91 (MQ=255)

                                                           aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                                 >  1:1030914/1‑91 (MQ=255)

                                                             cccTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtt                               <  2:2195245/91‑1 (MQ=255)

                                                              ccTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACttt                              <  2:2708704/91‑1 (MQ=255)

                                                               cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                             <  1:574484/91‑1 (MQ=255)

                                                               cTCAAATTAAAATAACCCCTATAGGGGCTATTTAAATGTTATTTTTTTAACATGATATTCTACTTCTTCATCAAAACTTAATTTAACtttt                             <  2:2313807/91‑1 (MQ=255)

                                                                  cAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCtt                          >  2:174847/1‑91 (MQ=255)

                                                                    aTAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                        >  2:1351982/1‑91 (MQ=255)

                                                                       aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg                     <  1:2881457/91‑1 (MQ=255)

                                                                       aaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTg                     >  1:523182/1‑91 (MQ=255)

                                                                         aaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca                   >  1:134061/1‑91 (MQ=255)

                                                                         aaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATTATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca                   >  2:3273643/1‑91 (MQ=255)

                                                                         aaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca                   >  1:159301/1‑91 (MQ=255)

                                                                          aTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTCATCTAACTTTTCTCAACCGcac                  >  2:636136/1‑91 (MQ=255)

                                                                          aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcac                  >  2:1482052/1‑91 (MQ=255)

                                                                           tAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcaca                 >  1:2704856/1‑91 (MQ=255)

                                                                            agccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt                >  2:1075496/3‑91 (MQ=255)

                                                                            agccccTATTGGGGTTATTTAAATGTTA‑TTTTTTATCATGATATTCTACTTTTTCATTCTAACTTAATTTTAACTTTCTTCACTGCACAtt               >  2:2595483/3‑91 (MQ=255)

                                                                            aaCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt                >  1:794402/1‑91 (MQ=255)

                                                                             gccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt               >  1:2092180/2‑91 (MQ=255)

                                                                              ccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa              >  2:1426407/1‑91 (MQ=255)

                                                                               cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa             >  2:1593587/1‑91 (MQ=255)

                                                                                ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa            <  1:362887/91‑1 (MQ=255)

                                                                                 cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa           <  2:3834844/91‑1 (MQ=255)

                                                                                  aaTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaaa          >  2:1045223/1‑91 (MQ=255)

                                                                                   aTAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAAc         >  2:2817505/1‑91 (MQ=255)

                                                                                      ggggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACaaa      >  2:3221959/1‑91 (MQ=255)

                                                                                       gggTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAat     >  2:2054582/1‑91 (MQ=255)

                                                                                       gggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAat     <  2:815786/91‑1 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATTATGTTATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  2:1214372/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  2:3697935/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  2:3154229/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  2:2917275/1‑91 (MQ=255)

                                                                                        ggCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAata    >  1:1207779/1‑91 (MQ=255)

                                                                                         gCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAatat   <  1:435719/91‑1 (MQ=255)

                                                                                         gCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAatat   <  2:1626281/91‑1 (MQ=255)

                                                                                         gCTATTTAAATGTGATTTTTTTATCATGATATTCTACTTCTTCGTCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAatat   <  1:3954327/91‑1 (MQ=255)

                                                                                          cTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATAtt  <  1:2155949/91‑1 (MQ=255)

GTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAAATATT  >  USA300TCH1516_ALE/690016‑690196

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