Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 954,787 T→C 7.2% intergenic (+441/+62) 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,7870TC 7.2% 772.6 / 33.9 237intergenic (+441/+62)gluD/glpQ_1NAD‑specific glutamate dehydrogenase/Glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (95/125);  new base C (8/9);  total (103/134)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 248 total aligned reads displayed.

GGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAG  >  USA300TCH1516_ALE/954697‑954877
                                                                                          |                                                                                          
ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                            <  1:2608751/91‑1 (MQ=255)
ggTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCt                                                                                            <  2:302685/91‑1 (MQ=255)
 gTGGGACGACGAAATAAATTTTGCGAAGATTTCATTTCTGTCCCACTCAAGAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                           <  2:3666978/91‑1 (MQ=21)
 gTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCtt                                                                                           <  1:3447530/91‑1 (MQ=255)
  tGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTa                                                                                          >  1:129749/1‑91 (MQ=255)
   gggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAc                                                                                         <  2:2895199/91‑1 (MQ=255)
    ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                        >  2:517704/1‑91 (MQ=255)
    ggACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAct                                                                                        <  1:2262156/91‑1 (MQ=255)
     gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                                       >  2:2268212/1‑91 (MQ=255)
     gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                                       >  1:2532139/1‑91 (MQ=255)
     gacgacGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActg                                                                                       <  2:632348/91‑1 (MQ=255)
      acgacgAAATAAATTTTGCGAAAATATCATTTTTGTCCCACTCCCAAAAAGACAGCAGTAAGATATTTCCAAATTTAAATTATCCTActgc                                                                                      <  1:3368200/91‑1 (MQ=12)
      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgc                                                                                      <  2:3212824/91‑1 (MQ=255)
      acgacgAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATTATCCTActgc                                                                                      <  2:3354672/91‑1 (MQ=14)
        gactaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAcggctg                                                                                    <  1:3118974/91‑1 (MQ=255)
        gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTActgctg                                                                                    <  2:922124/91‑1 (MQ=255)
        gacgaAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAAACAGCAGTAAGATATTTCCAAATTGAAATGATCCTActgctg                                                                                    <  1:3029723/91‑1 (MQ=9)
         acgaAAGAAATTTTGCGAAAATATCATTTCTGTCCCACTCACAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGt                                                                                   <  2:713705/91‑1 (MQ=255)
          cgaAATAAATGTTGCGAAAATATCATTTCTGTCCCACTCCAAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtt                                                                                  <  2:2852727/91‑1 (MQ=255)
              ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTAGCGCTGtttttt                                                                              >  2:3017930/1‑91 (MQ=255)
              ataaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGtttttt                                                                              >  2:2770935/1‑91 (MQ=255)
               taaatTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTta                                                                             >  2:412666/1‑91 (MQ=255)
                   ttttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGt                                                                         >  1:3725168/1‑91 (MQ=255)
                     ttGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTaa                                                                       <  1:3933243/91‑1 (MQ=255)
                       gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                     >  2:1614953/1‑91 (MQ=255)
                       gCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATa                                                                     <  1:2743018/91‑1 (MQ=255)
                        cGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAg                                                                    <  1:2936568/91‑1 (MQ=255)
                         gAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                                   >  2:2503723/1‑91 (MQ=255)
                         gAAAATATCATTTCTGTCACACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGc                                                                   <  2:2082659/91‑1 (MQ=255)
                             atatCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCtata                                                               <  2:334720/91‑1 (MQ=255)
                               atCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtt                                                             >  2:36045/1‑91 (MQ=255)
                                 cATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAtttt                                                           >  2:138187/1‑91 (MQ=255)
                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                          >  2:1672046/1‑91 (MQ=255)
                                  aTTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATAttttt                                                          >  2:46237/1‑91 (MQ=255)
                                   tttCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTg                                                         <  1:430937/91‑1 (MQ=255)
                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                                       <  2:1797118/91‑1 (MQ=255)
                                     tCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTa                                                       <  1:1528672/91‑1 (MQ=255)
                                      cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                                      <  2:1717837/91‑1 (MQ=255)
                                      cTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAt                                                      <  1:1638215/91‑1 (MQ=255)
                                       tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                                     <  1:2032085/91‑1 (MQ=255)
                                       tGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTAtt                                                     <  2:3307025/91‑1 (MQ=255)
                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                                   <  1:1319529/91‑1 (MQ=255)
                                         tCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTa                                                   <  1:2703560/91‑1 (MQ=255)
                                          cccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTAt                                                  <  2:764551/91‑1 (MQ=255)
                                           ccACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATg                                                 <  1:2739267/91‑1 (MQ=255)
                                               tCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                                             <  1:1378715/91‑1 (MQ=255)
                                               acacAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGtttt                                             <  2:3583014/88‑1 (MQ=255)
                                                  cAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTa                                          <  2:3126395/91‑1 (MQ=255)
                                                   taaaaaaCAGCAGTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAg                                         >  2:3386727/2‑91 (MQ=14)
                                                    aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                        >  1:692766/1‑91 (MQ=255)
                                                    aaaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGt                                        >  1:628432/1‑91 (MQ=255)
                                                     aaaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGtt                                       >  1:2806627/1‑91 (MQ=255)
                                                      aaGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGttt                                      >  1:1731241/1‑91 (MQ=255)
                                                         aCAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAAc                                   >  1:737523/1‑91 (MQ=255)
                                                          cagcagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACa                                  >  2:1286982/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:2371250/1‑91 (MQ=255)
                                                            gcagTAAGATATTTCCAAATTGAAATTATCCTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACAtt                                >  1:19252/1‑91 (MQ=18)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               <  2:383899/91‑1 (MQ=255)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               <  1:3783890/91‑1 (MQ=255)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               >  2:2032595/1‑91 (MQ=255)
                                                             cagTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATta                               <  1:1125318/91‑1 (MQ=255)
                                                              agTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtac                              <  1:3803159/91‑1 (MQ=255)
                                                               gcaGGATAATTTCAATTTGGAAATATCTTACTGCTGTCTTTTCTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtact                             >  2:966704/3‑91 (MQ=255)
                                                                tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                            <  2:3117957/91‑1 (MQ=255)
                                                                tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                            >  1:3397434/1‑91 (MQ=255)
                                                                tAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATtacta                            >  1:1588291/1‑91 (MQ=255)
                                                                     tAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTaa                       >  2:391502/1‑91 (MQ=255)
                                                                      aaTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAAt                      >  2:1524751/1‑91 (MQ=255)
                                                                       aTTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATg                     >  1:3517148/1‑91 (MQ=255)
                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa                    >  2:1410984/1‑91 (MQ=255)
                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa                    <  1:1043163/91‑1 (MQ=255)
                                                                        tttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGa                    >  1:810595/1‑91 (MQ=255)
                                                                         ttCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAc                   <  2:2803126/91‑1 (MQ=255)
                                                                          tCAATTTGGAAATATCTTACTGCTGTTTTTTTTTTTTTTTGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct                  >  2:2446578/1‑91 (MQ=255)
                                                                          tCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGAct                  <  1:2678251/91‑1 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTGTAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                <  1:3172877/91‑1 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                <  1:3246080/91‑1 (MQ=255)
                                                                            aaTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtc                <  2:2969945/91‑1 (MQ=255)
                                                                             aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct               >  1:1782353/1‑91 (MQ=255)
                                                                             aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct               >  1:2358193/1‑91 (MQ=255)
                                                                             aTTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct               >  1:913339/1‑91 (MQ=255)
                                                                             aTTTGGAAATATCTTACTGCCGTATTTTATAGTAATAGCGATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGActtct               >  2:177936/1‑91 (MQ=255)
                                                                               ttGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCttt             <  2:1439194/91‑1 (MQ=255)
                                                                                tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta            <  2:3662844/91‑1 (MQ=255)
                                                                                tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta            >  1:256875/1‑91 (MQ=255)
                                                                                tGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTta            <  2:3181860/91‑1 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATCTTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           >  2:224206/1‑91 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTTTTTTTTTTTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           >  2:3376647/1‑91 (MQ=255)
                                                                                 ggAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTtat           >  2:254757/1‑91 (MQ=255)
                                                                                    aaTATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATAtt        >  1:18580/1‑91 (MQ=255)
                                                                                     atatCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATAttt       <  2:326445/91‑1 (MQ=255)
                                                                                       atCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTAt     >  1:3207176/1‑91 (MQ=255)
                                                                                        tCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATc    <  2:53905/91‑1 (MQ=255)
                                                                                         cTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCa   <  2:184744/91‑1 (MQ=255)
                                                                                         cTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCa   <  1:1071445/91‑1 (MQ=255)
                                                                                          ttACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTTGTGTGACATTACTACTTAATGACTTCTTTATATTTATCAg  <  1:1414439/91‑1 (MQ=255)
                                                                                          ttACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAg  >  1:738174/1‑91 (MQ=255)
                                                                                          |                                                                                          
GGTGGGACGACGAAATAAATTTTGCGAAAATATCATTTCTGTCCCACTCCCAAAAAGACAGCAGTAGGATAATTTCAATTTGGAAATATCTTACTGCTGTTTTTTATAGTAATAGCTATATTTTTGTATTTATGTTTTCTAGTTTAACATTACTACTTAATGACTTCTTTATATTTATCAG  >  USA300TCH1516_ALE/954697‑954877

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