Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,908 Δ1 bp 10.1% intergenic (+331/‑59) 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,9080A.10.1% 600.2 / 58.6 188intergenic (+331/‑59)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (64/105);  new base . (3/16);  total (67/121)
Fisher's exact test for biased strand distribution p-value = 7.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 193 total aligned reads displayed.

CTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGA  >  USA300TCH1516_ALE/576818‑576998
                                                                                          |                                                                                           
cTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACa                                                                                             <  1:2601166/91‑1 (MQ=255)
cTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCAGAAGCGCCTCTAACa                                                                                             <  2:336009/91‑1 (MQ=255)
  tAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAAt                                                                                           <  2:613042/91‑1 (MQ=255)
  tAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAAt                                                                                           <  1:271217/91‑1 (MQ=255)
   aGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATa                                                                                          <  1:3630828/91‑1 (MQ=255)
    gTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATAc                                                                                         >  1:3830418/1‑91 (MQ=255)
    gTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATAc                                                                                         >  2:2153430/1‑91 (MQ=255)
     tttCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACa                                                                                        >  1:3415269/1‑91 (MQ=255)
      ttCCTTGATAATGTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACaa                                                                                       <  2:394511/91‑1 (MQ=255)
        ccTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACt                                                                                     >  1:1178934/1‑91 (MQ=255)
         cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                    <  1:433465/91‑1 (MQ=255)
         cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                    <  2:1162804/91‑1 (MQ=255)
         cTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACtt                                                                                    <  2:2071840/91‑1 (MQ=255)
          ttGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTg                                                                                   <  1:603751/91‑1 (MQ=255)
            gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTa                                                                                 <  2:2902966/91‑1 (MQ=255)
            gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAt                                                                               <  2:3243142/91‑1 (MQ=255)
             aTACTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCACTAACATACAAGTTGTAtt                                                                              <  1:3189649/91‑2 (MQ=255)
             aTAATTTGGAAGCGCCAGCAATTTTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAtt                                                                              <  2:3187111/91‑2 (MQ=37)
             aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACATACAAGTTGTAtt                                                                              <  1:2082725/91‑2 (MQ=255)
               aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTCGTAttgt                                                                            <  2:743030/91‑4 (MQ=37)
               aaTTTGGAAGCACCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGt                                                                             <  2:1533526/91‑1 (MQ=255)
                 tttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTa                                                                           <  1:2128666/91‑1 (MQ=255)
                  ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                          >  1:2103008/1‑91 (MQ=255)
                  ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                          >  1:3279285/1‑91 (MQ=255)
                  ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                          <  1:1821144/91‑1 (MQ=255)
                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                         <  1:2174195/91‑1 (MQ=255)
                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                         <  1:922941/91‑1 (MQ=255)
                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTCGTATGTTAgt                                                                         >  2:1502656/1‑91 (MQ=255)
                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                         <  1:3434890/91‑2 (MQ=25)
                   tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                         <  1:1300394/91‑2 (MQ=25)
                    ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                        >  2:3527394/1‑91 (MQ=255)
                     gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                       <  2:957436/91‑1 (MQ=255)
                     gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                       <  2:928036/91‑1 (MQ=255)
                     gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                       <  2:928130/91‑1 (MQ=255)
                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                      <  1:16362/91‑1 (MQ=255)
                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                      >  1:2507435/1‑91 (MQ=255)
                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                      >  2:2832670/1‑91 (MQ=255)
                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                      >  2:1247159/1‑91 (MQ=255)
                      aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                      <  1:663743/91‑5 (MQ=21)
                      aaGCGCCAGCAATATTAGTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                      <  2:3142955/91‑1 (MQ=255)
                       aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                     >  1:1501355/1‑91 (MQ=255)
                       aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                     <  2:800678/91‑1 (MQ=255)
                       aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                     <  2:2701721/91‑1 (MQ=255)
                        gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                    >  1:1129613/1‑91 (MQ=255)
                        gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                    >  2:1481598/1‑91 (MQ=255)
                           ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCttc                                                                 <  1:3742098/91‑2 (MQ=17)
                            cAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttt                                                                <  1:360721/91‑1 (MQ=255)
                                aaTATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGa                                                            >  1:3036337/1‑91 (MQ=255)
                                 atatTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGat                                                           >  2:955628/1‑91 (MQ=255)
                                   attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                         <  1:992831/91‑1 (MQ=255)
                                    ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAtt                                                        <  1:2817289/91‑1 (MQ=255)
                                     tattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAttt                                                       <  1:2023544/91‑1 (MQ=255)
                                      attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTa                                                      >  2:2681635/1‑91 (MQ=255)
                                      attaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTa                                                      >  2:3433891/1‑91 (MQ=255)
                                       ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                     <  2:816269/91‑1 (MQ=255)
                                        taATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAt                                                    >  2:444671/1‑91 (MQ=255)
                                        taATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAt                                                    <  1:699566/91‑1 (MQ=255)
                                         aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAtt                                                   >  2:836717/1‑91 (MQ=255)
                                          aTGTTATTCGCTAAATTCAGAGTGGAACCGTTCGGAAGCGCCTCTAACAATACAACTTGTATGTTCGTGGTGCTTTTTTGATATTTAAttt                                                  >  2:938913/1‑91 (MQ=255)
                                             ttATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGc                                               <  2:1323175/91‑1 (MQ=255)
                                             ttATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGc                                               <  2:2893596/91‑1 (MQ=255)
                                              tATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCa                                              <  2:3901961/91‑1 (MQ=255)
                                                ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                            >  1:336849/1‑91 (MQ=255)
                                                 tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                           <  1:998068/91‑1 (MQ=255)
                                                  cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTa                                          >  1:3918082/1‑91 (MQ=255)
                                                  cGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTa                                          >  1:480488/1‑91 (MQ=255)
                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCTCAGGTAc                                         >  2:3500788/1‑91 (MQ=255)
                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                         <  2:3942823/91‑1 (MQ=255)
                                                   gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                         >  1:2375368/1‑91 (MQ=255)
                                                    cTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACt                                        <  1:1556296/91‑1 (MQ=255)
                                                      aaaTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATTCAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTc                                      <  2:1970291/91‑1 (MQ=255)
                                                        aTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCaa                                    <  1:395886/91‑1 (MQ=255)
                                                        aTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCaa                                    <  1:3357194/91‑1 (MQ=255)
                                                           cAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTa                                 >  2:3800257/1‑91 (MQ=255)
                                                            agagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTAt                                >  2:2485537/1‑91 (MQ=255)
                                                             gagTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATg                               >  1:2987780/1‑91 (MQ=255)
                                                              agTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTCTGTTAGTGGTGCTTTTTTTTTTTTTGTTTTCGCAGGTACTTCAATTCTGt                              >  2:3663946/1‑91 (MQ=255)
                                                                tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                            <  2:2386076/91‑1 (MQ=255)
                                                                 ggAACCGTGCGGAAGCGCGTCTAACAATACAACTTGTATGTTAGTGGTGATTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  2:1268119/91‑1 (MQ=255)
                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  1:2658823/91‑1 (MQ=255)
                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  1:2478123/91‑1 (MQ=255)
                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  2:525134/91‑1 (MQ=255)
                                                                 ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                           <  2:1106306/91‑1 (MQ=255)
                                                                    aCCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTAc                        <  2:1601878/91‑1 (MQ=255)
                                                                      cGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGt                      <  1:955628/91‑1 (MQ=255)
                                                                        tGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                    <  2:3874664/91‑1 (MQ=255)
                                                                         gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                   <  2:3120795/91‑1 (MQ=255)
                                                                         gCGGAAGCGCCTCTAACAATACAACTTGTATGTAAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                   <  2:3535247/91‑1 (MQ=255)
                                                                           ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                 <  1:153872/91‑1 (MQ=255)
                                                                           ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                 >  2:855723/1‑91 (MQ=255)
                                                                           ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                 <  1:1133931/91‑1 (MQ=255)
                                                                            gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                <  1:3489423/91‑1 (MQ=255)
                                                                            gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                <  1:3243589/91‑1 (MQ=255)
                                                                               gcgcCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATc             <  2:2056091/91‑1 (MQ=255)
                                                                                  ccactAACATACAAGTTGTATTGTTAGAGGCGCTTTTTTGATATTTAATTTCGCAGGTACTTCAACTATGTATTTACGTAGTTCAATCAtt          >  2:3739286/4‑91 (MQ=21)
                                                                                  ccTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAtt          >  2:2512320/1‑91 (MQ=255)
                                                                                    tctgaCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTga        <  1:3555492/87‑1 (MQ=255)
                                                                                     ctAACAATACAACTTGTATGTTAGTGCTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag       >  2:2241840/1‑91 (MQ=255)
                                                                                      tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg      <  1:3271655/91‑1 (MQ=255)
                                                                                        aCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgaggag    >  1:2326240/1‑91 (MQ=255)
                                                                                          |                                                                                           
CTTAGTTTCCTTGATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGA  >  USA300TCH1516_ALE/576818‑576998

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: