Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,677 G→A 11.2% intergenic (+23/‑306) 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,6770GA11.2% 708.2 / 52.8 245intergenic (+23/‑306)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (109/107);  new base A (4/24);  total (113/132)
Fisher's exact test for biased strand distribution p-value = 2.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

AACCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGT  >  USA300TCH1516_ALE/970587‑970767
                                                                                          |                                                                                          
aaCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAg                                                                                            >  1:3843013/1‑91 (MQ=255)
aaCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAg                                                                                            <  1:3601396/91‑1 (MQ=255)
aaCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAg                                                                                            <  1:2038140/91‑1 (MQ=255)
aaCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAg                                                                                            <  1:2806793/91‑1 (MQ=255)
aaCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAg                                                                                            <  2:1337582/91‑1 (MQ=255)
 aCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGt                                                                                           >  2:1206802/1‑91 (MQ=255)
 aCCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGt                                                                                           <  2:1676413/91‑1 (MQ=255)
  ccAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCTGTc                                                                                          >  1:3187459/1‑91 (MQ=255)
  ccAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTc                                                                                          >  1:1330317/1‑91 (MQ=255)
   cAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCt                                                                                         >  1:268733/1‑91 (MQ=255)
      gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                      <  1:2441699/91‑1 (MQ=255)
      gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                      <  1:3518339/91‑1 (MQ=255)
           gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                 >  1:53901/1‑91 (MQ=255)
           gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                 >  2:38921/1‑91 (MQ=255)
            aaGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGa                                                                                >  2:1886434/1‑91 (MQ=255)
            aaGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGa                                                                                >  1:1390388/1‑91 (MQ=255)
              gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                              >  2:1451829/1‑91 (MQ=255)
              gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                              >  1:2998457/1‑91 (MQ=255)
               ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGactg                                                                             <  2:882885/91‑4 (MQ=37)
                 aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                           <  2:2805044/91‑1 (MQ=255)
                  aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                          >  2:2785694/1‑91 (MQ=255)
                  aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                          <  1:2822343/91‑1 (MQ=255)
                   aGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTa                                                                         <  1:864309/91‑1 (MQ=255)
                    gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                        <  1:351578/91‑1 (MQ=37)
                    gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAg                                                                        <  2:387889/91‑1 (MQ=25)
                     tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                       <  1:3726400/91‑1 (MQ=255)
                     tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAGAGACTGGTTAgc                                                                       <  2:2193292/91‑2 (MQ=37)
                     tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgc                                                                       <  1:3709693/91‑2 (MQ=25)
                      gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAGTTATCATTTAAAAATCTAACCAATCTATATAAAGACTGGTTAgct                                                                      <  2:323482/91‑3 (MQ=25)
                       aCTATTGATAATATTGGAACGCTGAGAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     <  1:3688378/91‑1 (MQ=255)
                       aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                     >  2:1098666/1‑91 (MQ=255)
                         tATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttttt                                                                   <  1:2206747/91‑1 (MQ=255)
                         tATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAGctttt                                                                   <  2:1316832/91‑1 (MQ=25)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                  >  2:3237363/1‑91 (MQ=255)
                          aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                  <  1:287658/91‑1 (MQ=21)
                           ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                 >  1:3407630/1‑91 (MQ=255)
                           ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                 <  2:821260/91‑2 (MQ=21)
                            tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                <  1:250578/91‑1 (MQ=255)
                            tGATAATATTAGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttttttaa                                                                <  2:395885/91‑3 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTGAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                              >  1:3946142/1‑91 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                              >  1:1300989/1‑91 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                              >  1:937154/1‑91 (MQ=255)
                              ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                              >  2:3454975/1‑91 (MQ=255)
                                 atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                           <  2:114354/91‑1 (MQ=255)
                                 atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                           <  2:2771394/91‑1 (MQ=255)
                                 atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                          <  2:59580/91‑1 (MQ=255)
                                  tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                          <  1:2145381/91‑1 (MQ=255)
                                   atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtttt                                                         <  1:1564823/91‑1 (MQ=255)
                                   atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGAtttt                                                         >  2:2687988/1‑91 (MQ=255)
                                    ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                        >  1:3351548/1‑91 (MQ=255)
                                    ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                        <  2:1511584/91‑1 (MQ=255)
                                     aggAACGCTGACAACATATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCt                                                       <  2:3341897/90‑1 (MQ=255)
                                      ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                      >  2:2866153/1‑91 (MQ=255)
                                       gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                     >  1:1002549/1‑91 (MQ=255)
                                        aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                    >  1:1979676/1‑91 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   >  2:2069944/1‑91 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   >  1:94610/1‑91 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   <  1:505257/91‑1 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   <  1:2556317/91‑1 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   <  1:1642124/91‑1 (MQ=255)
                                         aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                   <  2:1536748/91‑1 (MQ=255)
                                          cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                  <  1:2107590/91‑1 (MQ=255)
                                           gCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAg                                                 >  2:400696/1‑91 (MQ=255)
                                            cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                >  2:94599/1‑91 (MQ=255)
                                                aaaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGACTGGTTAGTTTTTCCTAGCTTTTCTACAAAGGTGtt                                            >  2:3378102/2‑91 (MQ=255)
                                                 aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                           <  2:1841127/91‑1 (MQ=255)
                                                     tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                       <  2:2290477/91‑1 (MQ=255)
                                                      atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                      <  2:1698476/91‑1 (MQ=255)
                                                      atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                      >  2:834661/1‑91 (MQ=255)
                                                      atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                      >  1:174406/1‑91 (MQ=255)
                                                       tatCGCTAAATAATTATCATTTAAAAAGCTAACCTGTCTTAATATAGATTGGATAGTTTTTTCTTGATTTTCTAAAAAGGTGTTAAAGATa                                     >  2:3678730/1‑91 (MQ=255)
                                                       tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAATGATa                                     <  1:3126275/91‑1 (MQ=255)
                                                          cGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAat                                  <  1:2654916/91‑1 (MQ=255)
                                                            cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                >  1:2339294/1‑91 (MQ=255)
                                                            cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                >  2:3705745/1‑91 (MQ=255)
                                                               aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                             <  1:1273176/91‑1 (MQ=255)
                                                               aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                             <  1:3495047/91‑1 (MQ=255)
                                                                 taatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttatttat                           >  1:1897694/1‑91 (MQ=255)
                                                                   atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                         >  1:1135145/1‑91 (MQ=255)
                                                                    ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                        <  1:1925059/91‑1 (MQ=255)
                                                                      aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                      <  2:1248528/91‑1 (MQ=255)
                                                                       tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                     <  2:295816/91‑1 (MQ=255)
                                                                       tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                     >  1:2006019/1‑91 (MQ=255)
                                                                         aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                   <  2:721578/91‑1 (MQ=255)
                                                                          tttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAcc                  >  1:738960/1‑91 (MQ=255)
                                                                          tttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAcc                  <  2:725566/91‑1 (MQ=255)
                                                                           ttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCa                 >  1:664809/1‑91 (MQ=255)
                                                                             aaaaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtt               >  1:690741/1‑91 (MQ=255)
                                                                               aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt             <  1:3462404/91‑1 (MQ=255)
                                                                               aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt             >  2:485642/1‑91 (MQ=255)
                                                                                  gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          <  1:292693/91‑1 (MQ=255)
                                                                                  acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag          >  1:2618320/2‑91 (MQ=37)
                                                                                   cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga         >  2:1130037/1‑91 (MQ=37)
                                                                                     aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg       <  1:1586018/91‑1 (MQ=255)
                                                                                      acaaGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATCAATTATTTATAATGTTACCCTTTTTAGATGa      >  2:1073120/4‑91 (MQ=255)
                                                                                       ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTATCATTTTGAGAGGaa     >  2:1104516/1‑91 (MQ=255)
                                                                                        ccgTCTTCATATAGATTGGCTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaaa    >  2:2492179/3‑91 (MQ=255)
                                                                                        cAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATTaaa    >  2:673090/1‑91 (MQ=255)
                                                                                          gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt  <  2:917727/91‑1 (MQ=255)
                                                                                          gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt  <  2:3830335/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACCAGGTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGT  >  USA300TCH1516_ALE/970587‑970767

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: