Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 413,758 G→A 5.7% I167I (ATC→ATT 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,7580GA 5.7% 731.7 / 18.9 222I167I (ATC→ATTmetICystathionine gamma‑synthase/O‑acetylhomoserine (thiol)‑lyase
Reads supporting (aligned to +/- strand):  ref base G (114/94);  new base A (6/7);  total (121/101)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

TCCGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATG  >  USA300TCH1516_ALE/413668‑413848
                                                                                          |                                                                                          
tCCGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCg                                                                                            >  2:263398/1‑91 (MQ=255)
tCCGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCg                                                                                            <  1:1251743/91‑1 (MQ=255)
tCCGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCg                                                                                            >  2:3319101/1‑91 (MQ=255)
 ccGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGa                                                                                           >  2:2404529/1‑91 (MQ=255)
  cGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGAt                                                                                          <  1:1394240/91‑1 (MQ=255)
  cGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGAt                                                                                          >  1:25687/1‑91 (MQ=255)
     cAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAAt                                                                                       >  2:1782524/1‑91 (MQ=255)
       atatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTg                                                                                     >  2:2792751/1‑91 (MQ=255)
       atatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTg                                                                                     <  2:284652/91‑1 (MQ=255)
        tatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGa                                                                                    >  2:1682184/1‑91 (MQ=255)
        tatatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGa                                                                                    >  1:3339788/1‑91 (MQ=255)
          tatTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCAATTATCGACa                                                                                  <  2:3545075/91‑1 (MQ=255)
           atTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACaa                                                                                 >  2:3026911/1‑91 (MQ=255)
           atTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCAATTATCGACaa                                                                                 <  1:954010/91‑1 (MQ=255)
               cGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAg                                                                             <  2:2193498/91‑1 (MQ=255)
               cGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCAATTATCGACAATact                                                                             <  2:710740/91‑3 (MQ=255)
                  ggCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCt                                                                          >  2:1712212/1‑91 (MQ=255)
                   gCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCtt                                                                         <  2:1860934/91‑1 (MQ=255)
                   gCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCtt                                                                         <  2:2486258/91‑1 (MQ=255)
                      gAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAg                                                                      >  1:421465/1‑91 (MQ=255)
                       aaTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGc                                                                     <  2:3407772/91‑1 (MQ=255)
                        aTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCt                                                                    <  2:2719855/91‑1 (MQ=255)
                        aTGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCt                                                                    <  1:1038784/91‑1 (MQ=255)
                         tGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCtt                                                                   <  2:1464935/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)
                             aGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACa                                                               <  2:1147178/91‑1 (MQ=255)
                              gACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACaa                                                              >  1:2417122/1‑91 (MQ=255)
                              caCTATATCAGCACCTTCTGCTAGCGGTGTTGAAAGCTAAGGTGTTAAAAAAGTATTGTCGATGATTGACAATAAGCCTTGAGCTTTACaa                                                              >  2:761121/2‑91 (MQ=255)
                               aCTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTCAGCTTTACaaa                                                             >  1:3522417/1‑91 (MQ=255)
                                cTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAg                                                            <  1:2550863/91‑1 (MQ=255)
                                cTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAg                                                            >  1:261162/1‑91 (MQ=255)
                                 tataTCAGCACCTTCTGCTAGTGGTGTGGAAAGATAAGATGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGt                                                           <  1:108969/91‑1 (MQ=255)
                                  atatCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGtt                                                          >  1:3785765/1‑91 (MQ=255)
                                  atatCAGCACCTTCTACTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGAAAATTGACAATAAGCCTTTAGCTTTACAAAGtt                                                          >  1:3306664/1‑91 (MQ=255)
                                    atCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGa                                                        >  2:3164163/1‑91 (MQ=255)
                                     tCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGAt                                                       >  2:2712665/1‑91 (MQ=255)
                                      cAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATa                                                      >  2:1042628/1‑91 (MQ=255)
                                       aGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAg                                                     <  2:1140752/91‑1 (MQ=255)
                                       aGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAg                                                     >  1:309591/1‑91 (MQ=255)
                                        gCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGt                                                    <  2:2954301/91‑1 (MQ=255)
                                           ccTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATg                                                 <  2:3735627/91‑1 (MQ=255)
                                            cTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATgg                                                >  2:300272/1‑91 (MQ=255)
                                            cTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATgg                                                <  2:2627727/91‑1 (MQ=255)
                                             ttCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGc                                               >  1:2201408/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                                                >  2:2747106/1‑88 (MQ=255)
                                               cTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCtt                                             >  1:50581/1‑91 (MQ=255)
                                               cTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCtt                                             <  1:2298332/91‑1 (MQ=255)
                                                 gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                           <  2:3632564/91‑1 (MQ=255)
                                                 gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                           >  2:2704151/1‑91 (MQ=255)
                                                 gCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTa                                           >  2:2436800/1‑91 (MQ=255)
                                                     gtggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATc                                       >  2:2547951/1‑91 (MQ=255)
                                                      tggtgTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCa                                      >  1:1595790/1‑91 (MQ=255)
                                                        gtgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAAt                                    >  2:1043614/1‑91 (MQ=255)
                                                        gtgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAAt                                    <  2:2029909/91‑1 (MQ=255)
                                                         tgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATa                                   <  2:2062976/91‑1 (MQ=255)
                                                         tgtTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATa                                   <  1:2494791/91‑1 (MQ=255)
                                                           ttGAAAGATGAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGc                                 <  1:3714046/91‑1 (MQ=255)
                                                            tGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCa                                <  1:1732705/91‑1 (MQ=255)
                                                            tGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCa                                <  1:3226001/91‑1 (MQ=255)
                                                            tGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCa                                <  1:1486658/91‑1 (MQ=255)
                                                             gAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCaa                               <  2:1583631/91‑1 (MQ=255)
                                                             gAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCaa                               >  1:614216/1‑91 (MQ=255)
                                                              aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGCTAGTATGGCTTTACATCAATAGCAAt                              <  2:1012729/91‑1 (MQ=255)
                                                              aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAAt                              <  1:1055611/91‑1 (MQ=255)
                                                              aaaGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAAt                              >  2:2116639/1‑91 (MQ=255)
                                                                  aTAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAtt                          >  1:3758603/1‑91 (MQ=255)
                                                                   tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                         >  1:1464171/1‑91 (MQ=255)
                                                                   tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                         >  2:2990276/1‑91 (MQ=255)
                                                                   tAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCAttt                         <  1:678032/91‑1 (MQ=255)
                                                                    aaGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTtg                        >  2:1224941/1‑91 (MQ=255)
                                                                      ggTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTtgtg                      >  1:2871792/1‑91 (MQ=255)
                                                                      ggTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTtgtg                      >  1:2412366/1‑91 (MQ=255)
                                                                       gtgtTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTgg                     >  1:3467674/1‑91 (MQ=255)
                                                                       gtgtTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATAAATAGCAATCATTTGTgg                     >  1:2752303/1‑91 (MQ=255)
                                                                         gtTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATTGTATGGCTTTACATCAATAGCAATCATTTGTGGGt                   >  2:684298/1‑91 (MQ=255)
                                                                         gtTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGt                   <  1:1335762/91‑1 (MQ=255)
                                                                           tAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTa                 <  2:2382048/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)
                                                                                aaGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGatat            <  2:311324/91‑1 (MQ=255)
                                                                                   tATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg         >  2:2824384/1‑91 (MQ=255)
                                                                                   tATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg         <  1:2792751/91‑1 (MQ=255)
                                                                                   tATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg         >  1:1875954/1‑91 (MQ=255)
                                                                                   tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTATATATTgg         >  1:1948724/1‑91 (MQ=255)
                                                                                   tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg         >  2:1826579/1‑91 (MQ=255)
                                                                                   tATTGTCAATTATCGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTgg         >  1:413095/1‑91 (MQ=255)
                                                                                    aTTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGt        <  2:3785765/91‑1 (MQ=255)
                                                                                     ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGtt       >  1:2580947/1‑91 (MQ=255)
                                                                                     ttGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGtt       >  1:674060/1‑91 (MQ=255)
                                                                                       gTCGATAATTGAGAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCa     >  1:2974029/1‑91 (MQ=255)
                                                                                       gTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCa     >  2:464869/1‑91 (MQ=255)
                                                                                       gTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCa     >  1:696574/1‑91 (MQ=255)
                                                                                         cGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAAt   >  2:3822752/1‑91 (MQ=255)
                                                                                          gATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATg  >  1:1125039/1‑91 (MQ=255)
                                                                                          |                                                                                          
TCCGCCAATATATTTCGTGGCTGAATGTAAGACTATATCAGCACCTTCTGCTAGTGGTGTTGAAAGATAAGGTGTTAAAAAAGTATTGTCGATAATTGACAATAAGCCTTTAGCTTTACAAAGTTGATAGTATGGCTTTACATCAATAGCAATCATTTGTGGGTTAGATATTGGTTCAATG  >  USA300TCH1516_ALE/413668‑413848

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: