Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,273,482 T→A 6.4% F136L (TTT→TTA mtlA → PTS system mannitol‑specific EIICB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,273,4820TA 6.4% 630.6 / 24.9 187F136L (TTT→TTAmtlAPTS system mannitol‑specific EIICB component
Reads supporting (aligned to +/- strand):  ref base T (91/84);  new base A (5/7);  total (96/91)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 211 total aligned reads displayed.

ATGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATAT  >  USA300TCH1516_ALE/2273393‑2273571
                                                                                         |                                                                                         
aTGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTta                                                                                          <  1:428663/91‑1 (MQ=255)
 tGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaa                                                                                         <  2:477664/91‑1 (MQ=255)
 tGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaa                                                                                         >  2:3071322/1‑91 (MQ=255)
 tGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaa                                                                                         >  2:170525/1‑91 (MQ=255)
  gATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaat                                                                                        <  1:365923/91‑1 (MQ=255)
  gATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaat                                                                                        <  2:3546429/91‑1 (MQ=255)
    ttATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaataa                                                                                      <  1:1106718/91‑1 (MQ=255)
    ttATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaataa                                                                                      <  2:1679674/91‑1 (MQ=255)
    ttATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaataa                                                                                      <  2:2634897/91‑1 (MQ=255)
      aTGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACt                                                                                    <  1:1102523/91‑1 (MQ=255)
      aTGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACt                                                                                    <  2:1963855/91‑1 (MQ=255)
       tGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAAAGTTATTAAATAAcat                                                                                   <  1:3152578/91‑3 (MQ=255)
          gTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAAAGTTATTAAATAACATTTc                                                                                <  1:820774/91‑1 (MQ=255)
           tCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTct                                                                               <  2:2906405/91‑1 (MQ=255)
            ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                              >  2:2055396/1‑91 (MQ=255)
            ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                              <  2:1506654/91‑1 (MQ=255)
            ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                              <  2:3727838/91‑1 (MQ=255)
             cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgc                                                                             <  2:809871/91‑1 (MQ=255)
             cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgc                                                                             <  1:1073639/91‑1 (MQ=255)
             cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgc                                                                             <  1:2115140/91‑1 (MQ=255)
             cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAAAGTTATTAAATAACATTTcaaa                                                                             <  1:258729/91‑4 (MQ=38)
              aCTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgct                                                                            <  1:1313133/91‑1 (MQ=255)
                ttgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTgg                                                                          <  2:328950/91‑1 (MQ=255)
                 tgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGt                                                                         >  1:216648/1‑91 (MQ=255)
                 tgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGt                                                                         >  2:2723640/1‑91 (MQ=255)
                  gttggttgGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTa                                                                        <  2:573787/91‑1 (MQ=255)
                   ttggttggTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAt                                                                       >  2:1576194/1‑91 (MQ=255)
                   ttggttggTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAt                                                                       <  2:3135388/91‑1 (MQ=255)
                   ttggttggTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAt                                                                       <  1:3449173/91‑1 (MQ=255)
                       ttggttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaataa                                                                                      <  2:2966317/72‑1 (MQ=255)
                       ttggttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTtaataa                                                                                      >  1:2966317/1‑72 (MQ=255)
                       ttggttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTa                                                                   <  1:2959022/91‑1 (MQ=255)
                         ggttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAgg                                                                 >  2:3876190/1‑91 (MQ=255)
                            tAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAtt                                                              >  2:1152944/1‑91 (MQ=255)
                              aTGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTa                                                            <  2:3182612/91‑1 (MQ=255)
                               tGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTAt                                                           >  2:283976/1‑91 (MQ=255)
                               tGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTAt                                                           <  2:3686350/91‑1 (MQ=255)
                                gaagaaGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATc                                                          <  1:723954/91‑1 (MQ=255)
                                      gaccgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACg                                                    <  1:390882/91‑1 (MQ=255)
                                       accgaccgATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGa                                                   <  1:1957045/91‑1 (MQ=255)
                                       accgaccAATTGATTCAACCTAGAACGCCGCAAGTTGTTGAAATGTTATTTAATAACTCTTCTGCTGGTATTTTAGGATTTATCATGACGa                                                   <  2:3804657/91‑1 (MQ=255)
                                        ccgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCAGCTGGTATTTTAGGATTTATCATGACGAt                                                  >  1:2161603/1‑91 (MQ=255)
                                          gaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTg                                                >  2:2189731/1‑91 (MQ=255)
                                           accAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGc                                               >  1:299397/1‑91 (MQ=255)
                                           accAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGc                                               <  2:2782443/91‑1 (MQ=255)
                                            ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                              >  2:1793097/1‑91 (MQ=255)
                                            ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                              >  2:1357976/1‑91 (MQ=255)
                                            ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                              <  2:2814763/91‑1 (MQ=255)
                                             cAATTGATTCAACCTAGAGCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAg                                             <  2:2751597/91‑1 (MQ=255)
                                             cAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAg                                             >  2:2144865/1‑91 (MQ=255)
                                               aTTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGa                                           >  2:1785430/1‑91 (MQ=255)
                                                 tGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATTACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAtt                                         <  1:3078395/91‑1 (MQ=255)
                                                 tGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAtt                                         <  2:3535520/91‑1 (MQ=255)
                                                  gATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAttt                                        >  1:2938806/1‑91 (MQ=255)
                                                  gATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAttt                                        <  1:3750177/91‑1 (MQ=255)
                                                   aTTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTa                                       <  1:782277/91‑1 (MQ=255)
                                                   aTTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTa                                       <  1:5179/91‑1 (MQ=255)
                                                     tCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaa                                     <  1:2297382/91‑1 (MQ=255)
                                                      cAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaaa                                    >  2:1300251/1‑91 (MQ=255)
                                                       aaCCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAt                                   >  1:765621/1‑91 (MQ=255)
                                                        aCCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTCGGATTTATCATGACGATTGCAGGATTTAAAAtt                                  >  1:3183792/1‑91 (MQ=255)
                                                        aCCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAtt                                  >  2:509803/1‑91 (MQ=255)
                                                           tAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTa                               >  1:3575789/1‑91 (MQ=255)
                                                           tAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTa                               >  1:2671058/1‑91 (MQ=255)
                                                            aGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAg                              <  1:394394/91‑1 (MQ=255)
                                                             gAACGCCGCAAGGTTTTGAAATGTTATTTAATAACGTTTCTGCTGGTATTTTAGGCTTTATCATGCCGATTGCAGGATTTACAATCTTAGc                             >  1:1399308/1‑91 (MQ=255)
                                                              aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                            >  1:1063264/1‑91 (MQ=255)
                                                              aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                            <  1:1647193/91‑1 (MQ=255)
                                                              aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                            >  1:927060/1‑91 (MQ=255)
                                                              aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                            <  2:2440484/91‑1 (MQ=255)
                                                                cgccgcACGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacc                          >  1:2731657/1‑91 (MQ=255)
                                                                cgccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacc                          >  1:157354/1‑91 (MQ=255)
                                                                cgccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCGGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacc                          >  1:2513088/1‑91 (MQ=255)
                                                                 gccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacca                         >  1:9964/1‑91 (MQ=255)
                                                                  ccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccac                        >  1:1218462/1‑91 (MQ=255)
                                                                  ccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccac                        <  2:3592753/91‑1 (MQ=255)
                                                                  ccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccac                        >  2:2705735/1‑91 (MQ=255)
                                                                   cgcATGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACt                       <  2:3576023/91‑1 (MQ=255)
                                                                     tctggTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCTCCACTTa                     >  2:3451743/4‑91 (MQ=255)
                                                                     cAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                     >  2:62475/1‑91 (MQ=255)
                                                                     cAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                     <  2:3500568/91‑1 (MQ=255)
                                                                     cAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAGTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                     >  1:3025755/1‑91 (MQ=255)
                                                                     cAAGGTTTCGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                     <  2:3371982/91‑1 (MQ=255)
                                                                      aaGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTAt                    >  2:1888203/1‑91 (MQ=255)
                                                                        ggTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGa                  <  1:2761590/91‑1 (MQ=255)
                                                                         gTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaa                 <  1:1981204/91‑1 (MQ=255)
                                                                          ttttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaaa                <  1:2250575/91‑1 (MQ=255)
                                                                          ttttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaaa                <  2:764286/91‑1 (MQ=255)
                                                                           tttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAt               >  2:3301204/1‑91 (MQ=255)
                                                                           tttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAt               <  1:170525/91‑1 (MQ=255)
                                                                            ttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAtt              <  2:1204250/91‑1 (MQ=255)
                                                                            ttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAtt              >  2:3160452/1‑91 (MQ=255)
                                                                              gAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCTCTTATGAAATtta            >  2:411658/1‑91 (MQ=255)
                                                                              gAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATtta            >  2:2614285/1‑91 (MQ=255)
                                                                                 aagTTATTAAATAACATTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATttatta         >  2:2455811/3‑91 (MQ=255)
                                                                                  tGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTtattat        >  1:863262/1‑91 (MQ=255)
                                                                                   gTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATg       <  1:1813308/91‑1 (MQ=255)
                                                                                   gTTATTAAATAACATTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATg       >  1:3282430/1‑91 (MQ=255)
                                                                                    ttATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGc      <  2:406657/91‑1 (MQ=255)
                                                                                        ttAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCatat  >  2:3401442/1‑91 (MQ=255)
                                                                                         |                                                                                         
ATGATTATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATAT  >  USA300TCH1516_ALE/2273393‑2273571

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: