Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,328,872 G→C 12.0% F264L (TTC→TTG lacE ← PTS system lactose‑specific EIICB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,328,8720GC12.0% 661.2 / 70.0 216F264L (TTC→TTGlacEPTS system lactose‑specific EIICB component
Reads supporting (aligned to +/- strand):  ref base G (95/95);  new base C (10/16);  total (105/111)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 223 total aligned reads displayed.

TGTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGAATAATGCA  >  USA300TCH1516_ALE/2328782‑2328962
                                                                                           |                                                                                          
tGTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTg                                                                                            >  2:1776042/1‑91 (MQ=255)
tGTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTg                                                                                            >  2:2191118/1‑91 (MQ=255)
tGTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTg                                                                                            <  2:3217352/91‑1 (MQ=255)
 gTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGa                                                                                           <  1:3046060/91‑1 (MQ=255)
 gTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGa                                                                                           <  2:1813518/91‑1 (MQ=255)
 gTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGa                                                                                           >  1:2647810/1‑91 (MQ=255)
 cttGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGa                                                                                           <  2:3861858/90‑1 (MQ=255)
    gCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGt                                                                                        <  2:753638/91‑1 (MQ=255)
     cACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGt                                                                                        <  2:67303/91‑1 (MQ=255)
     cACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGt                                                                                        <  2:519649/91‑1 (MQ=255)
     cACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGt                                                                                        <  1:2311738/91‑1 (MQ=255)
      ccccGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGtt                                                                                       <  2:3563775/90‑1 (MQ=255)
      aCCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGtt                                                                                       <  2:201041/91‑1 (MQ=255)
       cccGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGttg                                                                                      <  1:3651600/91‑2 (MQ=255)
         cGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGttgct                                                                                    <  1:2518466/91‑4 (MQ=255)
         cGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCt                                                                                   <  1:1160837/91‑1 (MQ=255)
         cGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCt                                                                                   <  1:601358/91‑1 (MQ=255)
          gTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCGAACTTCAAGttgctt                                                                                   <  1:2570140/91‑5 (MQ=255)
           tACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTc                                                                                 <  2:3082757/91‑1 (MQ=255)
           tACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTc                                                                                 <  1:3809995/91‑1 (MQ=255)
           tACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTc                                                                                 <  1:189625/91‑1 (MQ=255)
            aCCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCg                                                                                >  1:3722177/1‑91 (MQ=255)
            aCCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCg                                                                                <  2:3701958/91‑1 (MQ=255)
            aCCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCg                                                                                <  2:373437/91‑1 (MQ=255)
             ccTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGa                                                                               <  1:1039344/91‑1 (MQ=255)
              cTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGat                                                                              >  2:1897727/1‑91 (MQ=255)
                 cAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATAtt                                                                           >  2:1221405/1‑91 (MQ=255)
                  aaaTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTc                                                                          >  2:2394906/1‑91 (MQ=255)
                   aaTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCg                                                                         >  1:3163294/1‑91 (MQ=255)
                   aaTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCg                                                                         <  1:1500538/91‑1 (MQ=255)
                   aaTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCg                                                                         >  1:87573/1‑91 (MQ=255)
                    aTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGc                                                                        <  1:2094901/91‑1 (MQ=255)
                    aTGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGGTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGc                                                                        <  2:2942829/91‑1 (MQ=255)
                      gTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCat                                                                      <  2:3342524/91‑1 (MQ=255)
                      gTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCat                                                                      <  2:2263924/91‑1 (MQ=255)
                      gTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCat                                                                      <  1:3770446/91‑1 (MQ=255)
                       tAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCata                                                                     >  1:4920/1‑91 (MQ=255)
                         aCGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATg                                                                   <  2:1679498/91‑1 (MQ=255)
                         aCGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATg                                                                   >  1:2442044/1‑91 (MQ=255)
                          cGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGt                                                                  >  1:989964/1‑91 (MQ=255)
                          cGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGt                                                                  <  1:3120963/91‑1 (MQ=255)
                           gATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTa                                                                 >  1:769559/1‑91 (MQ=255)
                             tagaCATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAAt                                                               <  2:1654028/88‑1 (MQ=255)
                             tAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAAt                                                               <  1:2610853/91‑1 (MQ=255)
                             tAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAAt                                                               >  2:1951004/1‑91 (MQ=255)
                             tAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAAt                                                               <  1:1376521/91‑1 (MQ=255)
                                 cATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAGCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCt                                                           <  2:3317036/91‑1 (MQ=255)
                                 cATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCt                                                           >  2:907708/1‑91 (MQ=255)
                                 cATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTAGAAGAAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCt                                                           <  2:127299/91‑1 (MQ=255)
                                   tttGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGc                                                         >  1:1000587/1‑91 (MQ=255)
                                   tttGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGc                                                         >  1:272788/1‑91 (MQ=255)
                                      gtgtACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCTACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAAt                                                      <  2:2988220/91‑1 (MQ=255)
                                      gtgtACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAAt                                                      <  2:59625/91‑1 (MQ=255)
                                       tgtACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAAtt                                                     >  2:2653154/1‑91 (MQ=255)
                                         tACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGc                                                   <  2:1101559/91‑1 (MQ=255)
                                           ccTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTg                                                 <  1:1564069/91‑1 (MQ=255)
                                               aTGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGctc                                             <  1:1142832/91‑1 (MQ=255)
                                                 gTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTa                                           <  2:810803/91‑1 (MQ=255)
                                                  taataaTTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTAc                                          >  1:2046303/1‑91 (MQ=255)
                                                    ataatTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACaa                                        <  2:2305623/91‑1 (MQ=255)
                                                    ataatTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACaa                                        >  1:3850184/1‑91 (MQ=255)
                                                      aatTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAAtt                                      >  1:1876802/1‑91 (MQ=255)
                                                      aatTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAAtt                                      >  2:2534365/1‑91 (MQ=255)
                                                      aatTTTATCTGCGTGTTCCCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACaact                                      >  2:1260990/1‑89 (MQ=255)
                                                        ttttATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGa                                    >  2:1436361/1‑91 (MQ=255)
                                                            aTCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGa                                >  1:2116665/1‑91 (MQ=255)
                                                             tCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTTGCTCTACAATTGACGGAc                               >  2:3544845/1‑91 (MQ=255)
                                                             tCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGAc                               <  1:2865255/91‑1 (MQ=255)
                                                              cTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCGTATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGAcc                              <  1:3367303/91‑1 (MQ=255)
                                                               tGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCa                             >  1:2340436/1‑91 (MQ=255)
                                                               tGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCa                             <  1:1218587/91‑1 (MQ=255)
                                                                 cGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATg                           >  1:2317646/1‑91 (MQ=255)
                                                                    gtTCTCCAGCTTGAAGCAACTTGGAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAAt                        <  2:732176/91‑1 (MQ=255)
                                                                    gtTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAAt                        >  1:1784971/1‑91 (MQ=255)
                                                                    gtTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAAt                        <  1:2673666/91‑1 (MQ=255)
                                                                      tctcCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATAc                      >  1:1743177/1‑91 (MQ=255)
                                                                       ctcCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATAcc                     <  2:182392/91‑1 (MQ=255)
                                                                         ccAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTa                   >  1:3573161/1‑91 (MQ=255)
                                                                           aGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACaa                                        >  2:1344343/1‑68 (MQ=255)
                                                                           aGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACa                 <  2:3017551/91‑1 (MQ=255)
                                                                            gCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACaa                >  1:3415182/1‑91 (MQ=255)
                                                                             cTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACCATTGACGGACCATGAATTCCTACaaa               >  2:3366563/1‑91 (MQ=255)
                                                                             cTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACaaa               >  1:3433029/1‑91 (MQ=255)
                                                                             cTTGAAGCAACTTGAAGTTCGCTACGATATTCGGATAGGTAATGGCTGCAATTGCTGGCACTACAATTGACGGACCATGAATACCTACaaa               >  1:3095388/1‑91 (MQ=255)
                                                                               tGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAAcc             <  1:2153055/91‑1 (MQ=255)
                                                                                gAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCa            >  1:1611929/1‑91 (MQ=255)
                                                                                gAAGCAACTTGAAGTTCGCCTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAAcaa            >  1:3205492/1‑89 (MQ=255)
                                                                                 aaGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAg           >  2:571044/1‑91 (MQ=255)
                                                                                    cAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaat        <  1:2935319/91‑1 (MQ=255)
                                                                                      aaCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaata       >  1:628619/1‑91 (MQ=255)
                                                                                      aaCTTCAAGTTGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      >  1:3802729/1‑91 (MQ=255)
                                                                                       aCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      <  2:2780441/91‑1 (MQ=255)
                                                                                       aCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      >  2:2030021/1‑91 (MQ=255)
                                                                                       aCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      <  2:1006281/91‑1 (MQ=255)
                                                                                       aCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      <  2:2114941/91‑1 (MQ=255)
                                                                                       aCTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataa      >  1:1931246/1‑91 (MQ=255)
                                                                                       aCTTCAAGTTGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataat     >  2:2075114/1‑91 (MQ=255)
                                                                                       aCTTCAAGTTGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataat     >  2:2791254/1‑91 (MQ=255)
                                                                                       aCTTCAAGTTGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGaataat     >  2:2076131/1‑91 (MQ=255)
                                                                                           gAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTCCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGAATAATGCa  >  2:2866798/1‑91 (MQ=255)
                                                                                           |                                                                                          
TGTTGCACCCGTACCTCCAAATGTAACGATAAACATTTGTGTACCTGATGTAATAATTTTATCTGCGTGTTCTCCAGCTTGAAGCAACTTGAAGTTCGCTTCGATATTCGCATATGTAATGGCTGCAATTGCTGGCTCTACAATTGACGGACCATGAATACCTACAAACCAGAATAATGCA  >  USA300TCH1516_ALE/2328782‑2328962

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: