Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,763 G→A 6.0% intergenic (+417/+86) 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,7630GA 6.0% 673.2 / 20.1 201intergenic (+417/+86)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (77/112);  new base A (3/9);  total (80/121)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 213 total aligned reads displayed.

AGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACT  >  USA300TCH1516_ALE/954673‑954853
                                                                                          |                                                                                          
 gCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAgg                                                                                           >  2:3025813/1‑91 (MQ=255)
 gCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAgg                                                                                           >  2:2549782/1‑91 (MQ=255)
 gCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAgg                                                                                           >  1:2887914/1‑91 (MQ=255)
  cgtACAATAGTGTGCAAGTTGGGGGGGGACGACGAAATAAATGTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGa                                                                                          <  1:770751/89‑1 (MQ=255)
     aCAATAGTGTGCAAGTTGGGGTGGGACGACGATATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATaa                                                                                       <  1:3302494/91‑1 (MQ=255)
      cAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAAt                                                                                      <  2:2166929/91‑1 (MQ=255)
       aaTAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAAtt                                                                                     <  2:3553797/91‑1 (MQ=255)
        aTAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAAttt                                                                                    <  1:1086803/91‑1 (MQ=255)
         tAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTc                                                                                   <  1:3432761/91‑1 (MQ=255)
           gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                 >  2:1095808/1‑91 (MQ=255)
           gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                 >  1:2064840/1‑91 (MQ=255)
           gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                 >  1:1920383/1‑91 (MQ=255)
            tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATACTTTCAAt                                                                                >  1:1943694/1‑91 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               <  2:399124/91‑1 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               <  1:3849380/91‑1 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               <  2:3818332/91‑1 (MQ=255)
              tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                              >  2:2883187/1‑91 (MQ=255)
               gCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTg                                                                             <  1:997148/91‑1 (MQ=255)
                cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                            >  1:648553/1‑91 (MQ=255)
                cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                            <  1:562692/91‑1 (MQ=255)
                  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                          <  1:3689637/91‑1 (MQ=255)
                     tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                       >  2:2689603/1‑91 (MQ=255)
                      ggggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAatat                                                                      <  2:1992523/91‑1 (MQ=255)
                        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                    <  1:1042262/91‑1 (MQ=255)
                         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                   <  1:695138/91‑1 (MQ=255)
                         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                   >  1:246843/1‑91 (MQ=255)
                           gggACGACGAACTAAATTTTGCTTAAATATCATTTCTGTCCCACTCCCAGAAAGACAGCGGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                 <  2:2955169/91‑1 (MQ=25)
                           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                 <  2:3359682/91‑1 (MQ=255)
                           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                 >  1:37163/1‑91 (MQ=255)
                           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                 <  1:1601507/91‑1 (MQ=255)
                            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                <  2:1308582/91‑1 (MQ=255)
                            ggACGACAAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                <  2:3638765/91‑1 (MQ=255)
                             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                               <  1:2504114/91‑1 (MQ=255)
                              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                              <  2:3428295/91‑1 (MQ=255)
                               cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                             <  2:830512/91‑1 (MQ=255)
                               cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgct                                                             <  1:3233157/91‑1 (MQ=17)
                                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                            <  1:1116217/91‑1 (MQ=14)
                                 ccgaAATAACTGTTGCGAAAATATCATTTCTGTCCCACTCCCAAAGAGACAGCAGTAGGTTAATTTCAATTTGGAAATATCTTACTGCTGt                                                           <  2:1386371/90‑1 (MQ=37)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                           <  1:3032216/91‑1 (MQ=255)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                           <  2:3645055/91‑1 (MQ=255)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                           <  1:3562936/91‑1 (MQ=255)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                           <  2:612840/91‑1 (MQ=255)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGt                                                           <  1:1647347/91‑1 (MQ=14)
                                  cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                          >  2:2789817/1‑91 (MQ=255)
                                  cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                          <  2:2638096/91‑1 (MQ=255)
                                   gaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttt                                                         >  1:1567309/1‑91 (MQ=255)
                                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                      >  1:2197213/1‑91 (MQ=255)
                                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                      >  1:2541517/1‑91 (MQ=255)
                                       taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                     >  1:2364622/1‑91 (MQ=255)
                                       taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                     <  2:3340608/91‑1 (MQ=255)
                                         aatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtata                                                   <  1:1202988/91‑1 (MQ=255)
                                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGAAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                             <  2:641521/91‑1 (MQ=255)
                                                cGAAAATATCATTTCTGTCCCACTACCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAg                                            <  2:3165808/91‑1 (MQ=255)
                                                 gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                           <  1:1403225/91‑1 (MQ=255)
                                                   aaaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCta                                         <  2:2064953/91‑1 (MQ=255)
                                                    aaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtat                                        >  1:664470/1‑91 (MQ=255)
                                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                       <  2:1263641/91‑1 (MQ=255)
                                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                       <  1:2456668/91‑1 (MQ=255)
                                                      tatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTatat                                      <  2:1602776/91‑1 (MQ=255)
                                                        tCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                    >  1:2899614/1‑91 (MQ=255)
                                                         cATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtttt                                   <  1:1621090/91‑1 (MQ=255)
                                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                  >  2:318555/1‑91 (MQ=255)
                                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                  >  2:2425043/1‑91 (MQ=255)
                                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                  >  2:2374176/1‑91 (MQ=255)
                                                           tttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTg                                 >  2:1476449/1‑91 (MQ=255)
                                                               tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                             <  1:1215232/91‑1 (MQ=255)
                                                                   ccACTCCCAAAAAGACTGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                         <  2:3727908/91‑1 (MQ=255)
                                                                    cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                        >  1:1168468/1‑91 (MQ=255)
                                                                     aCTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtt                       <  2:1966950/91‑1 (MQ=255)
                                                                       tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                     >  1:1602714/1‑91 (MQ=255)
                                                                        cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                    <  1:942276/91‑1 (MQ=255)
                                                                        cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                    >  2:1790864/1‑91 (MQ=255)
                                                                         ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                   >  1:1701717/1‑91 (MQ=255)
                                                                         ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                   >  2:2969917/1‑91 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATACTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  >  2:1546322/1‑91 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  >  1:2220698/1‑91 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  <  2:338824/91‑1 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  <  1:21627/91‑1 (MQ=255)
                                                                           taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                 >  1:34812/2‑91 (MQ=14)
                                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                >  1:3040801/1‑91 (MQ=255)
                                                                             aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt               >  1:2071504/1‑91 (MQ=255)
                                                                              aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt              <  1:2787559/91‑1 (MQ=255)
                                                                                gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa            <  2:1859851/91‑1 (MQ=255)
                                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                    >  2:1802271/1‑81 (MQ=255)
                                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa          >  1:1657947/1‑91 (MQ=255)
                                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa          >  2:776849/1‑91 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         <  2:927728/91‑1 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         <  2:1550731/91‑1 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         >  2:714914/1‑91 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         >  1:387382/1‑91 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         <  1:2571876/91‑1 (MQ=255)
                                                                                    gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt        >  2:2231714/1‑91 (MQ=255)
                                                                                    gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt        >  2:1917310/1‑91 (MQ=255)
                                                                                     cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta       <  1:296111/91‑1 (MQ=255)
                                                                                      agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac      <  1:1790864/91‑1 (MQ=255)
                                                                                      agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac      >  1:256211/1‑91 (MQ=255)
                                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact  >  1:958342/1‑91 (MQ=255)
                                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact  <  2:1177411/91‑1 (MQ=255)
                                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact  >  1:2592991/1‑91 (MQ=255)
                                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact  <  1:1123536/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACT  >  USA300TCH1516_ALE/954673‑954853

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: