Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,779 G→A 5.4% 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 5.4% 751.0 / 18.6 222intergenic (+433/+70)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (87/123);  new base A (3/9);  total (90/132)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 237 total aligned reads displayed.

CAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATA  >  USA300TCH1516_ALE/954689‑954869
                                                                                          |                                                                                          
cAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTgg                                                                                            <  1:562692/91‑1 (MQ=255)
  aGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGaa                                                                                          <  1:3689655/91‑1 (MQ=255)
       gggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATc                                                                                     >  1:2208651/1‑91 (MQ=255)
        ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                    <  1:1042262/91‑1 (MQ=255)
           gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                 <  1:1601507/91‑1 (MQ=255)
            ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                <  1:1470596/91‑1 (MQ=255)
            ggACGACAAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                <  2:3638765/91‑1 (MQ=255)
             gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                               <  1:2504114/91‑1 (MQ=255)
              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                              <  2:2430994/91‑1 (MQ=255)
              acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                              <  2:383223/91‑1 (MQ=255)
               cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                                             <  2:830512/91‑1 (MQ=255)
               cgacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgct                                                                             >  1:649387/1‑91 (MQ=255)
                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                                            <  1:1173277/91‑1 (MQ=14)
                gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgctg                                                                            <  1:1116217/91‑1 (MQ=14)
                 ccgaAATAACTGTTGCGAAAATATCATTTCTGTCCCACTCCCAAAGAGACAGCAGTAGGTTAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  2:1386371/90‑1 (MQ=37)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           >  2:1319537/1‑91 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  1:3562936/91‑1 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  1:3032216/91‑1 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  2:3645055/91‑1 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  1:171196/91‑1 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                           <  1:1239425/91‑1 (MQ=255)
                 acgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGt                                                                           <  1:1647347/91‑1 (MQ=14)
                     aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                       <  2:3356542/91‑1 (MQ=255)
                     aataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGttttt                                                                       <  1:1616603/91‑1 (MQ=255)
                      ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                      >  1:2541517/1‑91 (MQ=255)
                       taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                                     >  2:3248829/1‑91 (MQ=255)
                        aaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtat                                                                    <  2:2208651/91‑1 (MQ=255)
                         aatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTtata                                                                   <  1:1202988/91‑1 (MQ=255)
                              tGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAAt                                                              <  1:796148/91‑1 (MQ=255)
                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                             >  1:666141/1‑91 (MQ=255)
                               gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGAAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                             <  2:641521/91‑1 (MQ=255)
                                 gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                           <  1:1403225/91‑1 (MQ=255)
                                 gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                           <  2:317614/91‑1 (MQ=255)
                                 gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                           >  1:3394345/1‑91 (MQ=255)
                                   aaaTATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCta                                                         <  2:2064840/91‑1 (MQ=255)
                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                       <  1:2457017/91‑1 (MQ=255)
                                     atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                       <  2:1263641/91‑1 (MQ=255)
                                        tCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                                    >  1:2899614/1‑91 (MQ=255)
                                        tCATTTCTGTCACACTCCCAAAAAGACAGCAGTAGGATAATGTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttt                                                    <  2:1569757/91‑1 (MQ=255)
                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                  >  2:3512827/1‑91 (MQ=255)
                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                  >  2:424085/1‑91 (MQ=255)
                                          aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                  >  2:2374176/1‑91 (MQ=255)
                                            ttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGt                                                >  1:218494/1‑91 (MQ=255)
                                               tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                             <  2:1007042/91‑1 (MQ=255)
                                                 tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                           <  1:3018947/91‑1 (MQ=255)
                                                 tCCCACTACAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                           <  2:2698980/91‑1 (MQ=255)
                                                     aCTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtt                                       <  1:1895771/91‑1 (MQ=255)
                                                        cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                    >  2:1790864/1‑91 (MQ=255)
                                                        cccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                    <  1:942276/91‑1 (MQ=255)
                                                         ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                                   <  2:1980982/91‑1 (MQ=255)
                                                         ccAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCt                                   >  2:2969917/1‑91 (MQ=255)
                                                          cAAAAAGACAGCAGTAGGATACTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                  >  2:1546322/1‑91 (MQ=255)
                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                  <  2:1934941/91‑1 (MQ=255)
                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                  <  1:21627/91‑1 (MQ=255)
                                                          cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                  >  2:3101893/1‑91 (MQ=255)
                                                            aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                >  1:3040801/1‑91 (MQ=255)
                                                             aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                               >  1:2071504/1‑91 (MQ=255)
                                                              aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                              <  1:2787559/91‑1 (MQ=255)
                                                              aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                              <  2:2249663/91‑1 (MQ=255)
                                                              aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                              <  2:2830173/91‑1 (MQ=255)
                                                               aGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTa                             <  1:239192/91‑1 (MQ=255)
                                                                gACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTaa                            <  2:1859851/91‑1 (MQ=255)
                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTc                                    <  1:1802271/81‑1 (MQ=255)
                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                          >  1:1657947/1‑91 (MQ=255)
                                                                  cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                          >  2:776849/1‑91 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         >  2:714914/1‑91 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         >  2:2670256/1‑91 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         <  2:927728/91‑1 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         >  1:387382/1‑91 (MQ=255)
                                                                   agcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAt                         <  2:1550731/91‑1 (MQ=255)
                                                                    gcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                        >  1:3248190/1‑91 (MQ=255)
                                                                     cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                       >  1:634501/1‑91 (MQ=255)
                                                                      agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                      >  1:256211/1‑91 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  >  1:2592991/1‑91 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  <  1:1123536/91‑1 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  >  1:958342/1‑91 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  <  2:1177411/91‑1 (MQ=255)
                                                                          ggATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTactact                  >  2:2942857/1‑91 (MQ=255)
                                                                           gATAATTTCAATTTGGAAATATGTTACTGCTGTTTTTTATAGTAATAGCGATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACtt                 <  2:2952050/91‑1 (MQ=255)
                                                                            aTAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTa                >  1:1584935/1‑91 (MQ=255)
                                                                            aTAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTa                <  2:3759106/91‑1 (MQ=255)
                                                                              aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt              >  1:272457/1‑91 (MQ=255)
                                                                              aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt              <  2:256231/91‑1 (MQ=255)
                                                                                 ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc           >  1:2898843/1‑91 (MQ=255)
                                                                                 ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc           <  2:552398/91‑1 (MQ=255)
                                                                                   cAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActt         <  1:2067133/91‑1 (MQ=255)
                                                                                     aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct       <  1:2451481/91‑1 (MQ=255)
                                                                                     aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct       >  1:3503792/1‑91 (MQ=255)
                                                                                     aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct       >  2:572815/1‑91 (MQ=255)
                                                                                      tttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtctt      <  1:3194132/91‑1 (MQ=255)
                                                                                       ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTTGTTTAACATTACTACTTAATGACTTCttt     <  2:3381659/91‑1 (MQ=255)
                                                                                       ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCttt     >  1:3765888/1‑91 (MQ=255)
                                                                                       ttGGAAATATCTTACTGCTGTTTCTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTTGTTTAACATTACTACTTAATGACTTCttt     <  2:3381797/91‑1 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:2106122/1‑91 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   <  1:3835510/91‑1 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:2742868/1‑91 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   <  2:2028801/91‑1 (MQ=255)
                                                                                         ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat   >  2:3394144/1‑91 (MQ=255)
                                                                                          gAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtata  >  1:1177380/1‑91 (MQ=255)
                                                                                          gAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtata  >  2:1822699/1‑91 (MQ=255)
                                                                                          |                                                                                          
CAAGTTGGGGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATA  >  USA300TCH1516_ALE/954689‑954869

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: