Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,172,078 T→A 9.2% F43I (TTT→ATT)  leuB → 3‑isopropylmalate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,172,0780TA 9.2% 610.6 / 39.5 195F43I (TTT→ATT) leuB3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (83/94);  new base A (8/10);  total (91/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 209 total aligned reads displayed.

TCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTG  >  USA300TCH1516_ALE/2171989‑2172168
                                                                                          |                                                                                           
tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                            <  1:606141/91‑1 (MQ=255)
tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                            <  2:2796912/91‑1 (MQ=255)
tCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAAtt                                                                                            <  2:1196235/91‑1 (MQ=255)
  ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                                          >  1:2736319/1‑91 (MQ=255)
  ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                                          >  1:2956257/1‑91 (MQ=255)
  ggTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtg                                                                                          <  1:1989475/91‑1 (MQ=255)
      cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                      >  1:1426460/1‑91 (MQ=255)
      cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                      >  1:2092645/1‑91 (MQ=255)
      cAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTtggtgg                                                                                      <  2:3206403/91‑1 (MQ=255)
       aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATATGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                     <  2:2165708/91‑1 (MQ=255)
       aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                     >  2:1250231/1‑91 (MQ=255)
       aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                     >  1:1393923/1‑91 (MQ=255)
       aGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTggtggt                                                                                     <  2:268603/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)
         aaaTTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGc                                                                                   <  1:415889/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)
             tttGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCctct                                                                               >  1:167681/1‑91 (MQ=255)
               tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                                  >  1:3872826/1‑86 (MQ=255)
               tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGcc                                                                                  <  2:3872826/86‑1 (MQ=255)
               tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                            <  1:2328963/91‑1 (MQ=255)
               tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                            <  2:341452/91‑1 (MQ=255)
               tGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAt                                                                            <  1:3407095/91‑1 (MQ=255)
                gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                           >  2:3795928/1‑91 (MQ=255)
                gAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTAtt                                                                           >  1:3642801/1‑91 (MQ=255)
                  aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                         >  1:3843760/1‑91 (MQ=255)
                  aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                         <  1:469949/91‑1 (MQ=255)
                  aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                         <  2:2548430/91‑1 (MQ=255)
                  aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                         <  2:3557578/91‑1 (MQ=255)
                  aCGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGa                                                                         <  1:1628168/91‑1 (MQ=255)
                   cGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGAt                                                                        >  1:3730179/1‑91 (MQ=255)
                    ggATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATa                                                                       >  1:3178538/1‑91 (MQ=255)
                     gATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATAc                                                                      >  2:2647293/1‑91 (MQ=255)
                       tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                    <  1:177832/91‑1 (MQ=255)
                       tctcTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAt                                                                    <  2:1658039/91‑1 (MQ=255)
                        ctctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACAtt                                                                   >  1:1845637/1‑91 (MQ=255)
                         tctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTc                                                                  <  2:1288699/91‑1 (MQ=255)
                          ctATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCg                                                                 <  1:3229893/91‑1 (MQ=255)
                           tATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCgg                                                                <  2:2483418/91‑1 (MQ=255)
                             tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCg                                                              <  2:3168831/91‑1 (MQ=255)
                             tCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCg                                                              <  1:709014/91‑1 (MQ=255)
                               aTTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAg                                                            <  1:1212964/91‑1 (MQ=255)
                                ttGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGc                                                           <  1:570742/91‑1 (MQ=255)
                                 tGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGcc                                                          <  1:3285473/91‑1 (MQ=255)
                                    ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                                       >  2:3041442/1‑91 (MQ=255)
                                    ttGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCttt                                                       >  2:3824802/1‑91 (MQ=255)
                                     tGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTa                                                      >  2:1826202/1‑91 (MQ=255)
                                      gacaTTATAAGTAATAAATATAACTATAATTATCCAATAGAGCACGACTAATTAGGTGTTGCCACTATTGATACATCCGACGAGCCTTTat                                                     >  2:3569570/4‑90 (MQ=255)
                                       aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                    >  2:2362860/1‑91 (MQ=255)
                                       aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                    <  1:3784153/91‑1 (MQ=255)
                                       aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                    >  2:1315295/1‑91 (MQ=255)
                                       aaaTTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAAc                                                    <  2:1922229/91‑1 (MQ=255)
                                           tataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgag                                                >  2:1652817/1‑91 (MQ=255)
                                            ataAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTgaga                                               <  2:1362261/91‑1 (MQ=255)
                                              aaGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGaaa                                             >  1:1037347/1‑91 (MQ=255)
                                                 taataaATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAAcc                                          >  2:1170059/1‑91 (MQ=255)
                                                  aataaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                         >  2:648981/1‑91 (MQ=255)
                                                  aagaaagaTAACGTTAATTAGCAAATAGAGCAGCACGAATTGGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCt                                         <  2:3517550/88‑1 (MQ=255)
                                                    taaataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTATAACTGACAAAACCTTa                                       >  2:1521744/1‑91 (MQ=255)
                                                      aataTAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTaaa                                     <  2:921312/91‑1 (MQ=255)
                                                       atatAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAAt                                    >  1:1965501/1‑91 (MQ=255)
                                                        tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                   <  2:1810063/91‑1 (MQ=255)
                                                        tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                   <  2:1092534/91‑1 (MQ=255)
                                                        tataACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATg                                   >  2:1953843/1‑91 (MQ=255)
                                                          taACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCg                                 <  1:2263519/91‑1 (MQ=255)
                                                           aaCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgt                                <  1:2892455/91‑1 (MQ=255)
                                                            aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtg                               <  1:2439036/91‑1 (MQ=255)
                                                            aCTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtg                               >  2:2913077/1‑91 (MQ=255)
                                                             cTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCgtgt                              >  2:3848705/1‑91 (MQ=255)
                                                                    tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                       <  2:1719670/91‑1 (MQ=255)
                                                                    tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                       <  1:1739600/91‑1 (MQ=255)
                                                                    tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                       >  1:1264763/1‑91 (MQ=255)
                                                                    tATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAagag                       <  1:3659803/91‑1 (MQ=255)
                                                                     aTCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                      <  2:3301534/91‑1 (MQ=255)
                                                                     aTCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGc                      >  1:762036/1‑91 (MQ=255)
                                                                      tCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCa                     <  1:2963808/91‑1 (MQ=255)
                                                                       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)
                                                                             gagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa              <  1:1397593/91‑1 (MQ=255)
                                                                             gagCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTa              <  2:3235224/91‑1 (MQ=255)
                                                                              agCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAt             >  1:2969326/1‑91 (MQ=255)
                                                                                gCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtt            <  2:1393923/91‑1 (MQ=255)
                                                                                 caccacCAAATTCGTGGTGCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt          >  1:2807866/1‑91 (MQ=38)
                                                                                  accacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTAtttt          >  1:1812971/1‑91 (MQ=255)
                                                                                    cacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTAc        <  1:3831259/91‑1 (MQ=255)
                                                                                    cacGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTAc        <  1:988619/91‑1 (MQ=255)
                                                                                     acGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCATCTGCTCTTTTAca       <  2:1627641/91‑2 (MQ=38)
                                                                                      cGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAATAGCATCTGCTCTTTTAcac      <  2:2807429/91‑3 (MQ=38)
                                                                                          tttGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg  >  2:143087/1‑91 (MQ=255)
                                                                                          tttGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTg  <  1:3512749/91‑1 (MQ=255)
                                                                                          |                                                                                           
TCGGTCCAGAAATTTTGAACGGATCTCTATCATTGCTTGAAATTATAAGTAATAAATATAACTTTAATTATCAAATAGAGCACCACGAATTTGGTGGTGCCTCTATTGATACATTCGGCGAGCCTTTAACTGAGAAAACCTTAAATGCGTGTAAAAGAGCAGATGCTATTTTACTGGGTG  >  USA300TCH1516_ALE/2171989‑2172168

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: