Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,067:1 +G 6.7% coding (116/1047 nt) leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0671.G 6.7% 846.6 / 53.2 195coding (116/1047 nt)leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (83/99);  new base G (4/9);  total (87/108)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 204 total aligned reads displayed.

CTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAT  >  USA300TCH1516_ALE/2171977‑2172157
                                                                                           |                                                                                           
cTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATaga                                                                                              >  1:983884/1‑91 (MQ=255)
  ggTGATGGAATCCGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGc                                                                                           >  1:2973983/1‑91 (MQ=255)
   gTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGca                                                                                          >  2:673699/1‑91 (MQ=255)
     gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcacc                                                                                        >  1:1719670/1‑91 (MQ=255)
     gATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcacc                                                                                        <  2:1351150/91‑1 (MQ=255)
      aTGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcacca                                                                                       <  2:2868725/91‑1 (MQ=255)
      aTGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcacca                                                                                       <  2:1027863/91‑1 (MQ=255)
       tGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcaccac                                                                                      >  2:2724606/1‑91 (MQ=255)
       tGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGcaccac                                                                                      <  1:913920/91‑1 (MQ=255)
        ggAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACg                                                                                     <  1:547076/91‑1 (MQ=255)
          aaTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGaa                                                                                   <  2:1968373/91‑1 (MQ=255)
          aaTCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGaa                                                                                   <  1:227154/91‑1 (MQ=255)
            tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                 <  2:2796912/91‑1 (MQ=255)
            tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                 <  2:1196235/91‑1 (MQ=255)
            tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                 <  1:606141/91‑1 (MQ=255)
              ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                               >  1:2736319/1‑91 (MQ=255)
              ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                               <  1:1988434/91‑1 (MQ=255)
                  cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                           >  1:1426536/1‑91 (MQ=255)
                  cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                           >  1:2092645/1‑91 (MQ=255)
                   aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                          >  1:1394763/1‑91 (MQ=255)
                   aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                          >  1:1271953/1‑91 (MQ=255)
                   aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                          >  2:1250231/1‑91 (MQ=255)
                   aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                          <  1:1788269/91‑1 (MQ=255)
                    gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGgtggtg                                                                         >  2:2739348/1‑91 (MQ=255)
                    gAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGgtggtg                                                                         <  1:2282282/91‑1 (MQ=255)
                     aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGc                                                                        <  2:2146839/91‑1 (MQ=255)
                     aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACCAAATTCGTGGTGc                                                                        <  1:2495925/91‑1 (MQ=255)
                      aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATGGAGCACCACGAATTTGGTGGTGcc                                                                       >  2:746775/1‑91 (MQ=255)
                      aaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGc                                                                        <  1:3478555/91‑1 (MQ=37)
                       aTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCct                                                                      <  2:847230/91‑1 (MQ=255)
                        ttttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCctc                                                                     >  1:2401683/1‑91 (MQ=255)
                         tttGATCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGctct                                                                     <  2:3579061/91‑4 (MQ=25)
                         tttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGGCACCACCAAATTCGTGGTGctct                                                                     <  2:2345224/91‑4 (MQ=25)
                         tttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCctct                                                                    >  1:167681/1‑91 (MQ=255)
                           tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                       >  1:3872826/1‑86 (MQ=255)
                           tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                       <  2:3872826/86‑1 (MQ=255)
                           tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                 <  1:3407095/91‑1 (MQ=255)
                           tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                 <  1:2328963/91‑1 (MQ=255)
                           tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                 <  2:341452/91‑1 (MQ=255)
                            gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                >  1:3642801/1‑91 (MQ=255)
                            gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                >  2:3795928/1‑91 (MQ=255)
                              aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                              >  1:3842827/1‑91 (MQ=255)
                              aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                              <  1:469949/91‑1 (MQ=255)
                              aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                              <  2:3557578/91‑1 (MQ=255)
                               cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGAt                                                             >  1:3730179/1‑91 (MQ=255)
                               cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGAt                                                             <  1:1422822/91‑1 (MQ=255)
                                ggATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATa                                                            >  1:3178538/1‑91 (MQ=255)
                                   tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                         <  2:1658039/91‑1 (MQ=255)
                                    ctctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAtt                                                        >  1:1845637/1‑91 (MQ=255)
                                      ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCg                                                      >  2:3350654/1‑91 (MQ=255)
                                      ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCg                                                      <  1:3229893/91‑1 (MQ=255)
                                         tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCg                                                   <  2:3169117/91‑1 (MQ=255)
                                         tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCg                                                   <  1:709014/91‑1 (MQ=255)
                                           aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                 <  1:1212964/91‑1 (MQ=255)
                                             tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                               <  2:2529400/91‑1 (MQ=255)
                                             tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                               <  1:3285473/91‑1 (MQ=255)
                                               cTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCtt                                             >  2:2673222/1‑91 (MQ=255)
                                               cTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCtt                                             >  2:1437323/1‑91 (MQ=255)
                                                ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                            >  2:3824802/1‑91 (MQ=255)
                                                ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                            >  2:3041442/1‑91 (MQ=255)
                                                 tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTa                                           >  2:1826202/1‑91 (MQ=255)
                                                  gacaTTATAAGTAATAAATATAACTATAATTATCCAATAGAGCACGACTAATTAGGTGTTGCCACTATTGATACATCCGACGAGCCTTTat                                          >  2:3569570/4‑90 (MQ=255)
                                                   aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                         <  1:3784153/91‑1 (MQ=255)
                                                       tataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                     >  2:1652817/1‑91 (MQ=255)
                                                       tataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                     >  2:1592671/1‑91 (MQ=255)
                                                        ataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgaga                                    <  2:1362289/91‑1 (MQ=255)
                                                          aaGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaa                                  >  1:1037347/1‑91 (MQ=255)
                                                             taataaATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                               >  1:1644485/1‑91 (MQ=255)
                                                             taataaATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                               >  2:1170059/1‑91 (MQ=255)
                                                              aataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                              <  1:573236/91‑1 (MQ=255)
                                                              aataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                              >  2:648981/1‑91 (MQ=255)
                                                              aagaaagaTAACGTTAATTAGCAAATAGAGCAGCACGAATTGGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                              <  2:3517550/88‑1 (MQ=255)
                                                               ataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCtt                             <  1:2918628/91‑1 (MQ=255)
                                                                  aataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaaa                          <  2:921312/91‑1 (MQ=255)
                                                                    tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                        <  1:2770108/91‑1 (MQ=255)
                                                                    tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                        >  2:1953319/1‑91 (MQ=255)
                                                                    tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                        <  2:1810063/91‑1 (MQ=255)
                                                                    tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                        <  2:1092534/91‑1 (MQ=255)
                                                                       aaCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgt                     <  2:2580211/91‑1 (MQ=255)
                                                                        aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtg                    <  1:2439036/91‑1 (MQ=255)
                                                                         cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                   >  2:3848705/1‑91 (MQ=255)
                                                                              aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag              <  2:1301621/91‑1 (MQ=255)
                                                                              aTTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAag              >  1:40842/1‑91 (MQ=255)
                                                                                tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag            <  2:1719670/91‑1 (MQ=255)
                                                                                tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag            <  1:3659803/91‑1 (MQ=255)
                                                                                 aTCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGTGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc           >  2:3229273/1‑91 (MQ=255)
                                                                                 aTCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc           >  2:3289946/1‑91 (MQ=255)
                                                                                 aTCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc           >  1:762036/1‑91 (MQ=255)
                                                                                  tCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCa          <  2:2973983/91‑1 (MQ=255)
                                                                                  tCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa          <  1:1890860/91‑1 (MQ=255)
                                                                                  tCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa          >  2:2471090/1‑91 (MQ=255)
                                                                                  tCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa          <  1:2961689/91‑1 (MQ=255)
                                                                                  atcaaTAGAGGCACCACCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa          >  1:1748364/4‑91 (MQ=21)
                                                                                   tcaaTAGAGGCACCACCAAATTCGTGGTGCTCTATTTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa          >  1:2849582/3‑91 (MQ=17)
                                                                                   cAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAACGCGTGTAAAAGAGCAg         >  2:3703046/1‑91 (MQ=255)
                                                                                   cAAATAGAGCACCACGAATTTGGTGGTGCCACTAGTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg         >  2:3702413/1‑91 (MQ=255)
                                                                                   cAAATAGAGCACCACGAATTTGATGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAg         <  2:3121208/91‑1 (MQ=255)
                                                                                    caaTAGAGGCACCCCCAAATTCGTGGCGCTCTATTTGATACATTCGGCGAGCCTTTAAATGAGAAAACCTTAAATGCGTGTAAAAGAGCAg         >  2:868059/2‑91 (MQ=255)
                                                                                     aaTAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGAt       >  1:2346467/1‑91 (MQ=255)
                                                                                          agCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt  >  1:2969314/1‑91 (MQ=255)
                                                                                           |                                                                                           
CTGGTGATGGAATCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAT  >  USA300TCH1516_ALE/2171977‑2172157

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: