Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,763 G→A 7.8% 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 7.8% 789.8 / 37.8 245intergenic (+417/+86)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (88/137);  new base A (8/11);  total (97/148)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

Only 100 of 254 total aligned reads displayed.

AGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACT  >  USA300TCH1516_ALE/954673‑954853
                                                                                          |                                                                                          
aGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAg                                                                                            >  2:866658/1‑91 (MQ=255)
aGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGAAACAAAAAAAAAAAGACAGCAGTAg                                                                                            <  2:1792667/91‑1 (MQ=14)
 gCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAgg                                                                                           >  1:662337/1‑91 (MQ=255)
  cTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGa                                                                                          <  2:2421562/91‑1 (MQ=255)
   ttACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGAt                                                                                         <  1:3137506/91‑1 (MQ=255)
     aCAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATaa                                                                                       <  2:1925403/91‑1 (MQ=255)
     aCAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATaa                                                                                       <  1:20035/91‑1 (MQ=255)
        aTAGTGAGAAAGTTGGGGTGGGACGACGAAATAAATTGTGCGAAAATATCAGTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAAttt                                                                                    <  2:36924/91‑1 (MQ=18)
         tAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTATGATAATTTc                                                                                   >  2:1113/1‑91 (MQ=255)
           gtgtGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                 <  2:3459714/91‑1 (MQ=255)
           gggtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCaa                                                                                 <  2:3459628/89‑1 (MQ=255)
            tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                >  2:2707123/1‑91 (MQ=255)
            tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                <  2:2651071/91‑1 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               >  1:852116/1‑91 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               <  2:3862029/91‑1 (MQ=255)
             gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                               <  1:3093457/91‑1 (MQ=255)
              tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                              >  1:2340697/1‑91 (MQ=255)
                  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                          <  2:718058/91‑1 (MQ=255)
                  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                          >  2:2218529/1‑91 (MQ=255)
                  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                          >  1:3824862/1‑91 (MQ=255)
                  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTACCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                          <  2:1140862/91‑1 (MQ=255)
                     tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGCCAGCAGTAGGATAATTTCAATTTGGACata                                                                       >  2:3562838/1‑91 (MQ=255)
                     tGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAata                                                                       <  1:3887254/91‑1 (MQ=255)
                      ggggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAatat                                                                      <  1:803693/91‑1 (MQ=255)
                       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                     >  2:1345614/1‑91 (MQ=255)
                       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                     >  1:373744/1‑91 (MQ=255)
                        ggTGGGAGGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                    <  2:3509852/91‑2 (MQ=17)
                        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACCGCAGTAGGATAATTTCAATTAGGAAatatgt                                                                    >  1:2216962/1‑89 (MQ=255)
                        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                    <  1:2608751/91‑1 (MQ=255)
                        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                    <  2:302685/91‑1 (MQ=255)
                        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                    <  1:929512/91‑2 (MQ=17)
                         gTGGGACGACGAAATAAATTTTGCGAAGATTTCATTTCTGTCCCACTCAAGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                   <  2:3666978/91‑1 (MQ=21)
                         gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGCCCCACTCCCAAAAAGACAGCAGTAGGCTAATTTCAATTTGGAAATATccc                                                                   >  2:652858/1‑89 (MQ=255)
                          tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                  >  1:129749/1‑91 (MQ=255)
                            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                >  1:2566455/1‑91 (MQ=255)
                            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                <  1:2262156/91‑1 (MQ=255)
                             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTaatg                                                               <  1:1875280/91‑4 (MQ=255)
                             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                               >  2:2268212/1‑91 (MQ=255)
                             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                               <  2:632348/91‑1 (MQ=255)
                              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                              <  2:3212824/91‑1 (MQ=255)
                              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                              <  1:3155631/91‑1 (MQ=255)
                              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                              <  2:3354671/91‑1 (MQ=14)
                                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                            <  1:3374806/91‑1 (MQ=17)
                                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                            <  1:1851459/91‑1 (MQ=14)
                                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGt                                                           <  1:3217311/91‑1 (MQ=17)
                                    aaataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                        <  2:3368161/91‑1 (MQ=255)
                                    aaataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                        >  2:1639915/1‑91 (MQ=255)
                                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAGCGCTGtttttt                                                      >  2:3017930/1‑91 (MQ=255)
                                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                      >  2:2770935/1‑91 (MQ=255)
                                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                      >  1:13601/1‑91 (MQ=255)
                                       taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                     >  2:412666/1‑91 (MQ=255)
                                          atTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAg                                                  <  2:2957272/91‑1 (MQ=255)
                                           ttttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGt                                                 >  1:3725168/1‑91 (MQ=255)
                                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                             <  1:2743018/91‑1 (MQ=255)
                                                cGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAg                                            <  1:2936568/91‑1 (MQ=255)
                                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                       <  2:333117/91‑1 (MQ=255)
                                                       atCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtt                                     >  2:36045/1‑91 (MQ=255)
                                                        tCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                    >  2:1844426/1‑91 (MQ=255)
                                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                  >  2:1672047/1‑91 (MQ=255)
                                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                  >  2:46237/1‑91 (MQ=255)
                                                           tttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTg                                 <  1:430937/91‑1 (MQ=255)
                                                             tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                               <  1:1528672/91‑1 (MQ=255)
                                                             tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                               <  2:1797118/91‑1 (MQ=255)
                                                             tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                               >  2:2852494/1‑91 (MQ=255)
                                                             tCTGCCCCACTCCCAAAAAGACATAGGCAGGGTTATATCTGTTTGTAAATATCTTACTGCTGTGTTTTATAGGTATAGTTATATTTTTTTa                               <  2:2884583/91‑1 (MQ=255)
                                                              cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                              <  1:1638215/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)
                                                               tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                             <  1:2032085/91‑1 (MQ=255)
                                                                 tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                           <  1:1319529/91‑1 (MQ=255)
                                                                 tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                           <  1:1071995/91‑1 (MQ=255)
                                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                          <  2:1540299/91‑1 (MQ=255)
                                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                          <  1:3447669/91‑1 (MQ=255)
                                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                          <  1:3539409/91‑1 (MQ=255)
                                                                  cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                          <  2:764551/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)
                                                                    cACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGt                        <  2:2521274/91‑1 (MQ=255)
                                                                       acacAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                     <  2:3583014/88‑1 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  <  2:3126395/91‑1 (MQ=255)
                                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                  <  2:1437726/91‑1 (MQ=255)
                                                                           taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                 >  2:3386727/2‑91 (MQ=14)
                                                                           taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTCTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                 >  1:2748415/2‑91 (MQ=0)
                                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                >  1:628432/1‑91 (MQ=255)
                                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                <  1:2402992/91‑1 (MQ=255)
                                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                <  1:70232/91‑1 (MQ=255)
                                                                              aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt              >  1:1731241/1‑91 (MQ=255)
                                                                                 aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc           >  1:737700/1‑91 (MQ=255)
                                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa          >  2:556376/1‑91 (MQ=255)
                                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa          >  2:560304/1‑91 (MQ=255)
                                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt         >  1:2371154/1‑91 (MQ=255)
                                                                                     cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta       <  1:1125318/91‑1 (MQ=255)
                                                                                      gggaGGATAATTTCGATGTGGAGATATCTTACTGCTGTTTTTTATAGTAGTAGCTATATTGTTGTAGTTATGTGTTCTAGTTTAACATtac      <  2:98773/88‑1 (MQ=255)
                                                                                      agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac      <  1:3803159/91‑1 (MQ=255)
                                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta    <  2:3117957/91‑1 (MQ=255)
                                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta    >  1:3397434/1‑91 (MQ=255)
                                                                                        tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta    >  1:1588291/1‑91 (MQ=255)
                                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact  <  2:2065995/91‑1 (MQ=255)
                                                                                          |                                                                                          
AGCTTACAATAGTGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACT  >  USA300TCH1516_ALE/954673‑954853

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: