Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,920:1 +T 6.1% intergenic (+343/‑47) 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,9201.T 6.1% 484.3 / 23.3 115intergenic (+343/‑47)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (41/67);  new base T (4/3);  total (45/70)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 123 total aligned reads displayed.

TGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGT  >  USA300TCH1516_ALE/576829‑577009
                                                                                            |                                                                                         
tGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTa                                                                                            <  1:812603/91‑1 (MQ=255)
 gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                                          <  1:3953726/91‑1 (MQ=255)
  aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAtt                                                                                         <  1:1420747/91‑2 (MQ=255)
   tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATg                                                                                         <  1:928357/91‑1 (MQ=255)
    aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAttgt                                                                                       <  1:794550/91‑4 (MQ=37)
    aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAttgt                                                                                       <  1:796663/91‑4 (MQ=37)
    aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAttgt                                                                                       <  1:796894/91‑4 (MQ=37)
      tttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTa                                                                                      >  1:467994/1‑91 (MQ=255)
        tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                    >  1:251085/1‑91 (MQ=255)
        tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                    <  1:3898590/91‑1 (MQ=255)
        tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                                    <  1:530225/91‑2 (MQ=25)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   >  1:582646/1‑91 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   >  1:583843/1‑91 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   <  1:1877968/91‑1 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   >  1:1448250/1‑91 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   <  1:3017543/91‑1 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   <  1:3018120/91‑1 (MQ=255)
         ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                   >  1:313718/1‑91 (MQ=255)
          gAAGCGCTAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTTTGTTAgtgg                                                                                  <  1:2080180/91‑1 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  <  1:3943149/91‑1 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  <  1:3930069/91‑1 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  <  1:3943186/91‑1 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  >  1:1070593/1‑91 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  >  1:61863/1‑91 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                  >  1:1070572/1‑91 (MQ=255)
          gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTAGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgagg                                                                                  <  1:3313481/91‑4 (MQ=21)
           aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                 <  1:853616/91‑1 (MQ=255)
           aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                 <  1:514336/91‑1 (MQ=255)
            aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                >  1:941629/1‑91 (MQ=255)
            aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                >  1:1256827/1‑91 (MQ=255)
            aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                <  1:1737544/91‑1 (MQ=255)
              cgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCt                                                                              <  1:659582/91‑1 (MQ=255)
               gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                             <  1:3305641/91‑1 (MQ=255)
               gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACCTGTATGTTAGTGGTGCtt                                                                             <  1:1203474/91‑1 (MQ=255)
                ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCttt                                                                            >  1:3405361/1‑91 (MQ=255)
                 cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttt                                                                           <  1:1261554/91‑1 (MQ=255)
                 cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttt                                                                           <  1:509659/91‑1 (MQ=255)
                   gCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttttt                                                                         >  1:939927/1‑91 (MQ=255)
                     aaTATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGa                                                                       >  1:3172128/1‑91 (MQ=255)
                      atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                      >  1:905590/1‑91 (MQ=255)
                      atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                                      <  1:2905846/91‑1 (MQ=255)
                        attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                    <  1:2648095/91‑1 (MQ=255)
                        attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                    <  1:188949/91‑1 (MQ=255)
                        attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                    <  1:342201/91‑1 (MQ=255)
                         ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAtt                                                                   <  1:2096710/91‑1 (MQ=255)
                         ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAtt                                                                   <  1:2440183/91‑1 (MQ=255)
                              aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAtt                                                              <  1:2927493/91‑1 (MQ=255)
                               aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                             >  1:3278012/1‑91 (MQ=255)
                                tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                            <  1:3375016/91‑1 (MQ=255)
                                tGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTc                                                            <  1:1025110/91‑1 (MQ=255)
                                 gTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCg                                                           >  1:856035/1‑91 (MQ=255)
                                   tATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCa                                                         >  1:2621816/1‑91 (MQ=255)
                                   tATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCa                                                         >  1:2622758/1‑91 (MQ=255)
                                     ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       >  1:939040/1‑91 (MQ=255)
                                     ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                       <  1:1360387/91‑1 (MQ=255)
                                        gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                    >  1:1662852/1‑91 (MQ=255)
                                          tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACtt                                                  <  1:198474/91‑1 (MQ=255)
                                           aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                                 >  1:1863053/1‑91 (MQ=255)
                                            aaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCa                                                <  1:3312540/91‑1 (MQ=255)
                                               tCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATTACAACTTGTATGTTAGTGGTGCTTCTTTGATATTTAATTTCGCAGGTACTTCAAtt                                             >  1:3306414/1‑91 (MQ=255)
                                                  gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                          >  1:3887480/1‑91 (MQ=255)
                                                  gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                          >  1:3887525/1‑91 (MQ=255)
                                                  gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                                          >  1:809026/1‑91 (MQ=255)
                                                   agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                         >  1:3064817/1‑91 (MQ=255)
                                                   agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGt                                         <  1:3164412/91‑1 (MQ=255)
                                                    gTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTa                                        <  1:3909779/91‑1 (MQ=255)
                                                    gTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTa                                        >  1:857514/1‑91 (MQ=255)
                                                      ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                      <  1:3498737/91‑1 (MQ=255)
                                                       gAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAttt                                     >  1:464855/1‑91 (MQ=255)
                                                              gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                              <  1:1020825/91‑1 (MQ=255)
                                                              gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                              <  1:1149128/91‑1 (MQ=255)
                                                              gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                              <  1:2718008/91‑1 (MQ=255)
                                                                ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                            <  1:2035057/91‑1 (MQ=255)
                                                                 aaaGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                           >  1:3774997/2‑91 (MQ=18)
                                                                   aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                         >  1:2442233/1‑91 (MQ=255)
                                                                   aGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAt                         <  1:4065405/91‑1 (MQ=255)
                                                                    gcgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATc                        <  1:765379/91‑1 (MQ=255)
                                                                     cgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCa                       >  1:2194348/1‑91 (MQ=255)
                                                                      gcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAt                      >  1:2937132/1‑91 (MQ=255)
                                                                      gcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAt                      >  1:2487034/1‑91 (MQ=255)
                                                                      gcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAt                      >  1:4079752/1‑91 (MQ=255)
                                                                          ctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag                  <  1:2216500/91‑1 (MQ=255)
                                                                           tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg                 <  1:1283693/91‑1 (MQ=255)
                                                                            aaCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagga                <  1:719043/91‑1 (MQ=255)
                                                                             aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag               <  1:3214688/91‑1 (MQ=255)
                                                                               aaTGCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAGCATTGAggagga             <  1:3207381/91‑1 (MQ=255)
                                                                               aaTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAggagga             <  1:3691281/91‑1 (MQ=255)
                                                                               aaTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAAGCATTGAggagga             <  1:3207581/91‑1 (MQ=255)
                                                                                 tACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa            >  1:2310622/1‑91 (MQ=255)
                                                                                 tACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa            >  1:2310875/1‑91 (MQ=255)
                                                                                 tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa           <  1:698411/91‑1 (MQ=255)
                                                                                  aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          <  1:2224792/91‑1 (MQ=255)
                                                                                  aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          <  1:2642359/91‑1 (MQ=255)
                                                                                  aCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAAt          >  1:4093997/1‑91 (MQ=255)
                                                                                  aCAACTTGTATGTTAGTGGTGCTTTTTTGATATGTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaacc          <  1:3399690/91‑3 (MQ=255)
                                                                                     aCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGAt       <  1:3675995/91‑1 (MQ=255)
                                                                                      cTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATc      <  1:1795207/91‑1 (MQ=255)
                                                                                        tGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCtt    <  1:675064/91‑1 (MQ=255)
                                                                                          tATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGt  <  1:394452/91‑1 (MQ=255)
                                                                                          tATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGt  <  1:445168/91‑1 (MQ=255)
                                                                                            |                                                                                         
TGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGT  >  USA300TCH1516_ALE/576829‑577009

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