Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,782 A→T 7.9% intergenic (+436/+67) 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,7820AT 7.9% 808.7 / 41.4 241intergenic (+436/+67)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (96/126);  new base T (8/11);  total (104/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 255 total aligned reads displayed.

GTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTT  >  USA300TCH1516_ALE/954692‑954872
                                                                                          |                                                                                          
  tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                                          <  1:3887254/91‑1 (MQ=255)
  tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                                          >  2:3563016/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)
    gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                        >  1:2421566/1‑91 (MQ=255)
    gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTAAAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                        <  2:3450174/91‑1 (MQ=21)
     ggTGGGAGGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                                       <  2:3509852/91‑2 (MQ=17)
     ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                       <  2:302685/91‑1 (MQ=255)
     ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                                       <  1:929512/91‑2 (MQ=17)
      gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                      <  1:1685248/91‑1 (MQ=255)
      gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                      <  1:3447530/91‑1 (MQ=255)
       tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                                     >  1:129965/1‑91 (MQ=255)
        gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                    <  2:2895199/91‑1 (MQ=255)
          gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTaatg                                                                                  <  1:1875280/91‑4 (MQ=255)
          gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                                  <  1:759807/91‑1 (MQ=255)
          gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                                  >  2:2268212/1‑91 (MQ=255)
           acgacgAAATAAATTTTGCGAAAATATCATTTTTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTTAAATTATCCTActgc                                                                                 <  1:3368200/91‑1 (MQ=12)
           acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                                 <  2:3212824/91‑1 (MQ=255)
            cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                                                <  1:728301/91‑1 (MQ=255)
             gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATGATCCTActgctg                                                                               <  1:3029723/91‑1 (MQ=9)
               cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                             >  1:3486307/1‑91 (MQ=255)
               cgaAATAAATGTTGCGAAAATATCATTTCTGTCCCACTCCAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                             <  2:2852727/91‑1 (MQ=255)
                 aaataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                                           >  2:1639915/1‑91 (MQ=255)
                 aaataaatTTTGCGAAAATATCATTTCTGTACAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                                           <  2:3367894/91‑1 (MQ=37)
                  aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                          <  1:2148865/91‑1 (MQ=255)
                   ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                         >  2:2770935/1‑91 (MQ=255)
                   ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                         >  1:13601/1‑91 (MQ=255)
                   ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTAGAAATATCTTAGTGCTGtttttt                                                                         >  2:3018014/1‑91 (MQ=255)
                    taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                                        >  2:412323/1‑91 (MQ=255)
                     aaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtat                                                                       <  1:2231995/91‑1 (MQ=255)
                          ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                                  >  1:3731351/1‑91 (MQ=255)
                            gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                >  2:1922343/1‑91 (MQ=255)
                              gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                              >  2:2503723/1‑91 (MQ=255)
                              gAAAATATCATTTCTGTCACACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                              <  2:2082659/91‑1 (MQ=255)
                                  atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                          <  2:334720/91‑1 (MQ=255)
                                     tCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                                       >  2:1844426/1‑91 (MQ=255)
                                      cATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtttt                                                      >  2:138187/1‑91 (MQ=255)
                                      cATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtttt                                                      <  2:3302073/91‑1 (MQ=255)
                                       aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                     >  2:1672047/1‑91 (MQ=255)
                                       aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                     >  2:46237/1‑91 (MQ=255)
                                       aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                     >  2:656114/1‑91 (MQ=255)
                                        tttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTg                                                    <  1:430937/91‑1 (MQ=255)
                                          tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                                  <  2:1796419/91‑1 (MQ=255)
                                           cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                                 <  2:1717916/91‑1 (MQ=255)
                                           cTGTCCCACTCCCAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                                <  2:3250474/91‑1 (MQ=255)
                                            tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                                >  2:910085/1‑91 (MQ=255)
                                              tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                              <  1:1071995/91‑1 (MQ=255)
                                              tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                              <  1:2672271/91‑1 (MQ=255)
                                               cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                             <  2:1540299/91‑1 (MQ=255)
                                               cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                             <  2:764551/91‑1 (MQ=255)
                                               cccACGTCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                             <  2:3638043/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                                           >  1:2696223/1‑91 (MQ=255)
                                                 cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                                           <  1:1406453/91‑1 (MQ=255)
                                                       cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                     <  2:3126252/91‑1 (MQ=255)
                                                       cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                     >  2:650340/1‑91 (MQ=255)
                                                        taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTCTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                                    >  1:2748415/2‑91 (MQ=0)
                                                         aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                   >  1:692766/1‑91 (MQ=255)
                                                              aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                              <  2:3584281/91‑1 (MQ=255)
                                                               cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                             >  2:1286982/1‑91 (MQ=255)
                                                               cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                             >  2:556376/1‑91 (MQ=255)
                                                               cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                             >  2:560304/1‑91 (MQ=255)
                                                               cagcagTAGGATAATTTCAATTTGGAAATATCTCACTACTGTTTTTCATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                             >  1:3153435/1‑91 (MQ=255)
                                                               cagcagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                             >  2:2414620/1‑91 (MQ=17)
                                                                agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                            >  1:2371154/1‑91 (MQ=255)
                                                                 gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                           <  1:1249491/91‑1 (MQ=255)
                                                                 gcagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                           >  1:19252/1‑91 (MQ=18)
                                                                  cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                          <  1:3783890/91‑1 (MQ=255)
                                                                  cagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                          >  2:1462288/1‑91 (MQ=18)
                                                                   agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                         <  1:3803159/91‑1 (MQ=255)
                                                                    gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                        >  2:966704/3‑91 (MQ=255)
                                                                    gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                        <  1:2605759/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)
                                                                        gATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACtt                    <  1:3057867/91‑1 (MQ=255)
                                                                          tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa                  <  1:1585617/91‑1 (MQ=255)
                                                                             tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa               <  1:1043163/91‑1 (MQ=255)
                                                                             tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa               <  2:3842461/91‑1 (MQ=255)
                                                                             tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa               >  2:1410984/1‑91 (MQ=255)
                                                                              ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc              <  2:2803126/91‑1 (MQ=255)
                                                                               tCAATTTGGAAATATCTTACTGCTGTTTTTTTTTTTTTTTGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct             >  2:2446578/1‑91 (MQ=255)
                                                                                cAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActt            >  1:3221584/1‑91 (MQ=255)
                                                                                 aaTTTGGAAATATCTTACTGCTGTTTTTTGTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc           <  1:3172843/91‑1 (MQ=255)
                                                                                 aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc           >  1:1751990/1‑91 (MQ=255)
                                                                                 aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc           <  2:2969945/91‑1 (MQ=255)
                                                                                  aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct          >  1:913339/1‑91 (MQ=255)
                                                                                  aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct          >  1:2358193/1‑91 (MQ=255)
                                                                                     tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta       <  2:3662900/91‑1 (MQ=255)
                                                                                     tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta       <  2:3181860/91‑1 (MQ=255)
                                                                                      ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATCTTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat      >  2:224206/1‑91 (MQ=255)
                                                                                      ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTTTTTTTTTTTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat      >  2:3376647/1‑91 (MQ=255)
                                                                                      ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat      >  1:2552718/1‑91 (MQ=255)
                                                                                      ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat      >  1:44962/1‑91 (MQ=255)
                                                                                      ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat      >  2:254757/1‑91 (MQ=255)
                                                                                       gAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTTTATTTATTTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtatt     >  2:475000/1‑90 (MQ=255)
                                                                                       gAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtata     <  1:2895137/91‑1 (MQ=255)
                                                                                          |                                                                                          
GTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTT  >  USA300TCH1516_ALE/954692‑954872

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: