Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 413,761 A→T 6.0% I166I (ATT→ATA metI ← Cystathionine gamma‑synthase/O‑acetylhomoserine (thiol)‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE413,7610AT 6.0% 716.1 / 26.6 217I166I (ATT→ATAmetICystathionine gamma‑synthase/O‑acetylhomoserine (thiol)‑lyase
Reads supporting (aligned to +/- strand):  ref base A (111/91);  new base T (6/7);  total (119/98)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 220 total aligned reads displayed.

GCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATGAAT  >  USA300TCH1516_ALE/413671‑413851
                                                                                          |                                                                                          
 ccAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATaa                                                                                           <  1:3514357/91‑1 (MQ=255)
  cAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAAt                                                                                          >  2:1782524/1‑91 (MQ=255)
    atatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTg                                                                                        >  2:2792751/1‑91 (MQ=255)
     tatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGa                                                                                       >  1:3339788/1‑91 (MQ=255)
     tatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGa                                                                                       >  1:1868060/1‑91 (MQ=255)
     tatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGa                                                                                       >  2:1682184/1‑91 (MQ=255)
        atTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCAATTATCGACaa                                                                                    <  1:954010/91‑1 (MQ=255)
          ttCGTGGCTGAGTGTAAGACTATGTCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGGGTTAAAAAAGTATTGTCAATTATCGACAATa                                                                                  <  2:2965605/91‑1 (MQ=255)
          ttCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATa                                                                                  >  1:2879444/1‑91 (MQ=255)
           tCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATaa                                                                                 <  1:56026/91‑1 (MQ=255)
            cGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAg                                                                                <  2:2193498/91‑1 (MQ=255)
            cGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCAATTATCGACAATact                                                                                <  2:710462/91‑3 (MQ=255)
             gTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGc                                                                               >  1:2190318/1‑91 (MQ=255)
               ggCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCt                                                                             >  2:1712212/1‑91 (MQ=255)
                gCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCtt                                                                            <  2:1861159/91‑1 (MQ=255)
                gCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCtt                                                                            >  2:2397313/1‑91 (MQ=255)
                   gAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAg                                                                         <  2:894020/91‑1 (MQ=255)
                   gAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAg                                                                         <  1:2272126/91‑1 (MQ=255)
                    aaTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGc                                                                        <  2:3407772/91‑1 (MQ=255)
                    aaTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGc                                                                        <  1:286303/91‑1 (MQ=255)
                     aTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCt                                                                       <  2:2719688/91‑1 (MQ=255)
                     aTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCt                                                                       <  2:2719874/91‑1 (MQ=255)
                     aTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCt                                                                       <  1:1038784/91‑1 (MQ=255)
                      tGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCtt                                                                      <  1:3382982/91‑1 (MQ=255)
                        tAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTa                                                                    <  2:1296441/91‑1 (MQ=255)
                        tAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTa                                                                    <  2:191849/91‑1 (MQ=255)
                         aaGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTAc                                                                   >  1:485739/1‑91 (MQ=255)
                          aGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACa                                                                  <  2:1147178/91‑1 (MQ=255)
                           gACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACaa                                                                 >  1:2416795/1‑91 (MQ=255)
                             cTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAg                                                               <  1:2550863/91‑1 (MQ=255)
                              tataTCAGCACCTTCTGCTAGTGGTGTGGAAAGATAAGATGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGt                                                              <  1:108969/91‑1 (MQ=255)
                                 atCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGa                                                           >  2:3164163/1‑91 (MQ=255)
                                  tCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGAt                                                          >  2:2773595/1‑91 (MQ=255)
                                   cAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATa                                                         >  2:1042764/1‑91 (MQ=255)
                                    aGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAg                                                        <  2:1140752/91‑1 (MQ=255)
                                       aCCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTAt                                                     <  1:659703/91‑1 (MQ=255)
                                        ccTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATg                                                    >  2:752202/1‑91 (MQ=255)
                                        ccTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATg                                                    <  2:3737051/91‑1 (MQ=255)
                                         cTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATgg                                                   >  2:300272/1‑91 (MQ=255)
                                         cTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATgg                                                   <  2:2627727/91‑1 (MQ=255)
                                          ttCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGc                                                  >  1:1383829/1‑91 (MQ=255)
                                          ttCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGc                                                  <  2:1381614/91‑1 (MQ=255)
                                           tCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCt                                                 >  2:2080309/1‑91 (MQ=255)
                                           tCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCt                                                 <  2:2603745/91‑1 (MQ=255)
                                            cTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATgg                                                   <  1:2747106/88‑1 (MQ=255)
                                            cTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCtt                                                <  1:2507203/91‑1 (MQ=255)
                                             tGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAACGTTGATAGTATGGCttt                                               >  1:3320513/1‑91 (MQ=255)
                                             tGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCttt                                               >  1:3765608/1‑91 (MQ=255)
                                             tGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCttt                                               >  2:556898/1‑91 (MQ=255)
                                              gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                              >  2:2436800/1‑91 (MQ=255)
                                              gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                              <  2:3632564/91‑1 (MQ=255)
                                              gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                              <  1:1470790/91‑1 (MQ=255)
                                                  gtggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATc                                          >  2:2547149/1‑91 (MQ=255)
                                                  gtggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATc                                          >  2:2548007/1‑91 (MQ=255)
                                                   tggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCa                                         >  1:1595790/1‑91 (MQ=255)
                                                    ggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCaa                                        >  1:794939/1‑91 (MQ=255)
                                                     gtgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAAt                                       <  2:2029909/91‑1 (MQ=255)
                                                      tgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATa                                      <  2:2062976/91‑1 (MQ=255)
                                                         tGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCa                                   <  1:1732705/91‑1 (MQ=255)
                                                         tGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCa                                   <  1:1486658/91‑1 (MQ=255)
                                                          gAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCaa                                  <  2:1583631/91‑1 (MQ=255)
                                                          gAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCaa                                  <  1:3883737/91‑1 (MQ=255)
                                                           aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGCTAGTATGGCTTTACATCAATAGCAAt                                 <  2:1012729/91‑1 (MQ=255)
                                                           aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAAt                                 >  2:2116639/1‑91 (MQ=255)
                                                           aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAAt                                 <  1:1055611/91‑1 (MQ=255)
                                                           aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAAt                                 <  1:2395807/91‑1 (MQ=255)
                                                            aaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATc                                >  1:3125558/1‑91 (MQ=255)
                                                             aGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCa                               <  2:2911704/91‑1 (MQ=255)
                                                             aGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCa                               >  2:2831016/1‑91 (MQ=255)
                                                             aGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCa                               <  1:374817/91‑1 (MQ=255)
                                                              gATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAt                              >  2:303724/1‑91 (MQ=255)
                                                                tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATg                                                    <  1:2778060/67‑1 (MQ=255)
                                                                tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                            >  2:2990276/1‑91 (MQ=255)
                                                                tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                            <  1:677626/91‑1 (MQ=255)
                                                                tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                            <  1:3678064/91‑1 (MQ=255)
                                                                 aaGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTtg                           >  2:3719217/1‑91 (MQ=255)
                                                                 aaGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTtg                           >  2:1224941/1‑91 (MQ=255)
                                                                    gtgtTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTgg                        >  2:1454010/1‑91 (MQ=255)
                                                                        tAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTa                    <  2:2382048/91‑1 (MQ=255)
                                                                         aaaaaaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAg                   <  1:1206553/91‑1 (MQ=255)
                                                                           aaaaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGat                 <  2:1149540/91‑1 (MQ=255)
                                                                           aaaaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGat                 >  2:1413298/1‑91 (MQ=255)
                                                                            aaaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGata                >  2:3670633/1‑91 (MQ=255)
                                                                            aaaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGata                <  2:2665608/91‑1 (MQ=255)
                                                                             aaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGatat               <  2:311324/91‑1 (MQ=255)
                                                                              cttATTGTCAATTATCGACAATAAGCCTTTAGCTTTATAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATAtt              >  1:3277452/3‑91 (MQ=255)
                                                                                tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg            >  2:1826579/1‑91 (MQ=255)
                                                                                tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg            >  1:413239/1‑91 (MQ=255)
                                                                                tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTACATATTgg            >  1:1948703/1‑91 (MQ=255)
                                                                                 aTTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGt           >  2:1060340/1‑91 (MQ=255)
                                                                                  ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTCCATCAATAGCAATCATTTGTGGGTTAGATATTGGtt          >  2:3506602/1‑91 (MQ=255)
                                                                                  ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGtt          >  1:2580947/1‑91 (MQ=255)
                                                                                  ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGtt          >  1:2463251/1‑91 (MQ=255)
                                                                                  ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGtt          >  1:684139/1‑91 (MQ=255)
                                                                                    gTCGATAATTGAGAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCa        >  1:2974029/1‑91 (MQ=255)
                                                                                    gTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCa        >  2:464869/1‑91 (MQ=255)
                                                                                      cGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAAt      <  2:2704938/91‑1 (MQ=255)
                                                                                      cGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAAt      >  2:3821857/1‑91 (MQ=255)
                                                                                       gATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATg     >  1:1125039/1‑91 (MQ=255)
                                                                                        aTAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATGa    <  1:3594040/91‑1 (MQ=255)
                                                                                          |                                                                                          
GCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATGAAT  >  USA300TCH1516_ALE/413671‑413851

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: