Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,273,489 T→A 5.6% F139I (TTT→ATT)  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,4890TA 5.6% 660.2 / 22.5 198F139I (TTT→ATT) mtlAPTS system mannitol‑specific EIICB component
Reads supporting (aligned to +/- strand):  ref base T (94/92);  new base A (9/2);  total (104/94)
Fisher's exact test for biased strand distribution p-value = 6.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 205 total aligned reads displayed.

ATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTATCAG  >  USA300TCH1516_ALE/2273399‑2273579
                                                                                          |                                                                                          
aTGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACt                                                                                            <  2:1963855/91‑1 (MQ=255)
aTGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACt                                                                                            >  1:2832332/1‑91 (MQ=255)
    gTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTc                                                                                        <  2:3550492/91‑1 (MQ=255)
     tCCCCTTGGTGGTTGGTTTATGAAGCAGACCGAGCAATTGATGCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTct                                                                                       <  1:3182918/91‑1 (MQ=255)
     tCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTct                                                                                       <  2:2906405/91‑1 (MQ=255)
      ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                                      <  2:3727838/91‑1 (MQ=255)
      ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                                      >  2:2055396/1‑91 (MQ=255)
      ccACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctg                                                                                      <  2:1506654/91‑1 (MQ=255)
       cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgc                                                                                     <  2:809871/91‑1 (MQ=255)
       cACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgc                                                                                     <  1:1073639/91‑1 (MQ=255)
        aCTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTctgct                                                                                    <  1:1312952/91‑1 (MQ=255)
          ttgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACTCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTgg                                                                                  <  2:2948725/91‑1 (MQ=255)
          ttgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTgg                                                                                  <  2:328950/91‑1 (MQ=255)
           tgttgGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGt                                                                                 >  1:216648/1‑91 (MQ=255)
            gttggttgGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTa                                                                                <  2:573787/91‑1 (MQ=255)
             ttggttggTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAt                                                                               >  2:1576194/1‑91 (MQ=255)
             ttggttggTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAt                                                                               <  2:3135417/91‑1 (MQ=255)
              tggttggtTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTAtt                                                                              <  1:2950628/91‑1 (MQ=255)
                   ggttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAgg                                                                         >  1:121041/1‑91 (MQ=255)
                     ttAATGAAGAAGACCGACCAATTGATTCAACCTAGGACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAt                                                                       <  1:3414074/91‑1 (MQ=255)
                     ttAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAt                                                                       >  2:2732875/1‑91 (MQ=255)
                      tAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAtt                                                                      >  2:1152944/1‑91 (MQ=255)
                      tAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAtt                                                                      >  2:2639818/1‑91 (MQ=255)
                      tAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGAtt                                                                      <  1:25546/91‑1 (MQ=255)
                        aTGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTa                                                                    <  2:3182612/91‑1 (MQ=255)
                         tGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTAt                                                                   >  2:283976/1‑91 (MQ=255)
                         tGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTAt                                                                   <  2:3686350/91‑1 (MQ=255)
                          gaagaaGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATc                                                                  >  2:2404729/1‑91 (MQ=255)
                            agaagaCCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCAt                                                                >  1:3648208/1‑91 (MQ=255)
                                gaccgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACg                                                            <  1:390882/91‑1 (MQ=255)
                                 accgaccgATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGa                                                           <  1:1957045/91‑1 (MQ=255)
                                 accgaccAATTGATTCAACCTAGAACGCCGCAAGTTGTTGAAATGTTATTTAATAACTCTTCTGCTGGTATTTTAGGATTTATCATGACGa                                                           <  2:3804657/91‑1 (MQ=255)
                                 accgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACgg                                                           >  2:3844774/1‑90 (MQ=255)
                                 accgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGa                                                           <  1:1578057/91‑1 (MQ=255)
                                  ccgaccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCAGCTGGTATTTTAGGATTTATCATGACGAt                                                          >  1:2161603/1‑91 (MQ=255)
                                     accAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGc                                                       <  2:2782443/91‑1 (MQ=255)
                                      ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                                      >  1:1358460/1‑91 (MQ=255)
                                      ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                                      <  2:2814763/91‑1 (MQ=255)
                                      ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                                      >  2:1357976/1‑91 (MQ=255)
                                      ccAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCa                                                      >  2:1793097/1‑91 (MQ=255)
                                       cAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAg                                                     >  2:2144865/1‑91 (MQ=255)
                                       cAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAg                                                     >  2:3079741/1‑91 (MQ=255)
                                         aTTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGa                                                   >  2:1785430/1‑91 (MQ=255)
                                          ttGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAt                                                  >  1:2651180/1‑91 (MQ=255)
                                          ttGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAt                                                  <  1:3025402/91‑1 (MQ=255)
                                           tGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATTACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAtt                                                 <  1:3078395/91‑1 (MQ=255)
                                           tGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAtt                                                 <  2:3535520/91‑1 (MQ=255)
                                            gATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGAttt                                                <  1:3750177/91‑1 (MQ=255)
                                             aTTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTa                                               <  1:5179/91‑1 (MQ=255)
                                               tCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaa                                             <  1:2297382/91‑1 (MQ=255)
                                               tCAACCTAGAACGCCGCAAGGTTTTGAAATGCTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaa                                             >  1:3211800/1‑91 (MQ=255)
                                                cAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaaa                                            >  2:1300251/1‑91 (MQ=255)
                                                cAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaaa                                            >  2:399028/1‑91 (MQ=255)
                                                cAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTaaaa                                            <  2:2715634/91‑1 (MQ=255)
                                                 aaCCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAt                                           >  1:766428/1‑91 (MQ=255)
                                                  aCCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAtt                                          >  2:509803/1‑91 (MQ=255)
                                                   ccTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAttt                                         >  1:3572248/1‑91 (MQ=255)
                                                    cTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAAtttt                                        >  2:408915/1‑91 (MQ=255)
                                                      aGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAg                                      >  1:929199/1‑91 (MQ=255)
                                                      aGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAg                                      <  1:393763/91‑1 (MQ=255)
                                                        aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                                    >  1:1063264/1‑91 (MQ=255)
                                                        aaCGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGca                                    <  2:2440504/91‑1 (MQ=255)
                                                          cgccgcACGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacc                                  >  1:2731657/1‑91 (MQ=255)
                                                          cgccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCGGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacc                                  >  1:2513088/1‑91 (MQ=255)
                                                           gccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacca                                 >  1:2967043/1‑91 (MQ=255)
                                                           gccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcacca                                 >  1:9964/1‑91 (MQ=255)
                                                           gccgcAAGGTTTTGAAATGTTATTTAATAACTATTCTGCGGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccc                                 >  2:2467971/1‑90 (MQ=255)
                                                            ccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccac                                >  1:1104650/1‑91 (MQ=255)
                                                            ccgcAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGcaccac                                <  2:3592753/91‑1 (MQ=255)
                                                             cgcATGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACt                               <  2:3576023/91‑1 (MQ=255)
                                                               tctggTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCTCCACTTa                             >  2:3451743/4‑91 (MQ=255)
                                                               cAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                             <  2:3500568/91‑1 (MQ=255)
                                                               cAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                             >  1:3025897/1‑91 (MQ=255)
                                                               cAAGGTTTCGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTa                             <  2:3371982/91‑1 (MQ=255)
                                                                aaGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTAt                            >  2:1888203/1‑91 (MQ=255)
                                                                  ggTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGa                          >  2:277209/1‑91 (MQ=255)
                                                                  ggTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGa                          <  1:2761590/91‑1 (MQ=255)
                                                                   gTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaa                         <  1:1981204/91‑1 (MQ=255)
                                                                    ttttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaaa                        <  2:764286/91‑1 (MQ=255)
                                                                    ttttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGaaa                        <  1:2250575/91‑1 (MQ=255)
                                                                     tttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAt                       <  1:170525/91‑1 (MQ=255)
                                                                     tttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAt                       >  2:3301204/1‑91 (MQ=255)
                                                                      ttGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAtt                      >  2:1715274/1‑91 (MQ=255)
                                                                       tGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAAttt                     <  1:801352/91‑1 (MQ=255)
                                                                        gAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCTCTTATGAAATtta                    >  2:411674/1‑91 (MQ=255)
                                                                        gAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATtta                    >  1:401734/1‑91 (MQ=255)
                                                                        gAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATtta                    >  2:2614285/1‑91 (MQ=255)
                                                                            tGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTtattat                <  1:1429543/91‑1 (MQ=255)
                                                                            tGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTtattat                >  1:863867/1‑91 (MQ=255)
                                                                            agTTATTAAATAACATTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTtattat                >  1:8223/2‑91 (MQ=255)
                                                                             gTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATg               <  1:1813308/91‑1 (MQ=255)
                                                                                 ttaaataaCATTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCata           >  2:1895782/4‑91 (MQ=255)
                                                                                  ttAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCatat          >  2:3401702/1‑91 (MQ=255)
                                                                                  taaataaCATTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCatat          >  2:2527817/3‑91 (MQ=255)
                                                                                     ataaCTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATAtttt       <  2:474202/91‑1 (MQ=255)
                                                                                      taaCTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTa      <  2:121041/91‑1 (MQ=255)
                                                                                       aaCTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTAt     >  1:641817/1‑91 (MQ=255)
                                                                                         cTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTATCa   <  1:2125178/91‑1 (MQ=255)
                                                                                         cTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTATCa   >  2:683113/1‑91 (MQ=255)
                                                                                          ttttCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTATCAg  >  1:727706/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATGGGTCCACTTGTTGGTTGGTTAATGAAGAAGACCGACCAATTGATTCAACCTAGAACGCCGCAAGGTTTTGAAATGTTATTTAATAACTTTTCTGCTGGTATTTTAGGATTTATCATGACGATTGCAGGATTTAAAATTTTAGCACCACTTATGAAATTTATTATGCATATTTTATCAG  >  USA300TCH1516_ALE/2273399‑2273579

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: