Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 702,560 T→A 7.4% I118N (ATT→AAT)  nhaK_1 → Sodium, potassium, lithium and rubidium/H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE702,5600TA 7.4% 808.2 / 41.4 244I118N (ATT→AAT) nhaK_1Sodium, potassium, lithium and rubidium/H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (109/115);  new base A (15/3);  total (124/120)
Fisher's exact test for biased strand distribution p-value = 5.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 251 total aligned reads displayed.

CATTAGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTCTTGA  >  USA300TCH1516_ALE/702470‑702650
                                                                                          |                                                                                          
cATTAGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAAt                                                                                            >  2:2211308/1‑91 (MQ=255)
 aTTAGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAAtt                                                                                           >  1:3219817/1‑91 (MQ=255)
  ttAGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATtg                                                                                          <  2:1248618/91‑1 (MQ=255)
    aGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATtgct                                                                                        >  2:1237078/1‑91 (MQ=255)
     ggATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATtgctg                                                                                       >  1:1267909/1‑91 (MQ=255)
     ggATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATtgctg                                                                                       >  2:270586/1‑91 (MQ=255)
     ggATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCAATTGCAAAtgctg                                                                                       <  1:2517527/91‑1 (MQ=255)
      gATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATtgctgc                                                                                      <  2:1155387/91‑1 (MQ=255)
        ttAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCa                                                                                    <  2:3061756/91‑1 (MQ=255)
        ttAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCa                                                                                    <  1:2337536/91‑1 (MQ=255)
        ttAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCa                                                                                    >  2:3783153/1‑91 (MQ=255)
         tAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCAt                                                                                   <  2:1844626/91‑1 (MQ=255)
            tCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCt                                                                                <  1:1705167/91‑1 (MQ=255)
             cATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCtt                                                                               <  2:697624/91‑1 (MQ=255)
              aTTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCttt                                                                              <  1:3609952/91‑1 (MQ=255)
               ttACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTg                                                                             >  1:1463021/1‑91 (MQ=255)
                 actactGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTc                                                                           >  1:2047114/1‑91 (MQ=255)
                 actactGTGATAGGTGCAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTc                                                                           <  2:3038816/91‑1 (MQ=255)
                  ctactGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATGGCTGCCATTCTTTGTcc                                                                          <  2:601019/91‑1 (MQ=255)
                    actGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTa                                                                        >  2:1061057/1‑91 (MQ=255)
                    actGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTa                                                                        >  2:1066886/1‑91 (MQ=255)
                      tgtgATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACt                                                                      <  1:1192754/91‑1 (MQ=255)
                      tgtgATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACt                                                                      >  2:1136718/1‑91 (MQ=255)
                        tgATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGa                                                                    >  2:356801/1‑91 (MQ=255)
                          aTAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATg                                                                  >  1:853364/1‑91 (MQ=255)
                           tAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTATCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGc                                                                 >  2:2299258/1‑91 (MQ=255)
                           tAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGc                                                                 <  1:935012/91‑1 (MQ=255)
                           tAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTCGAGCAGCATTTGCAATTGCTGCCATTCTTCGTCCTACTGATGc                                                                 >  1:3956517/1‑91 (MQ=255)
                            aGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCa                                                                <  2:2000776/91‑1 (MQ=255)
                            aGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCa                                                                >  1:1312427/1‑91 (MQ=255)
                              gtgtAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGt                                                              >  1:1026154/1‑91 (MQ=255)
                               tgtAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTa                                                             >  2:254107/1‑91 (MQ=255)
                               tgtAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTa                                                             >  1:1945562/1‑91 (MQ=255)
                                gtAGGTTTATTGATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAg                                                            <  1:3580012/91‑1 (MQ=255)
                                 tAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGc                                                           >  2:2868553/1‑91 (MQ=255)
                                   ggTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAg                                                         <  1:1671101/91‑1 (MQ=255)
                                    gTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGt                                                        <  2:1604805/91‑1 (MQ=255)
                                      ttatttatTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTgc                                                      <  1:153165/91‑1 (MQ=255)
                                        atttattCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTgcaa                                                    >  1:3931288/1‑91 (MQ=255)
                                        atttattCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTgcaa                                                    <  1:2326774/91‑1 (MQ=255)
                                        atttattCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTgcaa                                                    >  1:680319/1‑91 (MQ=255)
                                        atttattCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACGGATGCAGTAGCCGTgcaa                                                    >  2:3430009/1‑91 (MQ=255)
                                            attCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTgcaagcaa                                                <  1:3364416/91‑1 (MQ=255)
                                             ttCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAAt                                               >  2:2984035/1‑91 (MQ=255)
                                                aTTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCAc                                            <  2:3730422/91‑1 (MQ=255)
                                                 ttGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACt                                           >  1:177126/1‑91 (MQ=255)
                                                   ggATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaa                                         <  1:2624106/91‑1 (MQ=255)
                                                      tttGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaagg                                      <  1:98383/91‑1 (MQ=255)
                                                       ttGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaagga                                     <  1:1404744/91‑1 (MQ=255)
                                                        tGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaaggaa                                    >  2:2464427/1‑91 (MQ=255)
                                                        tGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaaggaa                                    >  1:1719725/1‑91 (MQ=255)
                                                        tGGCCAGATTTACCTATTGGAGCAGCATTTGCAAGTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaaggaa                                    <  2:1488185/91‑1 (MQ=255)
                                                          gCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTaaaggaaag                                  >  1:57736/1‑91 (MQ=255)
                                                             aGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTc                               >  2:899415/1‑91 (MQ=255)
                                                             aGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTc                               >  1:464596/1‑91 (MQ=255)
                                                              gATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCt                              <  1:3193626/91‑1 (MQ=255)
                                                                tttACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTg                            <  1:11558/91‑1 (MQ=255)
                                                                 ttACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGc                           <  2:1080883/91‑1 (MQ=255)
                                                                 ttACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGc                           <  1:2430317/91‑1 (MQ=255)
                                                                  tACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGcc                          <  2:1734549/91‑1 (MQ=255)
                                                                  tACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGcc                          <  1:2801907/91‑1 (MQ=255)
                                                                   aCCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCa                         >  1:2145671/1‑91 (MQ=255)
                                                                   aCCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCa                         >  2:3173391/1‑91 (MQ=255)
                                                                   aCCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCa                         >  1:820342/1‑91 (MQ=255)
                                                                    ccTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaa                        >  2:3348364/1‑91 (MQ=255)
                                                                    ccTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaa                        >  2:3075924/1‑91 (MQ=255)
                                                                    ccTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaa                        <  1:2842030/91‑1 (MQ=255)
                                                                     cTATTGGAGCAGCTTTTGCAATTGCTGCCCTTCTTTGTCCTACTGATGCAGTATCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaaa                       >  2:243291/1‑91 (MQ=255)
                                                                     cTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaaa                       >  2:2798599/1‑91 (MQ=255)
                                                                      tATTGGAGTAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCaaaa                      <  1:2946888/91‑1 (MQ=255)
                                                                      tATCGGAGCAGAATATGCAATTGCTGCCATTACTAGTCTTACTGATGCAGCCGCAGAACAACCAATCCCTAAAGGAAAGGTCCTGCCaaaa                      >  2:404918/1‑91 (MQ=255)
                                                                          ggAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAg                  <  1:904203/91‑1 (MQ=255)
                                                                          ggAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAg                  >  1:3772262/1‑91 (MQ=255)
                                                                             gcagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCaa               >  2:3027445/1‑91 (MQ=255)
                                                                             gcagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCaa               >  1:508431/1‑91 (MQ=255)
                                                                             gcagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCaa               >  2:3396483/1‑91 (MQ=255)
                                                                              cagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAAt              <  1:547458/91‑1 (MQ=255)
                                                                              cagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAAt              <  2:1933446/91‑1 (MQ=255)
                                                                              cagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAAt              >  1:1628365/1‑91 (MQ=255)
                                                                              cagcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAACAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAAt              >  1:2655729/1‑91 (MQ=255)
                                                                               agcaTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATg             <  2:810316/91‑1 (MQ=255)
                                                                                gcaattGCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGa            >  1:2063453/5‑91 (MQ=255)
                                                                                 caTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGAc           >  2:1446206/1‑91 (MQ=255)
                                                                                  aTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACa          <  1:413236/91‑1 (MQ=255)
                                                                                  aTTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACa          <  2:1742319/91‑1 (MQ=255)
                                                                                   attGCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACaa         >  2:3517157/2‑91 (MQ=255)
                                                                                    ttGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        <  2:1694362/91‑1 (MQ=255)
                                                                                    ttGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        <  2:3254098/91‑1 (MQ=255)
                                                                                    ttGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        >  2:492680/1‑91 (MQ=255)
                                                                                    ttGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        >  2:1951515/1‑91 (MQ=255)
                                                                                    ttGCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        >  2:1528762/1‑91 (MQ=255)
                                                                                    ttGCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        >  1:1057360/1‑91 (MQ=255)
                                                                                    ttGCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAt        >  1:256867/1‑91 (MQ=255)
                                                                                     tGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAtt       >  1:1918737/1‑91 (MQ=255)
                                                                                     tGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAtt       <  1:2175232/91‑1 (MQ=255)
                                                                                     tGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAAtt       <  2:2112185/91‑1 (MQ=255)
                                                                                      gCAAATGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTc      >  2:1412119/1‑91 (MQ=255)
                                                                                         aTTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTCTTg   <  2:3102319/91‑1 (MQ=255)
                                                                                          ttGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTCTTGa  >  2:387367/1‑91 (MQ=255)
                                                                                          ttGCTGCCATTCTTTGTCCTACTGATGCACTAGCAGTGCAAGCAATTACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTCTTGa  <  1:3167982/91‑1 (MQ=255)
                                                                                          |                                                                                          
CATTAGGATTAGTCATTACTACTGTGATAGGTGTAGGTTTATTTATTCATTGGATTTGGCCAGATTTACCTATTGGAGCAGCATTTGCAATTGCTGCCATTCTTTGTCCTACTGATGCAGTAGCAGTGCAAGCAATCACTAAAGGAAAGGTCTTGCCAAAAGGAGCAATGACAATTCTTGA  >  USA300TCH1516_ALE/702470‑702650

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: