Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,915 C→G 9.7% intergenic (+338/‑52) gltX → / → cysE Glutamate‑‑tRNA ligase/Serine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE576,9150CG 9.7% 707.5 / 49.2 206intergenic (+338/‑52)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (66/120);  new base G (4/16);  total (70/136)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 207 total aligned reads displayed.

CCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATC  >  USA300TCH1516_ALE/576826‑577005
                                                                                         |                                                                                           
ccTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACt                                                                                            <  1:1532971/91‑1 (MQ=255)
 cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                           <  2:2071840/91‑1 (MQ=255)
 cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                           <  1:433465/91‑1 (MQ=255)
  ttGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTg                                                                                          <  1:1853070/91‑1 (MQ=255)
    gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTa                                                                                        <  2:940041/91‑1 (MQ=255)
    gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                      <  2:3243142/91‑1 (MQ=255)
     aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAtt                                                                                     <  1:2082725/91‑2 (MQ=255)
     aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCACTAACATACAAGTTGTAtt                                                                                     <  1:261899/91‑2 (MQ=255)
      tAATTTGGAAGCGCCAGCACTATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGGGCGGAAGCGCCTCTAACAACACAACGTGTATg                                                                                     >  2:1538929/1‑91 (MQ=255)
       aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTCGTAttgt                                                                                   <  2:743030/91‑4 (MQ=37)
           tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                <  1:922941/91‑1 (MQ=255)
           tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                <  2:2817488/91‑1 (MQ=255)
           tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTCGTATGTTAgt                                                                                >  2:1502656/1‑91 (MQ=255)
           tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                                <  1:3434866/91‑2 (MQ=25)
           tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                                <  1:1300394/91‑2 (MQ=25)
            ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                               <  2:269150/91‑1 (MQ=255)
             gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                              <  2:3616388/91‑1 (MQ=255)
             gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                              <  2:928036/91‑1 (MQ=255)
             gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                              <  2:928130/91‑1 (MQ=255)
             gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                              <  2:956619/91‑1 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             >  2:1247159/1‑91 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             >  1:2507435/1‑91 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             >  2:3521978/1‑91 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             <  2:3534818/91‑1 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                             <  1:16362/91‑1 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCTCTAACAATACAACTTGTATGTTAgtggt                                                                             <  2:2549735/91‑1 (MQ=255)
              aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                             <  1:663743/91‑5 (MQ=21)
               aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                            <  2:800678/91‑1 (MQ=255)
                gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                           >  1:1129613/1‑91 (MQ=255)
                gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                           >  2:1481598/1‑91 (MQ=255)
                gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                           <  1:3251899/91‑1 (MQ=255)
                  gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                         <  1:3565018/91‑1 (MQ=255)
                   ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCttc                                                                        <  1:3742098/91‑2 (MQ=17)
                    cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttt                                                                       <  1:360736/91‑1 (MQ=255)
                     aGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttttt                                                                      >  2:1200458/1‑91 (MQ=255)
                     aGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttttt                                                                      >  1:3273264/1‑91 (MQ=255)
                      gCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttttt                                                                     >  2:1039102/1‑91 (MQ=255)
                      gCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttttt                                                                     <  1:2839272/91‑1 (MQ=255)
                       cAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTg                                                                    <  1:2311950/91‑1 (MQ=255)
                           attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                <  1:992831/91‑1 (MQ=255)
                            ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTCTGTTAGTGGTGCTTTTTTTATAtt                                                               >  2:3671883/1‑91 (MQ=255)
                               ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                            <  2:816269/91‑1 (MQ=255)
                                taATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAt                                                           <  1:699566/91‑1 (MQ=255)
                                 aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTCGTGGTGCTTTTTTGATATTTAAtt                                                          >  2:2837751/1‑91 (MQ=255)
                                 aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAtt                                                          >  2:836717/1‑91 (MQ=255)
                                  aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                         <  2:2212094/91‑1 (MQ=255)
                                  aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                         <  1:3521128/91‑1 (MQ=255)
                                  aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                         <  1:1155878/91‑1 (MQ=255)
                                     ttATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGc                                                      <  2:2893596/91‑1 (MQ=255)
                                     ttATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGc                                                      <  2:1323175/91‑1 (MQ=255)
                                       aTTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAg                                                    <  2:389352/91‑1 (MQ=255)
                                       aTTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAg                                                    >  2:3162847/1‑91 (MQ=255)
                                        ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                   >  1:336849/1‑91 (MQ=255)
                                         tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                                  <  1:998068/91‑1 (MQ=255)
                                          cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTa                                                 <  2:3954379/91‑1 (MQ=255)
                                          cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTa                                                 >  1:3918082/1‑91 (MQ=255)
                                           gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                <  2:3942823/91‑1 (MQ=255)
                                           gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                <  1:1562302/91‑1 (MQ=255)
                                           gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                >  1:2375368/1‑91 (MQ=255)
                                             tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTATTTTCGCAGGTACtt                                              >  2:3416686/1‑91 (MQ=255)
                                             tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACtt                                              <  1:1735459/91‑1 (MQ=255)
                                              aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATTCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                             <  2:1970291/91‑1 (MQ=255)
                                                aTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCaa                                           <  1:395886/91‑1 (MQ=255)
                                                    agagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTAt                                       <  2:3444844/91‑1 (MQ=255)
                                                     gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                      >  1:2987796/1‑91 (MQ=255)
                                                      agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                     >  1:2660418/1‑91 (MQ=255)
                                                         ggAACCGTGCGGAAGCGCGTCTAACAATACAACTTGTATGTTAGTGGTGATTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  2:1268119/91‑1 (MQ=255)
                                                         ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  1:2478123/91‑1 (MQ=255)
                                                         ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  2:1106598/91‑1 (MQ=255)
                                                         ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  2:525134/91‑1 (MQ=255)
                                                         ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                  <  1:3353857/91‑1 (MQ=255)
                                                              cGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                             <  1:955628/91‑1 (MQ=255)
                                                                 gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                          <  2:3120795/91‑1 (MQ=255)
                                                                   ggAAGCTCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                        >  1:3074820/1‑91 (MQ=255)
                                                                   ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                        >  2:855723/1‑91 (MQ=255)
                                                                       gcgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATc                    <  2:2056091/91‑1 (MQ=255)
                                                                        cgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCa                   >  1:276141/1‑91 (MQ=255)
                                                                         gcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAt                  <  2:383674/91‑1 (MQ=255)
                                                                          ccactAACATACAAGTTGTATTGTTAGAGGCGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAtt                 >  2:3739322/4‑91 (MQ=21)
                                                                             ctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTTTTATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag              >  2:2491433/1‑91 (MQ=255)
                                                                             ctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag              <  1:1540411/91‑1 (MQ=255)
                                                                             ctAACAATACAACTTGTATGTTAGTGCTGCTTTTTTTTTTTTTTATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag              >  2:2242023/1‑91 (MQ=255)
                                                                             ctAACATACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag              >  1:571765/1‑91 (MQ=38)
                                                                               aaCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagga            <  1:2134919/91‑1 (MQ=255)
                                                                                 cAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGaggagg          >  1:347914/1‑91 (MQ=255)
                                                                                   aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGTAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa        <  1:3135194/91‑1 (MQ=255)
                                                                                   aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa        <  1:1893737/91‑1 (MQ=255)
                                                                                   aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa        <  1:3066306/91‑1 (MQ=255)
                                                                                    tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa       <  2:485626/91‑1 (MQ=255)
                                                                                    tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa       >  1:287918/1‑91 (MQ=255)
                                                                                    tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa       <  1:1431258/91‑1 (MQ=255)
                                                                                     aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt      <  2:461464/91‑1 (MQ=255)
                                                                                     aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt      >  2:2460014/1‑91 (MQ=255)
                                                                                      cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg     <  2:3000976/91‑1 (MQ=255)
                                                                                      cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg     >  1:3845666/1‑91 (MQ=255)
                                                                                      cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg     <  2:2430961/91‑1 (MQ=255)
                                                                                       aaCTTGTATGTTAGTGGTGCTTTTTTTTTTTTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGa    >  2:3839971/1‑91 (MQ=255)
                                                                                       aaCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGa    <  1:1247046/91‑1 (MQ=255)
                                                                                        aCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGAt   <  1:2102889/91‑1 (MQ=255)
                                                                                        aCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATATCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGAt   <  1:2103297/91‑1 (MQ=255)
                                                                                         |                                                                                           
CCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATC  >  USA300TCH1516_ALE/576826‑577005

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