Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,771 T→C 7.8% intergenic (+425/+78) gluD → / ← glpQ_1 NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE954,7710TC 7.8% 770.9 / 39.8 243intergenic (+425/+78)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (93/131);  new base C (8/11);  total (101/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 256 total aligned reads displayed.

ATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACT  >  USA300TCH1516_ALE/954681‑954861
                                                                                          |                                                                                          
 tAGTGTGCAAGTTGGGGTGGGACGACTAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTc                                                                                           >  2:2238639/1‑91 (MQ=255)
  aGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCa                                                                                          >  2:1803999/1‑91 (MQ=255)
   gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                         <  2:3459713/91‑1 (MQ=255)
    tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                        >  2:2707123/1‑91 (MQ=255)
    tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                        <  2:2651226/91‑1 (MQ=255)
     gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                       <  1:3093457/91‑1 (MQ=255)
     gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                       <  2:3505339/91‑1 (MQ=255)
     gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                       <  1:2231623/91‑1 (MQ=255)
      tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                                      >  1:2340770/1‑91 (MQ=255)
      tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                                      >  1:858629/1‑91 (MQ=255)
        cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                                    >  1:3200782/1‑91 (MQ=255)
          aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                                  >  2:1039204/1‑91 (MQ=255)
          aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTACCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                                  <  2:1140862/91‑1 (MQ=255)
           gTTGGGGTGGGACGACGAAATAAATTTTGCAAAAATATCATTTCTGTACCAATCAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGaaa                                                                                 <  2:3057705/91‑1 (MQ=21)
             tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGCCAGCAGTAGGATAATTTCAATTTGGACata                                                                               >  2:3562838/1‑91 (MQ=255)
              ggggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAatat                                                                              <  1:803693/91‑1 (MQ=255)
               gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                             <  1:1826235/91‑1 (MQ=255)
               gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                             >  1:3543587/1‑91 (MQ=255)
               gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTAAAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                             <  2:3450174/91‑1 (MQ=21)
                ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                            >  2:558478/1‑91 (MQ=255)
                 gTGGGACGACGAAATAAATTTTGCGAAGATTTCATTTCTGTCCCACTCAAGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                           <  2:3666978/91‑1 (MQ=21)
                 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                           >  1:344150/1‑91 (MQ=255)
                 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                           <  2:2208700/91‑1 (MQ=255)
                 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGCCCCACTCCCAAAAAGACAGCAGTAGGCTAATTTCAATTTGGAAATATccc                                                                           >  2:652858/1‑89 (MQ=255)
                  tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                          >  1:129749/1‑91 (MQ=255)
                    ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                        >  2:517704/1‑91 (MQ=255)
                     gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTaatg                                                                       <  1:1875280/91‑4 (MQ=255)
                     gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                       >  1:2532139/1‑91 (MQ=255)
                      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                      <  2:3212824/91‑1 (MQ=255)
                      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                      <  2:1400834/91‑1 (MQ=17)
                        gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATGATCCTActgctg                                                                    <  1:3029723/91‑1 (MQ=9)
                         acgaAAGAAATTTTGCGAAAATATCATTTCTGTCCCACTCACAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                   <  2:713705/91‑1 (MQ=255)
                          cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                  >  1:3486307/1‑91 (MQ=255)
                          cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                  <  1:2270773/91‑1 (MQ=255)
                            aaatcaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATACTTTCAATTTGGAAATATCTTACTGCTGtgtt                                                                >  1:3821490/1‑91 (MQ=255)
                            aaataaatTTTGCGAAAATATCATTTCTGTACAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                                <  2:3367894/91‑1 (MQ=37)
                             aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                               <  1:2148865/91‑1 (MQ=255)
                              ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTAGAAATATCTTAGTGCTGtttttt                                                              >  2:3018014/1‑91 (MQ=255)
                               taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                             >  2:412666/1‑91 (MQ=255)
                                aaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtat                                                            <  1:334020/91‑1 (MQ=255)
                                     ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                       <  1:3933243/91‑1 (MQ=255)
                                     ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                       >  1:3731351/1‑91 (MQ=255)
                                       gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                     >  2:1614953/1‑91 (MQ=255)
                                        cGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAg                                                    <  1:2936568/91‑1 (MQ=255)
                                            aaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtat                                                <  2:392844/91‑1 (MQ=255)
                                             atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                               <  2:334720/91‑1 (MQ=255)
                                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                          >  1:825393/1‑91 (MQ=255)
                                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                          >  2:1672046/1‑91 (MQ=255)
                                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                       <  1:314284/91‑1 (MQ=255)
                                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                       >  2:2852494/1‑91 (MQ=255)
                                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                       <  1:1528672/91‑1 (MQ=255)
                                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                       <  2:1796419/91‑1 (MQ=255)
                                                     tCTGCCCCACTCCCAAAAAGACATAGGCAGGGTTATATCTGTTTGTAAATATCTTACTGCTGTGTTTTATAGGTATAGTTATATTTTTTTa                                       <  2:2884583/91‑1 (MQ=255)
                                                      cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                      <  2:1717837/91‑1 (MQ=255)
                                                      cTGTCCCACTCCCAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATGTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                     <  2:3250458/91‑1 (MQ=255)
                                                       tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                     <  1:1321507/91‑1 (MQ=255)
                                                       tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                     <  2:3307025/91‑1 (MQ=255)
                                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                   <  1:1071995/91‑1 (MQ=255)
                                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                   <  1:2672271/91‑1 (MQ=255)
                                                          cccAGTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                  <  2:3637860/91‑1 (MQ=255)
                                                          cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                  <  2:764551/91‑1 (MQ=255)
                                                           ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                 <  2:1854755/91‑1 (MQ=255)
                                                           ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                 >  2:3717954/1‑91 (MQ=255)
                                                           ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                 <  1:2739316/91‑1 (MQ=255)
                                                           ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                 <  1:2739039/91‑1 (MQ=255)
                                                            cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                <  2:2521274/91‑1 (MQ=255)
                                                            cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                >  1:2696223/1‑91 (MQ=255)
                                                               tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                             <  1:1378715/91‑1 (MQ=255)
                                                               tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                             <  2:572931/91‑1 (MQ=255)
                                                               acacAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                             <  2:3583014/88‑1 (MQ=255)
                                                               acaaaaaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                             <  2:2793902/87‑1 (MQ=255)
                                                                  cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                          <  2:1437726/91‑1 (MQ=255)
                                                                  cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                          <  2:3126395/91‑1 (MQ=255)
                                                                   taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                         >  2:3386626/2‑91 (MQ=14)
                                                                   taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTCTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                         >  1:2748415/2‑91 (MQ=0)
                                                                    aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                        >  1:628432/1‑91 (MQ=255)
                                                                    aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                        >  1:692766/1‑91 (MQ=255)
                                                                     aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                       <  1:162884/91‑1 (MQ=255)
                                                                      aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                      >  1:1731241/1‑91 (MQ=255)
                                                                        gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTCATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                    >  2:2473552/1‑91 (MQ=255)
                                                                         aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                   >  1:737523/1‑91 (MQ=255)
                                                                          cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                  >  2:1286982/1‑91 (MQ=255)
                                                                           agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                 <  2:2343135/91‑1 (MQ=255)
                                                                           agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                 >  1:2371250/1‑91 (MQ=255)
                                                                            gcagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                >  1:19252/1‑91 (MQ=18)
                                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta               <  1:1132892/91‑1 (MQ=255)
                                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta               <  1:3783890/91‑1 (MQ=255)
                                                                              agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTGATAGTGATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac              <  2:873108/91‑1 (MQ=255)
                                                                               gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact             >  2:966704/3‑91 (MQ=255)
                                                                               gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact             <  2:2965215/91‑1 (MQ=255)
                                                                               gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact             <  1:2605759/91‑1 (MQ=255)
                                                                                tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta            >  1:3397173/1‑91 (MQ=255)
                                                                                tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta            <  2:3117957/91‑1 (MQ=255)
                                                                                  ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact          >  1:3857819/1‑91 (MQ=255)
                                                                                  ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact          <  2:2065995/91‑1 (MQ=255)
                                                                                     tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa       <  1:1585617/91‑1 (MQ=255)
                                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt      >  2:1524751/1‑91 (MQ=255)
                                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt      <  2:2422790/91‑1 (MQ=255)
                                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa    <  2:3842461/91‑1 (MQ=255)
                                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa    <  1:1043163/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACT  >  USA300TCH1516_ALE/954681‑954861

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: