Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,049,713 G→T 5.3% Q457H (CAG→CAT menD → 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,049,7130GT 5.3% 745.4 / 20.9 208Q457H (CAG→CATmenD2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase
Reads supporting (aligned to +/- strand):  ref base G (114/83);  new base T (3/8);  total (117/91)
Fisher's exact test for biased strand distribution p-value = 6.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 210 total aligned reads displayed.

TGTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTG  >  USA300TCH1516_ALE/1049623‑1049803
                                                                                          |                                                                                          
 gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                           >  1:1080908/1‑91 (MQ=255)
 gtgCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGa                                                                                           <  2:2799187/91‑1 (MQ=255)
    cATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGa                                                                                        >  1:1516992/1‑91 (MQ=255)
       aaaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAata                                                                                     >  2:2279348/1‑91 (MQ=255)
       aaaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAata                                                                                     >  2:1011321/1‑91 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAatat                                                                                    <  1:1604209/91‑1 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  1:36295/91‑1 (MQ=255)
        aaCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAatat                                                                                    <  2:2848473/91‑1 (MQ=255)
         aCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAtt                                                                                   <  1:1450637/91‑1 (MQ=255)
         aCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATAtt                                                                                   >  1:575809/1‑91 (MQ=255)
         aCGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATAtt                                                                                   <  2:273069/91‑1 (MQ=255)
          cGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTg                                                                                  >  1:3288211/1‑91 (MQ=255)
            aaTAAGATTATTGATAGGTGTTTTGTCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTCAGTAATATTCATCTGAATAttatg                                                                                <  2:1556392/91‑4 (MQ=255)
            aaTAACATTAGTGATAGGTGATTTATCATTTTATCATGATATGACTGGACTATTAATGTCAAAATTAAATAATATTCATCTGAATAttatt                                                                                <  2:263351/91‑4 (MQ=255)
             aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtat                                                                               >  2:3107961/1‑91 (MQ=255)
             aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtat                                                                               >  2:2616105/1‑91 (MQ=255)
             aTAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtat                                                                               >  1:1281736/1‑91 (MQ=255)
              tAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatt                                                                              >  2:1492276/1‑91 (MQ=255)
              tAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatt                                                                              >  1:825426/1‑91 (MQ=255)
              tAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatt                                                                              >  1:2424878/1‑91 (MQ=255)
              tAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatt                                                                              >  2:3152358/1‑91 (MQ=255)
               aaCATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtatta                                                                             <  1:3614122/91‑1 (MQ=255)
                aCATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtattat                                                                            >  1:3440826/1‑91 (MQ=255)
                aCATTATTGATAAGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGtagtat                                                                            >  1:3061991/1‑91 (MQ=255)
                   ttattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGa                                                                         <  2:729412/91‑1 (MQ=255)
                   ttattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGa                                                                         <  1:1615633/91‑1 (MQ=255)
                   ttattGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAAGATTCAGATGAATATTGTATTATTGa                                                                         >  2:3419477/1‑91 (MQ=255)
                     attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                       <  2:1125574/91‑1 (MQ=255)
                     attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                       <  2:1811775/91‑1 (MQ=255)
                     attGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaac                                                                       >  1:2912386/1‑91 (MQ=255)
                       tGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaa                                                                     <  2:3187931/91‑1 (MQ=255)
                        gATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaac                                                                    >  1:1025944/1‑91 (MQ=255)
                        gATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGaacaac                                                                    >  1:225568/1‑91 (MQ=255)
                         aTAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACg                                                                   >  2:1135000/1‑91 (MQ=255)
                         aTAGGTGATTTATCATGTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACg                                                                   <  2:3771377/91‑1 (MQ=255)
                          tAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGa                                                                  >  1:477219/1‑91 (MQ=255)
                          tAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGa                                                                  >  1:3920678/1‑91 (MQ=255)
                           aGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAt                                                                 <  1:2561397/91‑1 (MQ=255)
                             gTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtgg                                                               <  1:1130115/91‑1 (MQ=255)
                                aTTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGAtggtgg                                                            >  2:402543/1‑91 (MQ=255)
                                 tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                           >  2:440189/1‑91 (MQ=255)
                                 tttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGc                                                           >  2:1901927/1‑91 (MQ=255)
                                  ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                          <  1:2551909/91‑1 (MQ=255)
                                  ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                          >  1:2675342/1‑91 (MQ=255)
                                  ttatcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                          >  2:886755/1‑91 (MQ=255)
                                  ttatcattttatcatGATATGAATGCACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCg                                                          <  1:3664057/91‑1 (MQ=255)
                                    atcattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGt                                                        >  1:2220503/1‑91 (MQ=255)
                                      cattttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAt                                                      >  1:2370101/1‑91 (MQ=255)
                                       attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtt                                                     >  1:3250225/1‑91 (MQ=255)
                                       attttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtt                                                     >  1:2757746/1‑91 (MQ=255)
                                        ttttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAttt                                                    >  2:2059677/1‑91 (MQ=255)
                                         tttatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttt                                                   >  1:3478198/1‑91 (MQ=255)
                                           tatcatGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTAtttttt                                                 <  1:2476327/91‑1 (MQ=255)
                                              catGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCat                                              <  1:1318524/91‑1 (MQ=255)
                                                 gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAtt                                           >  1:2958318/1‑91 (MQ=255)
                                                 gATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAtt                                           >  1:3008621/1‑91 (MQ=255)
                                                  atatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAttt                                          <  1:2173044/91‑1 (MQ=255)
                                                   tatGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTa                                         <  1:1656901/91‑1 (MQ=255)
                                                     tGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTAcc                                       <  1:1001232/91‑1 (MQ=255)
                                                      gAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTGTTTTTTCATATTTACca                                      <  2:2569385/91‑1 (MQ=255)
                                                      gAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACca                                      >  1:1492368/1‑91 (MQ=255)
                                                       aaTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTaaaac                                     <  2:2960768/91‑5 (MQ=255)
                                                       aaTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcac                                     <  2:3340300/91‑1 (MQ=255)
                                                       aaTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcac                                     <  1:1961744/91‑1 (MQ=255)
                                                        aTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcaca                                    >  2:913431/1‑91 (MQ=255)
                                                        aTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcaca                                    >  2:3306538/1‑91 (MQ=255)
                                                        aTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcaca                                    >  1:183013/1‑91 (MQ=255)
                                                        aTGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACcaca                                    >  1:2317980/1‑91 (MQ=255)
                                                          ggACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  <  2:45689/91‑1 (MQ=255)
                                                          ggACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACaaa                                  >  1:1727464/1‑91 (MQ=255)
                                                            aCTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAAAGATGGTGGCGGTATTTTTTCATATTTACCAcaacaa                                >  2:1340790/1‑91 (MQ=255)
                                                             cTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaag                               <  1:2000952/91‑1 (MQ=255)
                                                               aTTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaa                             >  2:3472870/1‑91 (MQ=255)
                                                                ttAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAaaagaaa                            >  1:3182625/1‑91 (MQ=255)
                                                                  aaTGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGt                          <  2:361462/91‑1 (MQ=255)
                                                                    tGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGc                        <  1:3727336/91‑1 (MQ=255)
                                                                     gTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCa                       >  2:1458989/1‑91 (MQ=255)
                                                                      tCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCaa                      <  1:1392078/91‑1 (MQ=255)
                                                                       cAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAAc                     <  2:3883684/91‑1 (MQ=255)
                                                                        aaaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACt                    <  2:2457137/91‑1 (MQ=255)
                                                                         aaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTg                   >  2:1228611/1‑91 (MQ=255)
                                                                         aaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTg                   <  2:1952219/91‑1 (MQ=255)
                                                                         aaaTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTg                   >  1:2029862/1‑91 (MQ=255)
                                                                           aTTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATAATAAACAAAAAAAGAAAGTGCAACTGAc                 <  2:1476858/91‑1 (MQ=255)
                                                                             tAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTa               >  1:3691025/1‑91 (MQ=255)
                                                                              aaaTAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAt              >  1:1852795/1‑91 (MQ=255)
                                                                               aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtt             >  2:3348483/1‑91 (MQ=255)
                                                                               aataatATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtt             >  1:2598908/1‑91 (MQ=255)
                                                                                ataataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt            <  2:2917233/91‑1 (MQ=255)
                                                                                ataataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAttt            >  2:841296/1‑91 (MQ=255)
                                                                                 taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCTTATTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  2:1424355/1‑91 (MQ=255)
                                                                                 taataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  2:3320182/1‑91 (MQ=255)
                                                                                 caataTTCATCTGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTAtttt           >  1:2782609/2‑91 (MQ=255)
                                                                                  aataTTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTg          >  2:1724376/1‑91 (MQ=255)
                                                                                   atatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGa         >  2:64954/1‑91 (MQ=255)
                                                                                    tatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGaa        <  1:886824/91‑1 (MQ=255)
                                                                                    tatTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGaa        >  2:1681454/1‑91 (MQ=255)
                                                                                      ttCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACg      <  1:2995034/91‑1 (MQ=255)
                                                                                      ttCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACg      <  2:2705727/91‑1 (MQ=255)
                                                                                        cAGATGATTATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCTTATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGgt    >  2:1875456/1‑91 (MQ=255)
                                                                                          |                                                                                          
TGTGCATAAACGAATAACATTATTGATAGGTGATTTATCATTTTATCATGATATGAATGGACTATTAATGTCAAAATTAAATAATATTCAGATGAATATTGTATTATTGAACAACGATGGTGGCGGTATTTTTTCATATTTACCACAAAAAGAAAGTGCAACTGACTATTTTGAACGGTTG  >  USA300TCH1516_ALE/1049623‑1049803

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: