Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,690 T→C 10.7% intergenic (+36/‑293) 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,6900TC10.7% 706.4 / 58.1 251intergenic (+36/‑293)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (111/113);  new base C (6/21);  total (117/134)
Fisher's exact test for biased strand distribution p-value = 7.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

AGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGATAT  >  USA300TCH1516_ALE/970600‑970780
                                                                                          |                                                                                          
 gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                                           >  2:1451829/1‑91 (MQ=255)
    aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                                        <  2:2805044/91‑1 (MQ=255)
    aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGt                                                                                        <  2:2422600/91‑1 (MQ=38)
    aaaGGGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGt                                                                                        <  1:506560/91‑1 (MQ=38)
     aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                                       >  2:2785694/1‑91 (MQ=255)
       gTGACTATTGATAATATTGGAACGCTGGCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                                     <  2:1654213/91‑1 (MQ=37)
       gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                                     >  1:335328/1‑91 (MQ=255)
       gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                                     >  1:12490/1‑91 (MQ=255)
       gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAg                                                                                     <  2:387889/91‑1 (MQ=25)
       gTGACTATTGATAATATTGGAACGCGGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                                     <  1:3666667/91‑1 (MQ=37)
        tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAGAGACTGGTTAgc                                                                                    <  2:2193292/91‑2 (MQ=37)
        tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgc                                                                                    <  1:3709680/91‑2 (MQ=25)
        cgACTATTAAAAATTTTGAAACGCTGACACCTTATGTCGCTAAATAGTTATCATTTAAAAAGCTAACCGATCTATATAGAGACTGGTTAgc                                                                                    <  2:1807975/90‑2 (MQ=255)
         gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAGTTATCATTTAAAAATCTAACCAATCTATATAAAGACTGGTTAgct                                                                                   <  2:323482/91‑3 (MQ=25)
          aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                                  >  2:1098666/1‑91 (MQ=255)
          aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                                  >  1:1278706/1‑91 (MQ=255)
            tATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAGctttt                                                                                <  2:1316832/91‑1 (MQ=25)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                               >  1:2757864/1‑91 (MQ=255)
             aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                               <  1:3716127/91‑1 (MQ=21)
              ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                              <  2:821260/91‑2 (MQ=21)
               tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                             <  1:2546247/91‑1 (MQ=255)
                gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                            <  1:2954204/91‑1 (MQ=255)
                 ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                           >  1:3900477/1‑91 (MQ=255)
                 ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                           >  2:3454975/1‑91 (MQ=255)
                    atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                                        <  2:2771394/91‑1 (MQ=255)
                    atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                                        <  2:114354/91‑1 (MQ=255)
                    atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                                        <  1:2110483/91‑1 (MQ=255)
                    atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                                       <  2:59580/91‑1 (MQ=255)
                      atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtttt                                                                      <  1:1564823/91‑1 (MQ=255)
                      atTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGAtttt                                                                      >  2:2687988/1‑91 (MQ=255)
                       ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                                     <  2:1511584/91‑1 (MQ=255)
                         ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                                   >  1:1576540/1‑91 (MQ=255)
                         ggAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTa                                                                   >  2:2866153/1‑91 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                                  >  1:539777/1‑91 (MQ=255)
                          gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                                  >  1:1002549/1‑91 (MQ=255)
                           aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                                 >  1:2383654/1‑91 (MQ=255)
                           aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTCTTCTaaa                                                                 >  1:1979356/1‑91 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                >  1:94610/1‑91 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                <  1:2017924/91‑1 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                <  2:1536748/91‑1 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                <  1:1642124/91‑1 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                >  2:2069944/1‑91 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                <  2:2650256/91‑1 (MQ=255)
                            aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                                >  1:2258841/1‑91 (MQ=255)
                              gCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAg                                                              >  2:400696/1‑91 (MQ=255)
                               cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                             >  2:94599/1‑91 (MQ=255)
                                  aCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtgt                                                          <  1:2905732/91‑1 (MQ=255)
                                    aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                                        <  1:813956/91‑1 (MQ=255)
                                        tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                    <  2:2290477/91‑1 (MQ=255)
                                         atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                                   >  1:174406/1‑91 (MQ=255)
                                         atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                                   >  2:834661/1‑91 (MQ=255)
                                         atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                                   <  2:1698476/91‑1 (MQ=255)
                                          tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAATGATa                                                  <  1:3126275/91‑1 (MQ=255)
                                             cGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAat                                               <  1:2654916/91‑1 (MQ=255)
                                               cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                             >  2:3705745/1‑91 (MQ=255)
                                                tAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattat                                            >  2:2003563/1‑91 (MQ=255)
                                                 aaaTAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattatt                                           >  2:3806918/1‑91 (MQ=255)
                                                  aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                                          <  1:3495047/91‑1 (MQ=255)
                                                    taatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttatttat                                        >  1:1897694/1‑91 (MQ=255)
                                                    taatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttatttat                                        <  2:884746/91‑1 (MQ=255)
                                                      atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTACAAAGGTGTTAAAGATAAATTATTTATaa                                      >  1:91068/1‑91 (MQ=255)
                                                       ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                                     >  2:3302691/1‑91 (MQ=255)
                                                       ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                                     <  1:1925059/91‑1 (MQ=255)
                                                          tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                                  >  1:2006019/1‑91 (MQ=255)
                                                          tCATGTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                                  <  1:344061/91‑1 (MQ=255)
                                                           cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                                 <  1:41483/91‑1 (MQ=255)
                                                            aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                                <  2:722183/91‑1 (MQ=255)
                                                             tttAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAcc                               <  2:725566/91‑1 (MQ=255)
                                                                aaaaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtt                            >  1:2155898/1‑91 (MQ=255)
                                                                  aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                          >  2:485642/1‑91 (MQ=255)
                                                                   aaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTg                         >  2:3400653/1‑91 (MQ=255)
                                                                     gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTCTTCAAGATAAATTATTTATAATGTTACCATTTTgag                       >  2:3919754/1‑91 (MQ=255)
                                                                     acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag                       >  1:2618320/2‑91 (MQ=37)
                                                                      cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                      >  1:3520440/1‑91 (MQ=255)
                                                                      cTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                      >  2:1130037/1‑91 (MQ=37)
                                                                       tAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGAt                     <  1:689032/91‑1 (MQ=255)
                                                                        aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg                    >  1:707616/1‑91 (MQ=255)
                                                                        aaCCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg                    >  1:2828213/1‑91 (MQ=37)
                                                                         acaaGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATCAATTATTTATAATGTTACCCTTTTTAGATGa                   >  2:1073120/4‑91 (MQ=255)
                                                                         aCCAGTCTCTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTGTGAGATGa                   >  2:1215734/1‑91 (MQ=255)
                                                                          ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTATCATTTTGAGAGGaa                  >  2:1104516/1‑91 (MQ=255)
                                                                          ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaa                  <  1:1322520/91‑1 (MQ=255)
                                                                           cAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATTaaa                 >  2:673090/1‑91 (MQ=255)
                                                                            aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg                <  1:376174/91‑1 (MQ=255)
                                                                             gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt               <  2:917727/91‑1 (MQ=255)
                                                                              tctatataaaGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg              >  1:3480615/5‑91 (MQ=255)
                                                                              tCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTg              <  2:3167425/91‑1 (MQ=255)
                                                                               cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa             >  2:1490331/1‑91 (MQ=255)
                                                                               cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa             <  1:3204294/91‑1 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAat          >  1:3257169/1‑91 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAat          >  2:2942327/1‑91 (MQ=255)
                                                                                  tatataGATTGGTTAGTTTTTTCTTGAATATATAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAat          <  2:1344853/91‑1 (MQ=255)
                                                                                   atataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAata         >  2:3963410/1‑91 (MQ=255)
                                                                                     agaGATTGGTTAGTTATTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTCTAATGTTACCATAGTGAGATGAAAGTGAAatatcg      >  2:3598307/3‑89 (MQ=255)
                                                                                       aGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGa     <  2:325203/91‑1 (MQ=255)
                                                                                       aGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTGATAATGTTACCATTTTGAGATGAAAGTGAAATATTGa     <  1:3342793/91‑1 (MQ=255)
                                                                                        gATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGat    <  1:306306/91‑1 (MQ=255)
                                                                                        gATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGat    <  1:1378484/91‑1 (MQ=255)
                                                                                         aTTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGata   >  1:2586345/1‑91 (MQ=255)
                                                                                          ttGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGatat  >  1:111827/1‑91 (MQ=255)
                                                                                          |                                                                                          
AGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTGATAT  >  USA300TCH1516_ALE/970600‑970780

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: