Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,550,126 G→A 7.0% I162I (ATC→ATT USA300HOU_RS07580 ← Glucosaminate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,550,1260GA 7.0% 254.8 / 6.4 88I162I (ATC→ATTUSA300HOU_RS07580Glucosaminate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (2/4);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

TGATTTACTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACATTCATA  >  USA300TCH1516_ALE/1549986‑1550266
                                                                                                                                            |                                                                                                                                            
tGATTTACTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCg                                                                                                                                              >  1:610250/1‑141 (MQ=255)
tGATTTACTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCg                                                                                                                                              >  1:103635/1‑141 (MQ=255)
     tACTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAAt                                                                                                                                         <  1:476522/141‑1 (MQ=255)
      aCTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCTATAAtt                                                                                                                                        >  4:16444/1‑141 (MQ=255)
      aCTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAAtt                                                                                                                                        <  2:443884/141‑1 (MQ=255)
      aCTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAAtt                                                                                                                                        >  4:332334/1‑141 (MQ=255)
       cTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCaataataa                                                                                                                                       >  2:56820/1‑139 (MQ=255)
        tAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTac                                                                                                                                      <  2:940674/141‑1 (MQ=255)
         aaTGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTaca                                                                                                                                     >  3:44317/1‑141 (MQ=255)
          aTGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTacaa                                                                                                                                    >  2:166131/1‑141 (MQ=255)
          aTGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTacaa                                                                                                                                    >  2:1000402/1‑141 (MQ=255)
          aTGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTacaa                                                                                                                                    >  2:734366/1‑141 (MQ=255)
               gtgAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCt                                                                                                                               >  4:312129/1‑141 (MQ=255)
                 gAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCttt                                                                                                                             >  4:233876/1‑140 (MQ=255)
                   aaTTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGAt                                                                                                                           >  3:440579/1‑141 (MQ=255)
                   aaTTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAAACGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGACTTCTTACCACCGATAATTACAACACCTTGAt                                                                                                                           >  2:686546/1‑141 (MQ=255)
                    aTTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATc                                                                                                                          <  1:757284/141‑1 (MQ=255)
                      ttGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCaa                                                                                                                        <  2:668715/141‑1 (MQ=255)
                      ttGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATATTGATCaa                                                                                                                        <  1:389440/141‑1 (MQ=255)
                            gTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACg                                                                                                                  <  2:324839/141‑1 (MQ=255)
                              aTCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGa                                                                                                                >  4:164297/1‑141 (MQ=255)
                                   cGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATTTGc                                                                                                           >  4:366238/1‑141 (MQ=255)
                                    ggTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGcc                                                                                                          <  1:843626/141‑1 (MQ=255)
                                       ttAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCctct                                                                                                       >  3:404953/1‑141 (MQ=255)
                                       ttAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCctct                                                                                                       <  1:271016/141‑1 (MQ=255)
                                        tAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCTAAAGCAGCATCGATAGCAGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCtt                                                                                                      >  4:280396/1‑141 (MQ=255)
                                        tAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCtt                                                                                                      <  2:254412/141‑1 (MQ=255)
                                        tAATGGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCtt                                                                                                      <  2:431371/141‑1 (MQ=255)
                                         aaTTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACTGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCttt                                                                                                     <  1:934732/141‑1 (MQ=255)
                                         aaTTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCttt                                                                                                     <  1:1085369/141‑1 (MQ=255)
                                          aTTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTa                                                                                                    >  3:380582/1‑141 (MQ=255)
                                            tGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTaaa                                                                                                  <  2:215594/141‑1 (MQ=255)
                                                  cGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATg                                                                                            <  2:57292/141‑1 (MQ=255)
                                                   gAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATgg                                                                                           <  1:1038428/141‑1 (MQ=255)
                                                    aataatCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGCTAGCCGAATTCTTACCACCGCTAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGa                                                                                          <  1:285136/141‑1 (MQ=255)
                                                    aataatCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGa                                                                                          <  1:749836/141‑1 (MQ=255)
                                                        atCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGGTCAAAATACGGATGTGCCTCTTTAAAATAATGGTacac                                                                                      <  1:1058065/141‑1 (MQ=255)
                                                        atCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAacac                                                                                      <  2:529322/141‑1 (MQ=255)
                                                          ctcCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACtt                                                                                    >  1:398787/1‑141 (MQ=255)
                                                            ccaccaCGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTa                                                                                  <  4:203493/141‑1 (MQ=255)
                                                            ccaccaCGATATAGAACCGGCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCACCACCGATGGCCGAATGCTGATCACCGATAATTACAACATGTTGATCAAAATACCGATGGTCCTCTTTAAAATAATGGAACACTGTa                                                                                  >  1:813049/1‑141 (MQ=255)
                                                               ccacGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGt                                                                               >  2:110836/1‑141 (MQ=255)
                                                                cacGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTa                                                                              <  2:210185/141‑1 (MQ=255)
                                                                  cGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTa                                                                                                    >  2:391909/1‑117 (MQ=255)
                                                                  cGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTa                                                                                                    <  1:391909/117‑1 (MQ=255)
                                                                   gATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAAt                                                                           >  3:99259/1‑141 (MQ=255)
                                                                    atatAGAGCCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATc                                                                          >  3:115076/1‑141 (MQ=255)
                                                                    atatAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGCACACTTTAGGTAAATc                                                                          >  2:70257/1‑141 (MQ=255)
                                                                    atatAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATc                                                                          <  4:264401/141‑1 (MQ=255)
                                                                         gAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCAc                                                                     >  3:445998/1‑141 (MQ=255)
                                                                         gAACCGTCACGTTAGCACCAACTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTAGATCAAAATACGGATGTGCCTCTTTAAAATAAATAAACACCTTAGGTAAAtacaaac                                                                     >  4:198199/1‑135 (MQ=255)
                                                                           aCCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCt                                                                   >  1:683253/1‑141 (MQ=255)
                                                                             cGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTc                                                                 <  1:600118/141‑1 (MQ=255)
                                                                              gTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCa                                                                <  4:380715/141‑1 (MQ=255)
                                                                              gTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCa                                                                >  2:923925/1‑141 (MQ=255)
                                                                               tCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGAAGTGTCTCTTTCATATAATGGAAAACTTTAGGTAAAACCGCAccattcc                                                               >  3:202511/1‑136 (MQ=255)
                                                                                   gTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTCAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTc                                                           >  4:335967/1‑141 (MQ=255)
                                                                                   gTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTc                                                           <  4:419789/141‑1 (MQ=255)
                                                                                   gTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTc                                                           <  1:544133/141‑1 (MQ=255)
                                                                                   gTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTc                                                           >  2:958796/1‑141 (MQ=255)
                                                                                          ccAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGta                                                    <  1:324766/141‑1 (MQ=255)
                                                                                             gCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGtatta                                                 >  1:527516/1‑141 (MQ=255)
                                                                                                ttttCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGc                                              >  4:398175/1‑141 (MQ=255)
                                                                                                ttttCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGc                                              <  1:359750/141‑1 (MQ=255)
                                                                                                  ttCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTg                                            >  2:187542/1‑141 (MQ=255)
                                                                                                    ccAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGAc                                          >  3:176145/1‑141 (MQ=255)
                                                                                                    ccAACTCCAAAGCAGCATCAATAGCCGAATACTTACCACCGATAATTACAACATCTTGATCCAAATACGGATGTGCCTCTCTAAAATAATGGAACACTTTAGGTAACTCTGCACCTTCATCTTCTAATGTATTATGCTGAc                                          <  4:137700/141‑1 (MQ=255)
                                                                                                    ccAACTCAAAAGCAGCATCGATAGCCGAATTCTTACCACCAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCTCCTTCAACTTCTAATGTATTATGCTGAc                                          <  4:259862/141‑1 (MQ=255)
                                                                                                        ttttttagcagcATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGCTGTGCCTCTTTTAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATTCTGACCata                                      >  2:776138/7‑141 (MQ=255)
                                                                                                         tCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCataa                                     <  1:529753/141‑1 (MQ=255)
                                                                                                                  gcaTCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATCCGGATGTGCCTCTTTAAAATAATGGAACACGTAAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATTGCCTGTc                            >  2:549847/1‑141 (MQ=255)
                                                                                                                  gcaTCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTATTTAAAATAATGGAACACTTTAGGTAAATCCGCTCCTTCAACTTCTAATGTATTATGCTGAACataatacaatgta                            >  4:283659/1‑134 (MQ=255)
                                                                                                                      cGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGa                        <  3:325428/141‑1 (MQ=255)
                                                                                                                        aTAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACTTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGAtt                      <  3:377180/141‑1 (MQ=255)
                                                                                                                          aGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATATGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGt                    >  2:681278/1‑141 (MQ=255)
                                                                                                                           gCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGtt                   <  1:40898/141‑1 (MQ=255)
                                                                                                                           gCCGAATTCTTACCACCAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCTCCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGtt                   <  1:56820/141‑1 (MQ=255)
                                                                                                                            ccGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTa                  >  2:682495/1‑141 (MQ=255)
                                                                                                                              gAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTaaa                <  1:601938/141‑1 (MQ=255)
                                                                                                                                aTTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTaaaaa              <  4:209315/141‑1 (MQ=255)
                                                                                                                                     tACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGAc         <  2:633886/141‑1 (MQ=255)
                                                                                                                                     tACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGAc         >  3:221383/1‑141 (MQ=255)
                                                                                                                                      accaccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACa        >  3:296500/1‑141 (MQ=255)
                                                                                                                                       ccaccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTTACCATAATAGCCTGTCGCGATTGTTAAAAATCGACAt       <  1:791543/141‑1 (MQ=255)
                                                                                                                                       ccaccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACAt       <  2:636056/141‑1 (MQ=255)
                                                                                                                                       ccaccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACAt       >  1:141169/1‑141 (MQ=255)
                                                                                                                                        caccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACAtt      <  4:380582/141‑1 (MQ=255)
                                                                                                                                        caccAATAATAACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCTCCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACAtt      >  1:344003/1‑141 (MQ=255)
                                                                                                                                          ccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACATTCa    >  3:266305/1‑141 (MQ=255)
                                                                                                                                          ccGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACATTCa    >  1:174883/1‑141 (MQ=255)
                                                                                                                                           cGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGAGAATGATAAAAATCGACATTCAt   <  2:1016109/141‑1 (MQ=255)
                                                                                                                                            gATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACATTCATa  <  4:350495/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                            
TGATTTACTAATGCTGTGAAATTTGGAAGTATCCACGGTTTAATTGAAGGCGAATAATCTCCACCACGATATAGAACCGTCACGTTAGCACCAGCTTTTTCCAACTCCAAAGCAGCATCGATAGCCGAATTCTTACCACCGATAATTACAACATCTTGATCAAAATACGGATGTGCCTCTTTAAAATAATGGAACACTTTAGGTAAATCCGCACCTTCAACTTCTAATGTATTATGCTGACCATAATAGCCTGTCGCGATTGTTAAAAATCGACATTCATA  >  USA300TCH1516_ALE/1549986‑1550266

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