Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,681 T→A 9.1% intergenic (+27/‑302) USA300HOU_RS04805 → / → USA300HOU_RS04810 Ureidoglycolate lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE970,6810TA 9.1% 471.8 / 31.0 143intergenic (+27/‑302)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (77/53);  new base A (5/8);  total (82/61)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 150 total aligned reads displayed.

GGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAA  >  USA300TCH1516_ALE/970592‑970771
                                                                                         |                                                                                          
ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                           >  1:169672/1‑91 (MQ=255)
ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                           >  1:1150621/1‑91 (MQ=255)
ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                           >  1:895066/1‑91 (MQ=255)
ggTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCttt                                                                                           >  1:3661961/1‑91 (MQ=255)
  tgatgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtat                                                                                         >  1:1329635/1‑91 (MQ=255)
   gatgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtata                                                                                        <  1:3098727/91‑1 (MQ=255)
   gatgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCtatata                                                                                        <  1:3300591/91‑1 (MQ=255)
   gatgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCtatata                                                                                        <  1:3841080/91‑1 (MQ=255)
    atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       <  1:3839928/91‑1 (MQ=255)
    atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       >  1:3152254/1‑91 (MQ=255)
    atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       >  1:2881539/1‑91 (MQ=255)
    atgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                       >  1:3633508/1‑91 (MQ=255)
     tgaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatata                                                                                      >  1:3437196/1‑91 (MQ=255)
      gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                     >  1:635472/1‑91 (MQ=255)
        aGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAt                                                                                   <  1:2281834/91‑1 (MQ=255)
        aGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAt                                                                                   <  1:1533622/91‑1 (MQ=255)
         gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                                  <  1:1530254/91‑1 (MQ=255)
         gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                                  >  1:3949865/1‑91 (MQ=255)
         gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                                  <  1:3394597/91‑1 (MQ=255)
           tAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATCATCATTTAAAAAGCTAACCAGTCTTTATATAGATTgg                                                                                <  1:3518595/91‑1 (MQ=255)
            aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                               <  1:2567558/91‑1 (MQ=255)
            aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                               <  1:2567476/91‑1 (MQ=255)
             aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                              <  1:1752346/91‑1 (MQ=255)
              aGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTa                                                                             >  1:740163/1‑91 (MQ=255)
              aGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTa                                                                             <  1:1071195/91‑1 (MQ=25)
               gtggcTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            <  1:3729002/87‑1 (MQ=255)
               gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGttta                                                                            >  1:2819413/1‑89 (MQ=255)
               gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            >  1:1785183/1‑91 (MQ=255)
               gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            >  1:2321748/1‑91 (MQ=255)
               gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                            >  1:2241683/1‑91 (MQ=255)
               gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAg                                                                            <  1:2545295/91‑1 (MQ=25)
                tNACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                           <  1:2424703/91‑1 (MQ=255)
                 gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtt                                                                          >  1:3587310/1‑91 (MQ=255)
                 gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtt                                                                          <  1:2236977/91‑1 (MQ=255)
                 gACTATAGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgct                                                                          <  1:1486021/91‑3 (MQ=21)
                  aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                         >  1:2721327/1‑91 (MQ=255)
                   cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                        >  1:3497940/1‑91 (MQ=255)
                   cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                        >  1:727457/1‑91 (MQ=255)
                   cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                        <  1:1115202/91‑1 (MQ=255)
                      ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                     >  1:2503455/1‑91 (MQ=255)
                      ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                     >  1:687544/1‑91 (MQ=255)
                      ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                     <  1:3634668/91‑2 (MQ=18)
                       tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                    <  1:2162058/91‑1 (MQ=255)
                        gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                   >  1:4050868/1‑91 (MQ=255)
                        gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                   >  1:1627906/1‑91 (MQ=255)
                        gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                   >  1:1495703/1‑91 (MQ=255)
                        gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGctttttaaa                                                                   <  1:2790408/91‑4 (MQ=18)
                         ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                  >  1:2000999/1‑91 (MQ=255)
                            atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                               <  1:1989374/91‑1 (MQ=255)
                             tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                              <  1:1853411/91‑1 (MQ=255)
                             tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                              >  1:936428/1‑91 (MQ=255)
                             tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                              <  1:1854385/91‑1 (MQ=255)
                              atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtttt                                                             >  1:220320/1‑91 (MQ=255)
                               ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                            >  1:1222190/1‑91 (MQ=255)
                               ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                            >  1:1221386/1‑91 (MQ=255)
                                tGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                           <  1:1745798/91‑1 (MQ=255)
                                 ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                          <  1:505193/91‑1 (MQ=255)
                                   aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                        >  1:1899867/1‑91 (MQ=255)
                                      gCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAg                                                     >  1:2828641/1‑91 (MQ=255)
                                         gACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtg                                                  <  1:2846493/91‑1 (MQ=255)
                                         gACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtg                                                  >  1:134135/1‑91 (MQ=255)
                                          aCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtgt                                                 >  1:3998968/1‑91 (MQ=255)
                                           cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                >  1:1081365/1‑91 (MQ=255)
                                           cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                >  1:1680310/1‑91 (MQ=255)
                                           cAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGtt                                                >  1:1081734/1‑91 (MQ=255)
                                            aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                               <  1:3605926/91‑1 (MQ=255)
                                              cTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTaaa                                             <  1:8235/91‑1 (MQ=255)
                                               ttATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAg                                            <  1:1547668/91‑1 (MQ=255)
                                                tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                           <  1:3148587/91‑1 (MQ=255)
                                                 atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                          >  1:3322540/1‑91 (MQ=255)
                                                      gCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatt                                     >  1:2468974/1‑91 (MQ=255)
                                                      gCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatt                                     <  1:2724400/91‑1 (MQ=255)
                                                      gCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatt                                     >  1:2885710/1‑91 (MQ=255)
                                                       cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                    >  1:3685638/1‑91 (MQ=255)
                                                          aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                                 >  1:1733496/1‑91 (MQ=255)
                                                                  tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                         <  1:2916717/91‑1 (MQ=255)
                                                                  tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                         >  1:3979061/1‑91 (MQ=255)
                                                                   cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                        <  1:3088848/91‑1 (MQ=255)
                                                                    aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                       <  1:2136222/91‑1 (MQ=255)
                                                                      ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                     >  1:2947073/1‑91 (MQ=255)
                                                                      ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                     <  1:2779505/91‑1 (MQ=255)
                                                                         aaaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAttt                  <  1:1917819/91‑1 (MQ=255)
                                                                          aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 <  1:1895792/91‑1 (MQ=255)
                                                                          aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 <  1:597612/91‑1 (MQ=255)
                                                                          aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 >  1:3184788/1‑91 (MQ=255)
                                                                          aaaGCTAACCAGACTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                 >  1:3183499/1‑91 (MQ=255)
                                                                            aGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTga               >  1:2487607/1‑91 (MQ=255)
                                                                              cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga             >  1:2000797/1‑91 (MQ=255)
                                                                                aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg           <  1:32952/91‑1 (MQ=255)
                                                                                aaCCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg           >  1:1075391/1‑91 (MQ=38)
                                                                                aaCCAATCTATATAAAGACTGGTTAGCTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg           >  1:1158278/1‑91 (MQ=25)
                                                                                 aCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGa          >  1:247996/1‑91 (MQ=255)
                                                                                  ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaa         >  1:3669772/1‑91 (MQ=255)
                                                                                  ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaa         >  1:2137997/1‑91 (MQ=255)
                                                                                  ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaa         >  1:1770155/1‑91 (MQ=255)
                                                                                    aGTCTTTATATAGATTGGTTAGTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt      <  1:2950836/91‑1 (MQ=255)
                                                                                     gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt      <  1:964802/91‑1 (MQ=255)
                                                                                      tCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg     >  1:1120809/1‑91 (MQ=255)
                                                                                       cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa    <  1:1667280/91‑1 (MQ=255)
                                                                                         ttATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGaaa  >  1:333364/1‑91 (MQ=255)
                                                                                         |                                                                                          
GGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAA  >  USA300TCH1516_ALE/970592‑970771

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: