Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,128,680 T→A 5.8% 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 5.8% 815.0 / 24.4 240Y286N (TAT→AAT) ctaB2Protoheme IX farnesyltransferase 2
Reads supporting (aligned to +/- strand):  ref base T (111/115);  new base A (7/7);  total (118/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 258 total aligned reads displayed.

TTAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGTT  >  USA300TCH1516_ALE/1128590‑1128769
                                                                                          |                                                                                         
ttAGGATAGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAatt                                                                                           <  1:1937435/91‑1 (MQ=255)
  aGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattat                                                                                         <  2:735266/91‑1 (MQ=255)
  aGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattat                                                                                         <  1:3258940/91‑1 (MQ=255)
  aGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattat                                                                                         <  2:3339492/91‑1 (MQ=255)
   ggatggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAattatt                                                                                        >  2:277914/1‑91 (MQ=255)
    gatggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTAttt                                                                                       >  2:56212/1‑91 (MQ=255)
     atggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTa                                                                                      <  2:963844/91‑1 (MQ=255)
      tggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAg                                                                                     <  2:3081624/91‑1 (MQ=255)
       ggatTGCATTAGGTTTAACGACATTTAAGAAGAATTCAGAGCGAACAAAATGGGCGACACAAATGTTTATATATTCACTAAATAATTTAGt                                                                                    <  2:1816324/91‑1 (MQ=255)
       ggatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGt                                                                                    >  2:862863/1‑91 (MQ=255)
       ggatTGAATTAGGTTTAAGAAAATTTAAGAAAAATTCAGATTAAAAAAAATGGGAAACAGAAATGTTTATATATTTACTAAATTATTTAGt                                                                                    >  2:358874/1‑91 (MQ=255)
        gatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTg                                                                                   >  2:1651041/1‑91 (MQ=255)
        gatTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTg                                                                                   <  2:1119817/91‑1 (MQ=255)
          ttGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGAt                                                                                 <  1:77809/91‑1 (MQ=255)
          ttGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTAGTGaa                                                                                 <  2:2149995/91‑2 (MQ=255)
           tGCATTAGGTTTTACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATc                                                                                <  1:1664221/91‑1 (MQ=255)
           tGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAAGTATTTAGTGATc                                                                                <  1:2996422/91‑1 (MQ=255)
            gCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCt                                                                               <  2:3113638/91‑1 (MQ=255)
            gCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATAATTTAGTGaata                                                                               <  1:3845346/91‑4 (MQ=255)
             cATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtt                                                                              <  1:2388667/91‑1 (MQ=255)
             cATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtt                                                                              >  1:2338475/1‑91 (MQ=255)
               ttAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtttt                                                                            <  1:2738436/91‑1 (MQ=255)
               ttAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtttt                                                                            <  2:1810187/91‑1 (MQ=255)
               ttAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCtttt                                                                            >  2:3872149/1‑91 (MQ=255)
                 aGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTc                                                                          >  2:1084995/1‑91 (MQ=255)
                  ggTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCg                                                                         >  2:3531772/1‑91 (MQ=255)
                     ttAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCgtgt                                                                      <  2:1043765/91‑1 (MQ=255)
                     ttAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCgtgt                                                                      <  1:956443/91‑1 (MQ=255)
                      tAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGtt                                                                     >  1:2214166/1‑91 (MQ=255)
                       aacaacATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATAGTCACTAAATTATTTAGTGATCTTTTTCGTGTTa                                                                    >  1:3183342/1‑91 (MQ=255)
                          aacaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCt                                                                 <  2:3493212/91‑1 (MQ=255)
                          aacaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCt                                                                 <  1:2849260/91‑1 (MQ=255)
                          aacaTTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCt                                                                 <  2:1028381/91‑1 (MQ=255)
                            caTTTAAGAAAAATTCAGATCAAACTAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCtgt                                                               <  1:3621025/91‑1 (MQ=255)
                              tttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGa                                                             >  1:868311/1‑91 (MQ=255)
                              tttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGa                                                             <  2:2820195/91‑1 (MQ=255)
                               ttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAt                                                            <  2:1328724/91‑1 (MQ=255)
                               ttAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAt                                                            <  2:1615099/91‑1 (MQ=255)
                                tAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAtt                                                           >  2:2161812/1‑91 (MQ=255)
                                tAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGAtt                                                           <  1:433466/91‑1 (MQ=255)
                                tAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGCTTATATATTCACTAAATTATTTAGTCGTCTTTTTCGTGTTAGCTGTGAtt                                                           >  2:3602835/1‑91 (MQ=255)
                                 aaGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCTCTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTg                                                          <  2:3196609/91‑1 (MQ=255)
                                   gAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGtt                                                        <  2:1253317/91‑1 (MQ=255)
                                   gAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGtt                                                        >  2:26280/1‑91 (MQ=255)
                                   gAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGtt                                                        >  2:1479224/1‑91 (MQ=255)
                                    aaaaaTTCAGATCAAACAAAATGGGCAACACAAATGTATATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGttt                                                       >  2:933158/1‑91 (MQ=255)
                                     aaaaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTc                                                      >  2:2385966/1‑91 (MQ=255)
                                       aaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAt                                                    <  1:278734/91‑1 (MQ=255)
                                       aaTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAt                                                    <  1:278111/91‑1 (MQ=255)
                                        aTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAtt                                                   >  2:1763989/1‑91 (MQ=255)
                                        aTTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAtt                                                   >  1:1243488/1‑91 (MQ=255)
                                           cAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttact                                                >  2:1724733/1‑91 (MQ=255)
                                           cAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttact                                                <  2:1711166/91‑1 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGt                                                                                    >  1:1177069/1‑53 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              >  2:2604604/1‑91 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              <  2:1949155/91‑1 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              <  2:1860898/91‑1 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              >  1:751601/1‑91 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              <  1:1471963/91‑1 (MQ=255)
                                             gATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttactta                                              >  2:505407/1‑91 (MQ=255)
                                              aTCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCAttacttac                                             <  2:916993/91‑1 (MQ=255)
                                                  aacaaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTg                                         <  2:1279410/91‑1 (MQ=255)
                                                   acaaaATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGa                                        >  2:2077664/1‑91 (MQ=255)
                                                     aaaaTGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATc                                      >  1:1720795/1‑91 (MQ=255)
                                                     aaaaTGGGCAACACAAATGTTTATATATTCAATAAATTATTTAGTGATCTTTTTCTTGTTAGCTGTGCTTGTTTCATTAGTTACTTTGATc                                      >  2:2969280/1‑91 (MQ=255)
                                                      aaaTGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCt                                     >  1:3564343/1‑91 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                 >  2:1220617/1‑91 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                 >  1:409629/1‑91 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                 <  2:3806912/91‑1 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGat                                 <  1:3216525/91‑1 (MQ=255)
                                                          gggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                <  2:862060/91‑1 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGttt                                                       >  2:301425/1‑68 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                <  1:2130238/91‑1 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                >  2:1493815/1‑91 (MQ=255)
                                                           ggCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatt                                >  1:342528/1‑91 (MQ=255)
                                                            gCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatta                               >  2:2721629/1‑91 (MQ=255)
                                                            gCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatta                               >  1:3739053/1‑91 (MQ=255)
                                                            gCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatta                               >  1:2604104/1‑91 (MQ=255)
                                                            gCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGatta                               <  2:2511729/91‑1 (MQ=255)
                                                                cacaAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGattaatta                           >  1:534574/1‑91 (MQ=255)
                                                                       gTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAgg                    <  1:1035407/91‑1 (MQ=255)
                                                                        tttATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGa                   >  1:2834242/1‑91 (MQ=255)
                                                                          tatataTTCACTAAATTATTTAGTGATCTTTTTCTTTTTTTCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGaat                 >  2:1178882/1‑89 (MQ=255)
                                                                           atatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGa                <  1:3203368/91‑1 (MQ=255)
                                                                           atatatTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGa                <  1:2554160/91‑1 (MQ=255)
                                                                               gatCACTAAATAATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTATTTTGATCTAGATTAATTAAGTTAGGATGaaaaa            >  2:3651298/3‑91 (MQ=255)
                                                                                ttCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAat           <  1:1210430/91‑1 (MQ=255)
                                                                                 tCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAata          <  2:1303410/91‑1 (MQ=255)
                                                                                  cACTAAATTATTTAGTGATCTTTTTCGTGTTATCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAatat         >  2:3445660/1‑91 (MQ=255)
                                                                                  cACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAatat         <  1:1743202/91‑1 (MQ=255)
                                                                                   aCTAAATAATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATg        >  1:2377074/1‑91 (MQ=255)
                                                                                    cTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATgg       >  2:1416194/1‑91 (MQ=255)
                                                                                      aaaTTATTTAGTGATCTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCg    <  2:2541148/91‑1 (MQ=255)
                                                                                       aaTTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCg    <  1:889882/91‑1 (MQ=255)
                                                                                       aaTTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCg    <  1:1597274/91‑1 (MQ=255)
                                                                                       aaTTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCg    >  1:1358352/1‑91 (MQ=255)
                                                                                        attattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGt   <  1:3032592/91‑1 (MQ=255)
                                                                                         ttattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGtt  <  2:3729161/91‑1 (MQ=255)
                                                                                         ttattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGtt  <  2:324655/91‑1 (MQ=255)
                                                                                         ttattTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGtt  <  1:3089485/91‑1 (MQ=255)
                                                                                          |                                                                                         
TTAGGATGGATTGCATTAGGTTTAACAACATTTAAGAAAAATTCAGATCAAACAAAATGGGCAACACAAATGTTTATATATTCACTAAATTATTTAGTGATCTTTTTCGTGTTAGCTGTGATTGTTTCATTACTTACTTTGATCTAGATTAATTAAGTTAGGATGAAAAATATGGGCGTT  >  USA300TCH1516_ALE/1128590‑1128769

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: