Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 576,920:1 +T 5.9% 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 5.9% 838.0 / 43.4 202intergenic (+343/‑47)gltX/cysEGlutamate‑‑tRNA ligase/Serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base . (73/117);  new base T (4/8);  total (77/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

Only 100 of 217 total aligned reads displayed.

GATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGTT  >  USA300TCH1516_ALE/576830‑577010
                                                                                           |                                                                                          
gATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTa                                                                                             <  2:940041/91‑1 (MQ=255)
 aTAATTTGGAAGCGCCAGCAATTTTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAtt                                                                                          <  2:3187111/91‑2 (MQ=37)
 aTAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCACTAACATACAAGTTGTAtt                                                                                          <  1:261899/91‑2 (MQ=255)
  tAATTTGGAAGCGCCAGCACTATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGGGCGGAAGCGCCTCTAACAACACAACGTGTATg                                                                                          >  2:1538929/1‑91 (MQ=255)
  tAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTAttg                                                                                         <  2:1504544/91‑3 (MQ=37)
   aaTTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTCGTAttgt                                                                                        <  2:743030/91‑4 (MQ=37)
     tttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTa                                                                                       <  1:2128666/91‑1 (MQ=255)
      ttGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAg                                                                                      <  1:471794/91‑1 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                     <  2:2817488/91‑1 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                     <  1:2174195/91‑1 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgt                                                                                     >  1:2084824/1‑91 (MQ=255)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                                     <  1:1300394/91‑2 (MQ=25)
       tGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAga                                                                                     <  1:3434866/91‑2 (MQ=25)
        ggAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtg                                                                                    <  2:269150/91‑1 (MQ=255)
         gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                   <  2:928130/91‑1 (MQ=255)
         gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                   <  2:928036/91‑1 (MQ=255)
         gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                   <  2:956619/91‑1 (MQ=255)
         gAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtgg                                                                                   <  2:3616388/91‑1 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  >  1:2507435/1‑91 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  <  2:3534818/91‑1 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  >  2:1247159/1‑91 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  >  2:3521978/1‑91 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  <  1:16362/91‑1 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCTCTAACAATACAACTTGTATGTTAgtggt                                                                                  <  2:2549735/91‑1 (MQ=255)
          aaGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAgaggc                                                                                  <  1:663777/91‑5 (MQ=21)
           aGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAgtggtg                                                                                 <  2:800678/91‑1 (MQ=255)
            gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                                >  1:1129613/1‑91 (MQ=255)
            gcgcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGc                                                                                <  2:1415258/91‑1 (MQ=255)
              gcCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtt                                                                              <  1:3565018/91‑1 (MQ=255)
               ccAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCACCACTAACATACAAGTTGTATTGTTAGAGGCGCttc                                                                             <  1:3742098/91‑2 (MQ=17)
                  gCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttttt                                                                          >  1:966276/1‑91 (MQ=255)
                  gCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCtttttt                                                                          >  2:1039102/1‑91 (MQ=255)
                   cAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTg                                                                         <  1:2311950/91‑1 (MQ=255)
                       attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                     <  1:3952845/91‑1 (MQ=255)
                       attattAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGatat                                                                     <  1:992934/91‑1 (MQ=255)
                        ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTCTGTTAGTGGTGCTTTTTTTATAtt                                                                    >  2:3671883/1‑91 (MQ=255)
                        ttattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAtt                                                                    <  1:2817289/91‑1 (MQ=255)
                         tattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAttt                                                                   <  1:2023544/91‑1 (MQ=255)
                         tattaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATAttt                                                                   <  1:941512/91‑1 (MQ=255)
                           ttaATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTaa                                                                 <  2:816682/91‑1 (MQ=255)
                            taATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAt                                                                >  2:444671/1‑91 (MQ=255)
                             aaTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTCGTGGTGCTTTTTTGATATTTAAtt                                                               >  2:2837751/1‑91 (MQ=255)
                              aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                              <  1:1155878/91‑1 (MQ=255)
                              aTGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAAttt                                                              <  1:3521128/91‑1 (MQ=255)
                                 ttATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGc                                                           <  2:2893596/91‑1 (MQ=255)
                                  tATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCa                                                          <  2:3901961/91‑1 (MQ=255)
                                   aTTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAg                                                         >  2:3162847/1‑91 (MQ=255)
                                    ttCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAgg                                                        >  1:336849/1‑91 (MQ=255)
                                     tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                                       <  1:836444/91‑1 (MQ=255)
                                     tCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGt                                                       <  1:3659647/91‑1 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                     >  1:2375368/1‑91 (MQ=255)
                                       gCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTAc                                                     <  1:1562302/91‑1 (MQ=255)
                                         tAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACtt                                                   <  1:1735459/91‑1 (MQ=255)
                                            aTTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCaa                                                <  2:2435363/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)
                                                    tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAt                                        <  2:2386076/91‑1 (MQ=255)
                                                    tGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCCTTTTTGATATTTAATCTCGCAGGTACTTCAATTATGTAt                                        >  2:2007604/1‑91 (MQ=255)
                                                     ggAACCGTGCGGAAGCGCGTCTAACAATACAACTTGTATGTTAGTGGTGATTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                       <  2:1268119/91‑1 (MQ=255)
                                                     ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                       <  2:525134/91‑1 (MQ=255)
                                                     ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                       <  1:3353857/91‑1 (MQ=255)
                                                     ggAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTAtt                                       <  2:1106598/91‑1 (MQ=255)
                                                            tGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAg                                <  2:3874664/91‑1 (MQ=255)
                                                             gCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGt                               <  1:3872037/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)
                                                               ggAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTc                             <  1:153872/91‑1 (MQ=255)
                                                                gAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCa                            <  1:3243589/91‑1 (MQ=255)
                                                                      ccactAACATACAAGTTGTATTGTTAGAGGCGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCAtt                      >  2:3739322/4‑91 (MQ=21)
                                                                        tctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTga                    >  2:1481603/1‑91 (MQ=255)
                                                                         ctAACAATACAACTTGTATGTTAGTGGTGCTTTTTTTTTATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag                   >  2:2491433/1‑91 (MQ=255)
                                                                         ctAACAATACAACTTGTATGTTAGTGCTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag                   >  2:2241840/1‑91 (MQ=255)
                                                                         ctAACATACAAGTTGTATTGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgag                   >  1:571765/1‑91 (MQ=38)
                                                                          tAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagg                  <  1:3271655/91‑1 (MQ=255)
                                                                           aaCAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTgagga                 >  2:2056501/1‑91 (MQ=255)
                                                                             cAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGaggagg               >  1:347914/1‑91 (MQ=255)
                                                                               aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGTAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa             <  1:3135194/91‑1 (MQ=255)
                                                                               aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa             <  1:1893712/91‑1 (MQ=255)
                                                                               aTACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaa             <  1:3066306/91‑1 (MQ=255)
                                                                                tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa            <  2:485626/91‑1 (MQ=255)
                                                                                tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa            >  1:287586/1‑91 (MQ=255)
                                                                                tACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGaaa            <  1:1431258/91‑1 (MQ=255)
                                                                                  cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg          <  2:3953200/91‑1 (MQ=255)
                                                                                  cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg          <  2:3000976/91‑1 (MQ=255)
                                                                                  cAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATg          >  1:3845666/1‑91 (MQ=255)
                                                                                   aaCTTGTATGTTAGTGGTGCTTTTTTGTTATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGa         >  2:3839831/1‑91 (MQ=255)
                                                                                   aaCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGa         <  1:1247046/91‑1 (MQ=255)
                                                                                    agttGTATGTTAGTGGTGCTTTTTTTTTTTTTTTTTTCGCAGGTACTTGCATTATGTATTTACGTAGTTCATTCATTGAGGACGATATgtt        >  2:3300004/3‑89 (MQ=255)
                                                                                    aCTTGTATGTTAGTGGTGCTTTTTTGATATTTAATATCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGAt        <  1:2103297/91‑1 (MQ=255)
                                                                                    aCTTGTATGTTAGTGGTGCTTTTTTGATATTGAATTTCGCAGGTACTTCAATTATGTATTTACATAGTTCAATCATTGAGGAGGAAATGAt        <  2:1229805/91‑1 (MQ=255)
                                                                                      ttGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCt      >  2:815359/1‑91 (MQ=255)
                                                                                      ttGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCt      <  2:1763375/91‑1 (MQ=255)
                                                                                       tGTATGTTCGTGGTGCTTTTTTGTTATTTAATTTCCCAGGTTCTTCAATTATGCATCTACGTAGTTCAATCATACAGGAGGAACTGATCtt     >  2:3273206/1‑91 (MQ=255)
                                                                                       tGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCtt     <  2:2470239/91‑1 (MQ=255)
                                                                                        gTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTg    <  1:2741733/91‑1 (MQ=255)
                                                                                        gTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTg    <  2:1697709/91‑1 (MQ=255)
                                                                                         tATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGt   <  1:2476355/91‑1 (MQ=255)
                                                                                         tATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACATAGTTCAATCATTGAGGAGGAAATGATCTTGt   >  2:3935068/1‑91 (MQ=255)
                                                                                           |                                                                                          
GATAATTTGGAAGCGCCAGCAATATTATTAATGTTATTCGCTAAATTCAGAGTGGAACCGTGCGGAAGCGCCTCTAACAATACAACTTGTATGTTAGTGGTGCTTTTTTGATATTTAATTTCGCAGGTACTTCAATTATGTATTTACGTAGTTCAATCATTGAGGAGGAAATGATCTTGTT  >  USA300TCH1516_ALE/576830‑577010

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: