BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,219,571	T→A	7.6%	Q300L (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,571	0	T	A	7.6%	643.2 / 28.8	197	Q300L (CAG→CTG)	USA300HOU_RS06070	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (94/87); new base A (5/10); total (99/98)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 208 total aligned reads displayed.

AATTCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCA  >  USA300TCH1516_ALE/1219481‑1219660

aaTTCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAAct                                                                                           >  1:1700074/1‑91 (MQ=255)

aaTTCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAAct                                                                                           >  2:1658069/1‑91 (MQ=255)

  ttCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAActgc                                                                                         >  1:3169862/1‑91 (MQ=255)

      ccATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGtt                                                                                     <  1:2368779/91‑1 (MQ=255)

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

        atatAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGa                                                                                   >  2:1614319/1‑91 (MQ=255)

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

         tataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTgta                                                                                  <  2:498101/91‑3 (MQ=255)

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

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

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

          ataACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAAc                                                                                 <  1:1276517/91‑1 (MQ=255)

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

           taACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACa                                                                                <  1:3484712/91‑1 (MQ=255)

            aaCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaa                                                                               <  1:3728588/91‑1 (MQ=38)

             aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACaaa                                                                              <  1:306551/91‑1 (MQ=38)

             aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaaa                                                                              >  1:2716200/1‑91 (MQ=255)

             aCTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACaaa                                                                              >  1:1400598/1‑91 (MQ=255)

              cTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTTCAACAGCAGTTGTACAAAt                                                                             <  1:773377/91‑1 (MQ=38)

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

               tGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATc                                                                            <  1:1134049/91‑1 (MQ=255)

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

                    aGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACt                                                                       <  1:1368493/91‑1 (MQ=255)

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

                      tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTTATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                     >  2:3784966/1‑91 (MQ=255)

                      tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                     <  1:2497743/91‑1 (MQ=255)

                      tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                     >  2:2596546/1‑91 (MQ=255)

                      tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                     <  2:328272/91‑1 (MQ=255)

                      tATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTaa                                                                     >  2:141162/1‑91 (MQ=255)

                          tttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCa                                                                 >  2:1502007/1‑91 (MQ=255)

                           ttATTTTTAGACTGCTCATTTTCTTCAATCAACTTCGCTAGTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCAt                                                                <  2:663657/91‑1 (MQ=255)

                           ttATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCAt                                                                <  1:2577310/91‑1 (MQ=255)

                             aTTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTGCAACTGCTGTTGAACAAATCGAACTAAATCATTc                                                              >  1:529543/1‑91 (MQ=255)

                             aTTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTc                                                              >  1:3832590/1‑91 (MQ=255)

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

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

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

                               ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                            >  2:1450484/1‑91 (MQ=255)

                               ttttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGc                                                            >  2:2764319/1‑91 (MQ=255)

                                tttAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCa                                                           <  1:3466227/91‑1 (MQ=255)

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

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

                                   aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                        <  2:1734617/91‑1 (MQ=255)

                                   aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                        >  2:12908/1‑91 (MQ=255)

                                   aGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgt                                                        >  2:1039290/1‑91 (MQ=255)

                                    gACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgtt                                                       <  1:1159561/91‑1 (MQ=255)

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

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

                                     aCTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGCTATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgttg                                                      >  2:3607050/1‑91 (MQ=255)

                                       tGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACgttgtt                                                    <  1:3085844/91‑1 (MQ=255)

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

                                               ttcttcAATCAACTTCGCTAATTTATCTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGt                                            >  2:2211878/1‑91 (MQ=255)

                                                 cttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTc                                          <  2:2864374/91‑1 (MQ=255)

                                                 cttcAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTc                                          <  1:952231/91‑1 (MQ=255)

                                                     aatcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGt                                      >  2:1952122/1‑91 (MQ=255)

                                                       tcaaCTTCGCTAATTTATTTTGATATTTGTGCAGCTGCTGTTGAACAAATCGAACTAAATCATTCGCATGTTGTTTAACGCGTTCACGTTc                                    <  2:3425040/91‑1 (MQ=255)

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

                                                       tcaaCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACCAAATCATTCGCACGTTGTTTAACGCGTTCACGTTc                                    >  2:213513/1‑91 (MQ=255)

                                                          aCTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGcc                                 <  1:1330222/91‑1 (MQ=255)

                                                           cTTCGCTAATTTATTTTGTTATTTGTTCAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                >  2:3002299/1‑91 (MQ=255)

                                                           cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                <  2:897581/91‑1 (MQ=255)

                                                           cTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCa                                <  2:369975/91‑1 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             >  2:1824444/1‑91 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             >  1:2595456/1‑91 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             >  1:1524520/1‑91 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             <  1:2570214/91‑1 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             >  2:1098357/1‑91 (MQ=255)

                                                              cGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCAcgc                             <  1:2139045/91‑1 (MQ=255)

                                                                cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                           >  2:3480844/1‑91 (MQ=255)

                                                                cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                           >  2:3558455/1‑91 (MQ=255)

                                                                cTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACgcgc                           >  1:1037236/1‑91 (MQ=255)

                                                                  aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                         <  1:2090486/91‑1 (MQ=255)

                                                                  aaTTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcat                         <  1:888442/91‑1 (MQ=255)

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

                                                                    tttatttTGATATTTGTGCAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatca                       <  1:3883689/91‑1 (MQ=255)

                                                                    tttatttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGcatca                       <  2:2509596/91‑1 (MQ=255)

                                                                      tatttTTATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATa                     <  1:3501417/91‑1 (MQ=255)

                                                                       atttTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaa                    >  2:1726640/1‑91 (MQ=255)

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

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

                                                                         tttGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATaaaa                  >  2:402569/1‑91 (MQ=255)

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

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

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

                                                                                 ttGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCaa          <  1:225879/91‑1 (MQ=255)

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

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

                                                                                  tgttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGGTGTTTACCGCGTTCACGTTCGCCACGCGCATCCTAAAAACGCTCAAg         >  2:1596980/5‑91 (MQ=255)

                                                                                  tGTTCAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAg         >  2:714499/1‑91 (MQ=255)

                                                                                   gttcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc        >  1:486446/4‑91 (MQ=255)

                                                                                   gTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGc        >  1:2809030/1‑91 (MQ=255)

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

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

                                                                                   gTACAACTGCTGTTGAACAAATCGAACTAAACCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCATGCGCATCATAATAACGATCAAGc        >  2:2966118/1‑91 (MQ=255)

                                                                                     tcAACAGCAGTTGTACAAATCGAACTAAATCATTCGCCCGTTGTTTAACGCGTTCACGTTCGCCACGCGCCTCATAAAAACGCCCAAGCaa      >  2:54808/2‑91 (MQ=255)

                                                                                     aCAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaa      <  1:2074435/91‑1 (MQ=255)

                                                                                      cAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa     <  1:2191554/91‑1 (MQ=255)

                                                                                      cAACTGCTGTTAAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCaaa     <  2:3536587/91‑1 (MQ=255)

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

                                                                                         ctgctgTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCa  >  1:3380084/1‑91 (MQ=255)

AATTCACCATATAACTGTTCAGTATCTTTATTTTTAGACTGCTCATATTCTTCAATCAACTTCGCTAATTTATTTTGATATTTGTACAACTGCTGTTGAACAAATCGAACTAAATCATTCGCACGTTGTTTAACGCGTTCACGTTCGCCACGCGCATCATAAAAACGATCAAGCAAATCA  >  USA300TCH1516_ALE/1219481‑1219660

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