Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,129,380 C→G 7.1% intergenic (+213/+25) hlb_2 → / ← USA300HOU_RS10870 Phospholipase C/putative leukocidin‑like protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,129,3800CG 7.1% 966.9 / 41.5 281intergenic (+213/+25)hlb_2/USA300HOU_RS10870Phospholipase C/putative leukocidin‑like protein 1
Reads supporting (aligned to +/- strand):  ref base C (133/127);  new base G (5/15);  total (139/142)
Fisher's exact test for biased strand distribution p-value = 3.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 288 total aligned reads displayed.

GTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTC  >  USA300TCH1516_ALE/2129290‑2129470
                                                                                          |                                                                                          
gTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTc                                                                                            <  2:2047600/91‑1 (MQ=255)
  cTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTa                                                                                          <  1:1347522/91‑1 (MQ=255)
    ttACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATg                                                                                        >  1:863434/1‑91 (MQ=255)
    ttACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATg                                                                                        >  1:3517360/1‑91 (MQ=255)
        tctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGc                                                                                    >  2:128948/1‑91 (MQ=255)
        gctACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGc                                                                                    <  2:664314/90‑1 (MQ=255)
         ctACACGCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCa                                                                                   >  1:3069809/1‑91 (MQ=255)
          tACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCCATGTAGCAg                                                                                  >  2:825787/1‑91 (MQ=255)
            cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                                <  2:3708042/91‑1 (MQ=255)
            cacTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGc                                                                                >  2:3367217/1‑91 (MQ=255)
             acTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCa                                                                               <  1:305667/91‑1 (MQ=255)
                cAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACt                                                                            <  1:2612728/91‑1 (MQ=255)
                cAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACt                                                                            <  1:3496187/91‑1 (MQ=255)
                 agaCAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCGATGTAGCAGGCAACtt                                                                           <  2:579948/89‑1 (MQ=255)
                 aaaCAAGATTCATACACTGCACGTTATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACtt                                                                           <  2:3681269/91‑1 (MQ=255)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTtta                                                                        <  1:959250/91‑1 (MQ=17)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTtta                                                                        <  1:1946037/91‑1 (MQ=18)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTtta                                                                        <  2:3194058/91‑1 (MQ=18)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTtta                                                                        >  1:2181313/1‑91 (MQ=255)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTtta                                                                        >  2:1254570/1‑91 (MQ=255)
                    cAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTtta                                                                        <  1:494761/91‑1 (MQ=255)
                     aaGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttat                                                                       >  2:2809317/1‑91 (MQ=255)
                      aGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttatt                                                                      >  1:808905/1‑91 (MQ=255)
                       gATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTttatta                                                                     >  1:488175/1‑91 (MQ=255)
                        aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTGCCTGCTACATAGAATGTAGTAGGTAACTTttattta                                                                    <  1:104816/91‑3 (MQ=17)
                        aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTAc                                                                    >  1:2075530/1‑91 (MQ=255)
                        aTTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTAc                                                                    >  2:3131723/1‑91 (MQ=255)
                          tCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACtt                                                                  >  1:3370192/1‑91 (MQ=255)
                          gcATACGCTGCACGTCGTAATATATCTATCTATTGAAATAGAATTAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTGCtt                                                                  <  2:1708713/90‑1 (MQ=255)
                           cATACACTGCCCGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                 <  2:223856/91‑1 (MQ=255)
                           cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                 >  2:2162688/1‑91 (MQ=255)
                           cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                 <  1:1280170/91‑1 (MQ=255)
                           cATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTa                                                                 <  2:2080087/91‑1 (MQ=255)
                            aTACACTGGTCGTCATAATAAATCTATCTATTCAACTATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTGTATTACTTAt                                                                <  2:2507239/91‑1 (MQ=255)
                            aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                                >  1:973745/1‑91 (MQ=255)
                            aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                                <  1:1764193/91‑1 (MQ=255)
                            aTACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAt                                                                <  2:2644929/91‑1 (MQ=255)
                              acacTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTAttt                                                              <  2:1112715/91‑1 (MQ=255)
                                 cTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCtt                                                           <  1:3577827/91‑1 (MQ=255)
                                 cTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCtt                                                           >  1:234957/1‑91 (MQ=255)
                                  tGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCttt                                                          <  1:93958/91‑1 (MQ=255)
                                  tGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCttt                                                          >  1:1390071/1‑91 (MQ=255)
                                    cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                                        <  2:900542/91‑1 (MQ=255)
                                    cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                                        <  1:1712957/91‑1 (MQ=255)
                                    cACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtc                                                        >  2:3529288/1‑91 (MQ=255)
                                       gTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcatt                                                     >  1:721085/1‑91 (MQ=255)
                                          ataataAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatc                                                  <  2:18470/91‑1 (MQ=255)
                                           taataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatca                                                 >  1:1465004/1‑91 (MQ=255)
                                           taataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatca                                                 <  2:1394273/91‑1 (MQ=255)
                                            aataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcat                                                <  2:2364810/91‑1 (MQ=255)
                                            aataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcat                                                >  1:3872239/1‑91 (MQ=255)
                                            aataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTacattatcat                                                >  1:3792180/1‑91 (MQ=255)
                                             ataaATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatt                                               >  2:1512833/1‑91 (MQ=255)
                                             ataaATCTATATATTAAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTtcattatcatt                                               <  1:2768130/91‑1 (MQ=255)
                                               aaaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTaa                                             <  2:2916428/91‑1 (MQ=255)
                                               aaaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTaa                                             >  1:1879843/1‑91 (MQ=255)
                                                aaTCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAg                                            >  1:3087101/1‑91 (MQ=255)
                                                 atctatctATTCTAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACCTTTATTACTTATTTCTTTTCATTATCATTAAGt                                           <  2:3292364/91‑1 (MQ=255)
                                                 atctatctATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGt                                           <  2:1677655/91‑1 (MQ=255)
                                                   ctatctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTAc                                         >  2:3896167/1‑91 (MQ=255)
                                                   ctatctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATGATCATTAAGTAc                                         <  1:3027516/91‑1 (MQ=255)
                                                   ctatctatTCAAATATAAATAAAAGTTACCTACGACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTAc                                         <  1:3917874/91‑1 (MQ=255)
                                                    tatctatTCAAATATAAATAAAAGTGACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACt                                        <  1:3676943/91‑1 (MQ=255)
                                                     atctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtt                                       >  1:1918300/1‑91 (MQ=255)
                                                       ctatTCACATAGAAATGAAAGTTACCTACTACAGTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtttt                                     <  2:2433771/91‑1 (MQ=255)
                                                       ctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtttt                                     >  1:288295/1‑91 (MQ=255)
                                                       ctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtttt                                     <  1:3703108/91‑1 (MQ=255)
                                                       ctatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACtttt                                     <  1:3613534/91‑1 (MQ=255)
                                                       ctatTCAAATATAAAAGAACGTTACCTACTACGTGCTATGTAGCAGGCCACTTTTATTACTTATTTCTTTTCATTATCAGTAAGTACtttt                                     <  2:1100531/91‑1 (MQ=255)
                                                        tatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTa                                    <  1:2962646/91‑1 (MQ=255)
                                                        tatTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTa                                    >  1:681460/1‑91 (MQ=255)
                                                        tatTCAAATATAAATAAAAGTTACATACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTa                                    <  2:3685667/91‑1 (MQ=255)
                                                         atTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTAc                                   <  2:3486946/91‑1 (MQ=255)
                                                         atTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTAc                                   <  2:3472694/91‑1 (MQ=255)
                                                          ttCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACa                                  <  2:1602298/91‑1 (MQ=255)
                                                          ttCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACa                                  >  1:1527731/1‑91 (MQ=255)
                                                           tCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACaa                                 <  1:5530/91‑1 (MQ=255)
                                                           tCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACaa                                 >  1:3294593/1‑91 (MQ=255)
                                                            cAAATATAAATAAAAGTTACCTACTACATTCCCTGTAGCAGGCAACTTTTATTACTTATTTCTTTTCGCTATCATTAAGCACTTTTACaaa                                >  2:3022224/1‑91 (MQ=255)
                                                                tataAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTc                            >  1:3470191/1‑91 (MQ=255)
                                                                 ataaataaAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTca                           <  2:242538/91‑1 (MQ=255)
                                                                  taaataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTcac                          <  1:1618282/91‑1 (MQ=255)
                                                                   aaataaaAGTTACCTCCTACATTCTATGTAGCAGGCAACTTTTATTACTTATTGCTTATCATTATCATGAAGTGCTTTTATAAACTTcacc                         >  2:3541576/1‑90 (MQ=255)
                                                                    aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt                        <  1:2127058/91‑1 (MQ=255)
                                                                    aataaaAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACAt                        >  2:3608818/1‑91 (MQ=255)
                                                                        aaaGTTGCCTGCTACATAGAATGTAGTAGGTAACTTTTATTACTTATTTCTTTTCATTACCATTAAGTACTTTTACAAACTTCCCATTAtg                    >  2:1343929/1‑91 (MQ=18)
                                                                        aaaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtg                    >  1:1813707/1‑91 (MQ=255)
                                                                         aaGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgt                   >  2:1767609/1‑91 (MQ=255)
                                                                          aGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTAtgtg                  >  1:2473740/1‑91 (MQ=255)
                                                                            ttACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTACTTCTTTTCATTTTCATTAAGTACTTTTACAAACTTCACATTATGTGTc                >  2:50419/1‑91 (MQ=255)
                                                                             tACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCt               >  1:2847437/1‑91 (MQ=255)
                                                                             tACCTACTACATTCTATGTAGCAGGCAACTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt              <  2:1407892/91‑1 (MQ=255)
                                                                              aCCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCtt              >  1:3212661/1‑91 (MQ=255)
                                                                               ccTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTc             >  2:362955/1‑91 (MQ=255)
                                                                                  actacATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCaa          <  1:3828556/91‑1 (MQ=255)
                                                                                   ctacATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAAt         <  1:2596006/91‑1 (MQ=255)
                                                                                    tacATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATc        <  2:3403589/91‑1 (MQ=255)
                                                                                       aTTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCACTCAAc     >  2:2819935/1‑91 (MQ=255)
                                                                                          cTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTc  >  1:1010619/1‑91 (MQ=255)
                                                                                          cTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTc  <  2:3375303/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTCTTTACTCTACACTCAAACAAGATTCATACACTGCACGTCATAATAAATCTATCTATTCAAATATAAATAAAAGTTACCTACTACATTCTATGTAGCAGGCAACTTTTATTACTTATTTCTTTTCATTATCATTAAGTACTTTTACAAACTTCACATTATGTGTCTTCCAATCAACTTC  >  USA300TCH1516_ALE/2129290‑2129470

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: