Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,775 T→A 7.9% intergenic (+429/+74) 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,7750TA 7.9% 794.6 / 43.9 241intergenic (+429/+74)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (94/128);  new base A (8/11);  total (102/139)
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.27e-01

Only 100 of 260 total aligned reads displayed.

TGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTT  >  USA300TCH1516_ALE/954685‑954865
                                                                                          |                                                                                          
tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                            >  2:2731884/1‑91 (MQ=255)
tgtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAt                                                                                            <  2:2650916/91‑1 (MQ=255)
 gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                           <  1:2231623/91‑1 (MQ=255)
 gtgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAtt                                                                                           <  1:156522/91‑1 (MQ=255)
  tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                                          >  1:2340770/1‑91 (MQ=255)
  tgCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAAttt                                                                                          >  1:858629/1‑91 (MQ=255)
   gCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTg                                                                                         <  1:1856282/91‑1 (MQ=255)
    cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                                        >  1:3200782/1‑91 (MQ=255)
       gTTGGGGTGGGACGACGAAATAAATTTTGCAAAAATATCATTTCTGTACCAATCAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGaaa                                                                                     <  2:3057705/91‑1 (MQ=21)
       acccgggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaaa                                                                                     >  2:1055755/5‑91 (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:373744/1‑91 (MQ=255)
           gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                 >  1:3543587/1‑91 (MQ=255)
           gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTAAAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                 <  2:3450174/91‑1 (MQ=21)
            ggTGGGAGGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTGAAATTATcc                                                                                <  2:3509852/91‑2 (MQ=17)
             gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                               <  1:3447530/91‑1 (MQ=255)
             gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGCCCCACTCCCAAAAAGACAGCAGTAGGCTAATTTCAATTTGGAAATATccc                                                                               >  2:652858/1‑89 (MQ=255)
              tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                              >  1:129965/1‑91 (MQ=255)
               gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                             <  2:2895199/91‑1 (MQ=255)
                ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                            >  1:2566455/1‑91 (MQ=255)
                ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                            >  2:517704/1‑91 (MQ=255)
                  acgacgAAATAAATTTTGCGAAAATATCATTTTTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTTAAATTATCCTActgc                                                                          <  1:3368200/91‑1 (MQ=12)
                  acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                          <  2:3212824/91‑1 (MQ=255)
                  acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                          <  1:3155631/91‑1 (MQ=255)
                  acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                          <  2:3354671/91‑1 (MQ=14)
                    gactaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAcggctg                                                                        <  1:3118974/91‑1 (MQ=255)
                    gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATCATCCTActgctg                                                                        <  1:3029696/91‑1 (MQ=9)
                     acgaAAGAAATTTTGCGAAAATATCATTTCTGTCCCACTCACAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                       <  2:713705/91‑1 (MQ=255)
                      cgaAATACATTTTGCGAAAATATCATTTCTGTCCCACTCCAAAAAAGCCAGCAGTAGGATAATTTCTCTTTGGAAACATCTTACTGCTGtt                                                                      >  2:3676616/1‑91 (MQ=25)
                      cgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                      <  1:2270773/91‑1 (MQ=255)
                      cgaAATAAATGTTGCGAAAATATCATTTCTGTCCCACTCCAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                      <  2:2852727/91‑1 (MQ=255)
                        aaatcaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATACTTTCAATTTGGAAATATCTTACTGCTGtgtt                                                                    >  1:3821490/1‑91 (MQ=255)
                        aaataaatTTTGCGAAAATATCATTTCTGTACAAATAAAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttt                                                                    <  2:3367894/91‑1 (MQ=37)
                          ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                  >  1:13601/1‑91 (MQ=255)
                          ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTAGAAATATCTTAGTGCTGtttttt                                                                  >  2:3018014/1‑91 (MQ=255)
                           taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                                 >  2:412666/1‑91 (MQ=255)
                                 ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                           <  1:3933243/91‑1 (MQ=255)
                                     gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                       >  2:2503723/1‑91 (MQ=255)
                                        aaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtat                                                    <  2:392844/91‑1 (MQ=255)
                                             cATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtttt                                               >  2:138187/1‑91 (MQ=255)
                                              aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                              >  2:1672046/1‑91 (MQ=255)
                                              aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                              >  1:825393/1‑91 (MQ=255)
                                                 tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                           <  2:1796435/91‑1 (MQ=255)
                                                  cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                          <  1:1638215/91‑1 (MQ=255)
                                                  cTGTCCCACTCCCAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATGTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                         <  2:3250458/91‑1 (MQ=255)
                                                   tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                         <  2:3307025/91‑1 (MQ=255)
                                                   tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                         <  1:2032085/91‑1 (MQ=255)
                                                     tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                       <  1:1071995/91‑1 (MQ=255)
                                                      cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                      <  2:764551/91‑1 (MQ=255)
                                                      cccACGTCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                      <  2:3638043/91‑1 (MQ=255)
                                                       ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                     <  1:2739267/91‑1 (MQ=255)
                                                       ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                     <  2:1854755/91‑1 (MQ=255)
                                                           tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                                 <  2:572931/91‑1 (MQ=255)
                                                            cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                <  1:144379/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:628432/1‑91 (MQ=255)
                                                                aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                            >  1:692766/1‑91 (MQ=255)
                                                                aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                            <  1:2402992/91‑1 (MQ=255)
                                                                 aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                           >  1:2806627/1‑91 (MQ=255)
                                                                    gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                        <  1:2615231/91‑1 (MQ=255)
                                                                     aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                       <  2:3584281/91‑1 (MQ=255)
                                                                     aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                       >  1:737523/1‑91 (MQ=255)
                                                                      cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                      >  2:1286982/1‑91 (MQ=255)
                                                                      cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                      >  2:556376/1‑91 (MQ=255)
                                                                       agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                     <  2:2343135/91‑1 (MQ=255)
                                                                       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)
                                                                         cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                   <  2:2659039/91‑1 (MQ=255)
                                                                         cagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                   >  2:1462288/1‑91 (MQ=18)
                                                                          agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTGATAGTGATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                  <  2:873108/91‑1 (MQ=255)
                                                                           gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                 >  2:966704/3‑91 (MQ=255)
                                                                           gTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                 <  2:1552644/91‑1 (MQ=255)
                                                                            tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                >  1:1588334/1‑91 (MQ=255)
                                                                            tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                <  2:3118012/91‑1 (MQ=255)
                                                                            tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                >  1:3397173/1‑91 (MQ=255)
                                                                              ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact              >  1:3857819/1‑91 (MQ=255)
                                                                              ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact              <  2:2065995/91‑1 (MQ=255)
                                                                               gATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACtt             <  2:2090399/91‑1 (MQ=255)
                                                                               gATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACtt             <  1:3057867/91‑1 (MQ=255)
                                                                                   aTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATg         >  1:3517148/1‑91 (MQ=255)
                                                                                   aTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATg         >  1:2713195/1‑91 (MQ=255)
                                                                                    tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa        <  1:1043163/91‑1 (MQ=255)
                                                                                    tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa        >  2:1410984/1‑91 (MQ=255)
                                                                                    tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa        <  2:3842461/91‑1 (MQ=255)
                                                                                     ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc       <  2:2803126/91‑1 (MQ=255)
                                                                                     ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc       <  2:2252974/91‑1 (MQ=255)
                                                                                      tCAATTTGGAAATATCTTACTGCTGTTTTTTTTTTTTTTTGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct      >  2:2446578/1‑91 (MQ=255)
                                                                                      tCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct      <  1:2678251/91‑1 (MQ=255)
                                                                                       cAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActt     >  1:3221584/1‑91 (MQ=255)
                                                                                        aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc    >  1:1751990/1‑91 (MQ=255)
                                                                                        aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc    <  2:2969945/91‑1 (MQ=255)
                                                                                         aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct   >  1:913607/1‑91 (MQ=255)
                                                                                         aTTTGGAAATATCTTACTGCCGTATTTTATAGTAATAGCGATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct   >  2:177936/1‑91 (MQ=255)
                                                                                          tttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtctt  <  2:1684951/91‑1 (MQ=255)
                                                                                          |                                                                                          
TGTGCAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTT  >  USA300TCH1516_ALE/954685‑954865

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: