BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,219,574	T→A	5.1%	Q299L (CAG→CTG)	USA300HOU_RS06070 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,219,574	0	T	A	5.1%	790.8 / 23.2	229	Q299L (CAG→CTG)	USA300HOU_RS06070	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (113/104); new base A (2/10); total (115/114)
Fisher's exact test for biased strand distribution p-value = 1.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

TCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTA  >  USA300TCH1516_ALE/1219484‑1219664

tCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAActgct                                                                                            >  1:2217671/1‑91 (MQ=255)

tCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAActgct                                                                                            >  2:2520419/1‑91 (MQ=255)

  aCCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGt                                                                                          <  1:4023/91‑1 (MQ=255)

  aCCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGt                                                                                          <  1:1578153/91‑1 (MQ=255)

  aCCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGt                                                                                          <  1:1914874/91‑1 (MQ=255)

   ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                         >  1:453306/1‑91 (MQ=255)

   ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                         >  2:167416/1‑91 (MQ=255)

   ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                         >  2:2972394/1‑91 (MQ=255)

    cATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTg                                                                                        <  1:2591953/91‑1 (MQ=255)

      tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTgta                                                                                      <  1:426639/91‑3 (MQ=255)

      tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGaa                                                                                      <  1:326865/91‑1 (MQ=255)

      tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGaa                                                                                      <  2:3608295/91‑1 (MQ=255)

       ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCATTCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                     >  2:2728037/1‑91 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACa                                                                                    <  1:509101/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACa                                                                                    <  1:2019612/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACa                                                                                    <  1:1419677/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACa                                                                                    <  1:1398741/91‑1 (MQ=255)

        taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACa                                                                                    >  1:292782/1‑91 (MQ=255)

         aaCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaa                                                                                   <  1:3255164/91‑1 (MQ=39)

         aaCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaa                                                                                   >  1:1682064/1‑91 (MQ=255)

          aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaaa                                                                                  <  2:1332152/91‑1 (MQ=255)

           cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGTACAAAt                                                                                 >  2:825420/1‑91 (MQ=255)

           cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAAt                                                                                 >  2:2178484/1‑91 (MQ=255)

           cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAAt                                                                                 >  1:1714051/1‑91 (MQ=255)

             gTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCg                                                                               >  2:3094751/1‑91 (MQ=255)

              ttCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGa                                                                              <  1:2473997/91‑1 (MQ=255)

              ttCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGa                                                                              >  2:314636/1‑91 (MQ=255)

               tCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGaa                                                                             >  1:1316715/1‑91 (MQ=255)

               tCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGaa                                                                             >  1:2228881/1‑91 (MQ=255)

                cAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAAc                                                                            <  1:2902785/91‑1 (MQ=255)

                cAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAAc                                                                            <  2:535850/91‑1 (MQ=255)

                cAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAAc                                                                            <  1:3481525/91‑1 (MQ=255)

                  gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                          <  2:1112923/91‑1 (MQ=255)

                  gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                          >  2:1243718/1‑91 (MQ=255)

                  gTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTa                                                                          <  2:3335011/91‑1 (MQ=255)

                   tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACtca                                                                         >  1:902428/1‑89 (MQ=255)

                       tttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCa                                                                     <  1:688359/91‑1 (MQ=255)

                       tttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCa                                                                     <  1:1274808/91‑1 (MQ=255)

                        ttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCAt                                                                    >  2:2917689/1‑91 (MQ=255)

                           tttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCg                                                                 >  1:1114681/1‑91 (MQ=255)

                           tttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCg                                                                 <  2:1013901/91‑1 (MQ=255)

                              ttAGTCTGCTCTTTTTCTTCTTTCAACTTCGCTAATTTACTTTGATATTTGTACAACTGCTGTTGAACAAATCGATCTAAATCATTCGCAc                                                              <  2:3861213/91‑1 (MQ=255)

                              ttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCAc                                                              >  2:2186376/1‑91 (MQ=255)

                              ttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCAc                                                              <  2:1563060/91‑1 (MQ=255)

                               tAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACg                                                             <  1:758858/91‑1 (MQ=255)

                                 gACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgtt                                                           >  1:1965818/1‑91 (MQ=255)

                                 gACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgtt                                                           >  1:3849165/1‑91 (MQ=255)

                                   cTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTGGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgttgt                                                         <  1:3140050/91‑1 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAAGTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGttt                                                       >  2:3483541/1‑91 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGttt                                                       >  1:3432649/1‑91 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGttt                                                       <  2:2373542/91‑1 (MQ=255)

                                     gCTCATATTCTTCAATCAACTTCGCTAATTTATCTTGATATTTGTACAACTGCTGTTGATCAAATCGAACTAAATCATTCGCACGTTGttt                                                       <  1:3512854/91‑1 (MQ=255)

                                      cTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTa                                                      <  2:3184450/91‑1 (MQ=255)

                                        cATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAc                                                    >  1:2107246/1‑91 (MQ=255)

                                           atTCTTCAATCGACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAAcgcg                                                 <  1:3407216/91‑1 (MQ=255)

                                            ttcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGt                                                <  1:729586/91‑1 (MQ=255)

                                            ttcttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGt                                                <  2:3778944/91‑1 (MQ=255)

                                                  aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGt                                          <  1:1875579/91‑1 (MQ=255)

                                                   atcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGtt                                         >  1:2419193/1‑91 (MQ=255)

                                                   atcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGtt                                         >  2:1887830/1‑91 (MQ=255)

                                                    tcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTc                                        >  1:2790401/1‑91 (MQ=255)

                                                     caaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCg                                       >  2:176111/1‑91 (MQ=255)

                                                        cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                    >  2:1131793/1‑91 (MQ=255)

                                                        cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                    >  2:1847300/1‑91 (MQ=255)

                                                           cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATTATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                                 >  2:1057749/1‑91 (MQ=255)

                                                            gCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgcg                                <  1:560375/91‑1 (MQ=255)

                                                            gCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgcg                                >  2:597423/1‑91 (MQ=255)

                                                             cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                               <  1:3440817/91‑1 (MQ=255)

                                                               aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                             >  1:2166630/1‑91 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            >  2:3809063/1‑91 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            <  1:2523271/91‑1 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            >  1:3231552/1‑91 (MQ=255)

                                                                atttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatc                            >  1:914723/1‑91 (MQ=255)

                                                                  ttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatcat                          <  2:491924/91‑1 (MQ=255)

                                                                  ttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatcat                          >  1:2841868/1‑91 (MQ=255)

                                                                   tatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATa                         <  2:585233/91‑1 (MQ=255)

                                                                    atttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaa                        <  1:1241349/91‑1 (MQ=255)

                                                                     ttttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaa                       <  1:1063256/91‑1 (MQ=255)

                                                                       ttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaaa                     <  1:2480805/91‑1 (MQ=255)

                                                                       ttGATACTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaaa                     >  1:3077236/1‑91 (MQ=255)

                                                                         gATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACg                   >  1:1453817/1‑91 (MQ=255)

                                                                           tatTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGTGTTCACGTTCGCCACGCGCATCATAAAAACGAt                 <  1:3029701/91‑1 (MQ=255)

                                                                           tatTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGAt                 >  2:1614166/1‑91 (MQ=255)

                                                                             tttGTACAGCTGCTGTTGAACAAATTTAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCa               <  2:3428980/91‑1 (MQ=255)

                                                                              ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa              >  2:2327632/1‑91 (MQ=255)

                                                                               tGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAg             >  1:2953336/1‑91 (MQ=255)

                                                                                gTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc            <  2:1377052/91‑1 (MQ=255)

                                                                                 ttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCa           >  1:1467381/3‑91 (MQ=255)

                                                                                 tACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCa           <  1:2683881/91‑1 (MQ=255)

                                                                                    aaCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAAt        >  2:1468699/1‑91 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  2:2066479/1‑91 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  2:1222834/1‑91 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  2:2620748/1‑91 (MQ=255)

                                                                                     aCTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATc       >  2:3958203/1‑91 (MQ=255)

                                                                                      ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa      <  1:812046/91‑1 (MQ=255)

                                                                                       tgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCAt     <  1:374740/91‑1 (MQ=255)

                                                                                        gctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCAtt    <  1:538850/91‑1 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  2:415648/91‑1 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  >  2:2932355/1‑91 (MQ=255)

                                                                                          tgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTa  <  2:2571125/91‑1 (MQ=255)

TCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCATTTA  >  USA300TCH1516_ALE/1219484‑1219664

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