Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,128,680 T→A 9.3% Y286N (TAT→AAT)  ctaB2 → Protoheme IX farnesyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,128,6800TA 9.3% 738.9 / 53.9 237Y286N (TAT→AAT) ctaB2Protoheme IX farnesyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (118/96);  new base A (14/8);  total (132/105)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 267 total aligned reads displayed.

TTAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGTTC  >  USA300TCH1516_ALE/1128590‑1128770
                                                                                          |                                                                                          
 tAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAatta                                                                                           <  2:966304/91‑1 (MQ=255)
 tAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAatta                                                                                           >  2:740540/1‑91 (MQ=255)
 tAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAatta                                                                                           <  1:2880422/91‑1 (MQ=255)
 tAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAatta                                                                                           <  2:3070117/91‑1 (MQ=255)
  aGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattat                                                                                          >  2:2546921/1‑91 (MQ=255)
  aGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattat                                                                                          <  1:3724810/91‑1 (MQ=255)
   ggatggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattatt                                                                                         >  1:1330231/1‑91 (MQ=255)
   ggatggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAataatt                                                                                         <  2:2898537/91‑1 (MQ=255)
    gatggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAAttt                                                                                        <  1:2886447/91‑1 (MQ=255)
     atggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTa                                                                                       >  1:703642/1‑91 (MQ=255)
     atggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTa                                                                                       <  1:1176846/91‑1 (MQ=255)
      tggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAg                                                                                      <  2:1817981/91‑1 (MQ=255)
      tggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAg                                                                                      >  2:484359/1‑91 (MQ=255)
      tggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAg                                                                                      >  2:62597/1‑91 (MQ=255)
      tggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTAg                                                                                      <  1:2392944/91‑1 (MQ=255)
       ggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGt                                                                                     >  1:3238424/1‑91 (MQ=255)
       ggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGt                                                                                     >  2:22787/1‑91 (MQ=255)
       ggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGt                                                                                     <  2:3487266/91‑1 (MQ=255)
        gatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTg                                                                                    >  1:533420/1‑91 (MQ=255)
         atTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGa                                                                                   >  2:1111297/1‑91 (MQ=255)
          ttGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGAt                                                                                  <  2:1999159/91‑1 (MQ=255)
          ttGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGAt                                                                                  <  2:1923197/91‑1 (MQ=255)
           tGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTAGTGaat                                                                                 <  1:1753548/91‑3 (MQ=255)
            gCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCt                                                                                <  1:861082/91‑1 (MQ=255)
             cATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtt                                                                               >  1:992736/1‑91 (MQ=255)
             cATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtt                                                                               >  1:3155813/1‑91 (MQ=255)
             cATTAGGCTTAACAAAATTTCAGAAAAATCTAGATCAAACACAATGGGCAACACAAATGTCTATATATTCACTAAATTACTTAGTGATCtt                                                                               >  2:1191376/1‑91 (MQ=255)
               ttAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtttt                                                                             <  1:2272393/91‑1 (MQ=255)
                  ggTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCg                                                                          >  2:355105/1‑91 (MQ=255)
                   gTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCgt                                                                         >  2:1965808/1‑91 (MQ=255)
                   gTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCgt                                                                         <  2:3809947/91‑1 (MQ=255)
                    tttAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCgtg                                                                        >  2:2387942/1‑91 (MQ=255)
                      tAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGtt                                                                      >  1:47848/1‑91 (MQ=255)
                      tAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGtt                                                                      >  2:506605/1‑91 (MQ=255)
                      tAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGtt                                                                      >  1:654026/1‑91 (MQ=255)
                         caacaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGc                                                                   <  1:2420516/91‑1 (MQ=255)
                           acaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTTTTAGCtt                                                                 >  2:2856443/1‑90 (MQ=255)
                           acaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCtg                                                                 <  1:75409/91‑1 (MQ=255)
                           acaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCtg                                                                 <  1:3468792/91‑1 (MQ=255)
                            caTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCtgt                                                                >  1:2105362/1‑91 (MQ=255)
                              tttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGa                                                              <  1:950969/91‑1 (MQ=255)
                               ttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTCGCTGTTAt                                                             >  2:1254986/1‑91 (MQ=255)
                               ttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCtgtggt                                                             >  1:2981464/1‑89 (MQ=255)
                                tAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAtt                                                            <  1:621755/91‑1 (MQ=255)
                                tAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAtt                                                            >  1:2680221/1‑91 (MQ=255)
                                 aaGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTg                                                           >  2:2582365/1‑91 (MQ=255)
                                  aGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGt                                                          <  2:2922362/91‑1 (MQ=255)
                                   gAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGtt                                                         >  1:545354/1‑91 (MQ=255)
                                    aaaaaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGttt                                                        >  2:330409/1‑91 (MQ=255)
                                     aaaaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTc                                                       <  1:3036825/91‑1 (MQ=255)
                                      aaaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCa                                                      <  1:2651366/91‑1 (MQ=255)
                                      aaaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCa                                                      <  2:1605371/91‑1 (MQ=255)
                                       aaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAt                                                     >  1:2568378/1‑91 (MQ=255)
                                       aaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAt                                                     <  1:609587/91‑1 (MQ=255)
                                         ttCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAtta                                                   >  1:1517227/1‑91 (MQ=255)
                                           cAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttact                                                 >  2:1853241/1‑91 (MQ=255)
                                           cAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttact                                                 <  2:1713544/91‑1 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                               >  2:2093241/1‑91 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                               >  1:819866/1‑91 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                               >  1:2793878/1‑91 (MQ=255)
                                              aTCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttacttac                                              <  1:1922201/91‑1 (MQ=255)
                                               tCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCTTTTTATCTGTGATTGTTGCATtacttcct                                             >  2:264618/1‑91 (MQ=255)
                                               tCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATtacttact                                             <  1:1381174/91‑1 (MQ=255)
                                               tCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATtacttact                                             >  1:1746291/1‑91 (MQ=255)
                                               tCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATtacttact                                             >  2:15002/1‑91 (MQ=255)
                                                caaacaaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTacttactt                                            <  2:2736707/91‑1 (MQ=255)
                                                caaacaaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTacttactt                                            <  2:2727167/91‑1 (MQ=255)
                                                    caaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGAt                                        <  1:788084/91‑1 (MQ=255)
                                                    caaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTACTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGAt                                        <  1:3135830/91‑1 (MQ=255)
                                                        aTGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAg                                    <  2:567245/91‑1 (MQ=255)
                                                        aTGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAg                                    <  1:2051521/91‑1 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                  >  1:2213930/1‑91 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                  >  2:268092/1‑91 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                 <  2:2768519/91‑1 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTGAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                 <  1:3333539/91‑1 (MQ=255)
                                                            gCAACACGAGTGTTTATATATTCACTCGATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTGACTTTGATCTAGatta                                <  2:3110289/91‑1 (MQ=255)
                                                            gCAACACAGATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatta                                <  2:428071/91‑1 (MQ=255)
                                                                cacaAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGattaatta                            <  1:3628963/91‑1 (MQ=255)
                                                                cacaAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGattaatta                            >  1:3487448/1‑91 (MQ=255)
                                                                    aaTGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGtt                        >  1:2479941/1‑91 (MQ=255)
                                                                      tGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAg                      >  2:1435460/1‑91 (MQ=255)
                                                                         ttATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTTGCTTTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGAt                   >  2:3582803/1‑91 (MQ=255)
                                                                         ttATATATTCACTAAATGATTTAGGGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGAt                   <  1:1403218/91‑1 (MQ=255)
                                                                           atatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGa                 <  2:1045250/91‑1 (MQ=255)
                                                                           atatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGa                 <  1:2167611/91‑1 (MQ=255)
                                                                           atatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGa                 >  1:2367191/1‑91 (MQ=255)
                                                                            tatatTCACTATATTATTTAGTGATCTTTTTTGTGTTAGCTGTGATTGTTTGATTACTTACTTTGAGCTAGATTAATTAAGTTAGGATGaa                <  2:3768370/91‑1 (MQ=255)
                                                                            tatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGaa                >  2:3782086/1‑91 (MQ=255)
                                                                            tatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGaa                >  2:2219285/1‑91 (MQ=255)
                                                                              tatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGaaaa              >  1:1230847/1‑91 (MQ=255)
                                                                               atTCACTAGATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGaaaaa             <  2:2492576/91‑1 (MQ=255)
                                                                               atTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGaaaaa             >  1:2882402/1‑91 (MQ=255)
                                                                                  cACTAAATAATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAatat          >  1:249194/1‑91 (MQ=255)
                                                                                  cACTAAATAATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAatat          >  2:1100229/1‑91 (MQ=255)
                                                                                   aCTAAATAATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATg         >  1:72862/1‑91 (MQ=255)
                                                                                     tAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATggg       <  2:1627610/91‑1 (MQ=255)
                                                                                        attattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGt    <  2:1333457/91‑1 (MQ=255)
                                                                                         ttattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGtt   >  1:2017148/1‑91 (MQ=255)
                                                                                         ttattTAGTGATCTGGTTCGTGTTAGCTGTGATTGTTTCATGACTTACTGTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGtt   <  2:3183843/91‑1 (MQ=255)
                                                                                          tattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAACTATGGGCGTTc  >  1:3575853/1‑91 (MQ=255)
                                                                                          |                                                                                          
TTAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGTTC  >  USA300TCH1516_ALE/1128590‑1128770

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: