Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,779 G→A 7.8% intergenic (+433/+70) 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,7790GA 7.8% 798.6 / 37.4 245intergenic (+433/+70)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (94/131);  new base A (8/11);  total (102/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 257 total aligned reads displayed.

CAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATA  >  USA300TCH1516_ALE/954689‑954869
                                                                                          |                                                                                          
  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                                          >  2:1039204/1‑91 (MQ=255)
  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                                          <  1:2782166/91‑1 (MQ=255)
   gTTGGGGTGGGACGACGAAATAAATTTTGCAAAAATATCATTTCTGTACCAATCAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGaaa                                                                                         <  2:3057705/91‑1 (MQ=21)
    ttGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAat                                                                                        <  2:3224275/91‑1 (MQ=255)
     tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                                       >  2:3563016/1‑91 (MQ=255)
       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                     >  2:1217216/1‑91 (MQ=255)
       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                     >  2:1345614/1‑91 (MQ=255)
       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                     <  1:1826235/91‑1 (MQ=255)
       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTAAAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                     <  2:3450174/91‑1 (MQ=21)
        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACCGCAGTAGGATAATTTCAATTAGGAAatatgt                                                                                    >  1:2216962/1‑89 (MQ=255)
        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                    <  1:2608751/91‑1 (MQ=255)
        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                    >  2:558478/1‑91 (MQ=255)
         gTGGGACGACGAAATAAATTTTGCGAAGATTTCATTTCTGTCCCACTCAAGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                   <  2:3666978/91‑1 (MQ=21)
         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                   <  1:1685248/91‑1 (MQ=255)
         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                   >  1:344150/1‑91 (MQ=255)
         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                   <  2:2208700/91‑1 (MQ=255)
         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGCCCCACTCCCAAAAAGACAGCAGTAGGCTAATTTCAATTTGGAAATATccc                                                                                   >  2:652858/1‑89 (MQ=255)
          tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                                  >  1:129965/1‑91 (MQ=255)
           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                 <  1:1766161/91‑1 (MQ=255)
           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                 <  2:2895199/91‑1 (MQ=255)
            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                <  1:2262156/91‑1 (MQ=255)
            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                >  2:517704/1‑91 (MQ=255)
             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                               >  1:2532139/1‑91 (MQ=255)
             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                               <  1:759807/91‑1 (MQ=255)
              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                              <  1:2358983/91‑1 (MQ=255)
              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                              <  2:3354672/91‑1 (MQ=14)
               cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                                             <  1:728301/91‑1 (MQ=255)
                gactaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAcggctg                                                                            <  1:3118974/91‑1 (MQ=255)
                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgctg                                                                            <  1:1924023/91‑1 (MQ=255)
                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATGATCCTActgctg                                                                            <  1:3029723/91‑1 (MQ=9)
                 acgaAAGAAATTTTGCGAAAATATCATTTCTGTCCCACTCACAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  2:713705/91‑1 (MQ=255)
                  cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                          >  1:3486307/1‑91 (MQ=255)
                    aaataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                                        <  2:3368161/91‑1 (MQ=255)
                     aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                       >  1:2016575/1‑91 (MQ=255)
                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                      >  1:13601/1‑91 (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:3731351/1‑91 (MQ=255)
                             ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                               <  1:3933243/91‑1 (MQ=255)
                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                             >  2:1614953/1‑91 (MQ=255)
                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                             >  2:1922343/1‑91 (MQ=255)
                                cGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAg                                                            <  1:2936568/91‑1 (MQ=255)
                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                       <  2:333117/91‑1 (MQ=255)
                                       atCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtt                                                     >  2:36045/1‑91 (MQ=255)
                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                  >  2:46237/1‑91 (MQ=255)
                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                  >  2:1672046/1‑91 (MQ=255)
                                             tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                               >  2:2852494/1‑91 (MQ=255)
                                             tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                               <  2:1797118/91‑1 (MQ=255)
                                              cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                              <  2:1717837/91‑1 (MQ=255)
                                              cTGTCCCACTCCCAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                             <  2:3250474/91‑1 (MQ=255)
                                               tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                             >  2:910085/1‑91 (MQ=255)
                                                 tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                           <  1:2703560/91‑1 (MQ=255)
                                                 tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                           <  1:1071995/91‑1 (MQ=255)
                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                          <  2:764551/91‑1 (MQ=255)
                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                          <  1:3539409/91‑1 (MQ=255)
                                                   ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                         <  1:2739267/91‑1 (MQ=255)
                                                    cACTCCCGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                        <  2:3517126/91‑1 (MQ=255)
                                                    cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                        <  1:1406453/91‑1 (MQ=255)
                                                       tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                                     <  2:572931/91‑1 (MQ=255)
                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                  <  2:3126395/91‑1 (MQ=255)
                                                           taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                                 >  2:3386727/2‑91 (MQ=14)
                                                           taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTCTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                                 >  1:2747602/2‑91 (MQ=0)
                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                <  1:2402992/91‑1 (MQ=255)
                                                                gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTCATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                            >  2:2473552/1‑91 (MQ=255)
                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTCACTACTGTTTTTCATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                          >  1:3153435/1‑91 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         <  2:2343135/91‑1 (MQ=255)
                                                                    gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                        >  1:634420/1‑91 (MQ=255)
                                                                    gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                        <  1:1249491/91‑1 (MQ=255)
                                                                     cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                       <  2:383899/91‑1 (MQ=255)
                                                                     cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                       <  1:3783890/91‑1 (MQ=255)
                                                                      agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                      <  1:3803159/91‑1 (MQ=255)
                                                                       gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                     >  2:966704/3‑91 (MQ=255)
                                                                       gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                     <  2:2965215/91‑1 (MQ=255)
                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                    <  2:3117957/91‑1 (MQ=255)
                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                    >  1:1588334/1‑91 (MQ=255)
                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                    >  1:3397173/1‑91 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  >  1:3857819/1‑91 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  <  2:2065995/91‑1 (MQ=255)
                                                                             tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa               <  1:1585617/91‑1 (MQ=255)
                                                                             tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa               >  2:391502/1‑91 (MQ=255)
                                                                              aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt              <  2:2422790/91‑1 (MQ=255)
                                                                              aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt              >  2:1524751/1‑91 (MQ=255)
                                                                                tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa            <  1:1043163/91‑1 (MQ=255)
                                                                                tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa            >  1:810595/1‑91 (MQ=255)
                                                                                tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa            >  2:1410904/1‑91 (MQ=255)
                                                                                 ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc           <  2:2803126/91‑1 (MQ=255)
                                                                                  tCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct          <  1:2678251/91‑1 (MQ=255)
                                                                                   cAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActt         >  1:3221584/1‑91 (MQ=255)
                                                                                    aaTTTGGAAATATCTTACTGCTGTTTTTTGTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc        <  1:3172877/91‑1 (MQ=255)
                                                                                     aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct       >  1:913607/1‑91 (MQ=255)
                                                                                     aTTTGGAAATATCTTACTGCCGTATTTTATAGTAATAGCGATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct       >  2:177936/1‑91 (MQ=255)
                                                                                       ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCttt     <  2:1439194/91‑1 (MQ=255)
                                                                                        tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta    <  2:3662844/91‑1 (MQ=255)
                                                                                        tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta    >  1:256875/1‑91 (MQ=255)
                                                                                        tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta    <  2:3181860/91‑1 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATCTTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:224206/1‑91 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  1:44962/1‑91 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:254757/1‑91 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:2721700/1‑91 (MQ=255)
                                                                                          |                                                                                          
CAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATA  >  USA300TCH1516_ALE/954689‑954869

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: