Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,681 T→A 11.0% 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,6810TA11.0% 746.5 / 64.6 247intergenic (+27/‑302)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (108/111);  new base A (4/23);  total (113/134)
Fisher's exact test for biased strand distribution p-value = 7.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

Only 100 of 249 total aligned reads displayed.

GTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAA  >  USA300TCH1516_ALE/970593‑970770
                                                                                        |                                                                                         
gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                         >  2:965334/1‑91 (MQ=255)
gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                         <  1:3518309/91‑1 (MQ=255)
gTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTta                                                                                         <  1:2441699/91‑1 (MQ=255)
     gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                    >  2:38921/1‑91 (MQ=255)
        gTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGAtt                                                                                 >  1:2998457/1‑91 (MQ=255)
          tAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTgg                                                                               <  2:1080964/91‑1 (MQ=37)
           aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGt                                                                              <  2:2805044/91‑1 (MQ=255)
           aaaGGGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGt                                                                              <  1:506560/91‑1 (MQ=38)
            aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                             <  1:2822343/91‑1 (MQ=255)
            aaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGtt                                                                             >  2:2785694/1‑91 (MQ=255)
              gTGACTATTGATAATATTGGAACGCTGGCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                           <  2:1654213/91‑1 (MQ=37)
              gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                           <  1:78529/91‑1 (MQ=255)
              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)
                gACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtt                                                                         <  2:922951/91‑1 (MQ=255)
                 aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                        >  2:1322520/1‑91 (MQ=255)
                  cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                       >  1:966611/1‑91 (MQ=255)
                  cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                       >  2:640482/1‑91 (MQ=255)
                    aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                     >  2:3237363/1‑91 (MQ=255)
                    aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                     >  1:2757864/1‑91 (MQ=255)
                    aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                     <  1:3716127/91‑1 (MQ=21)
                    aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttt                                                                     <  1:287658/91‑1 (MQ=21)
                     ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                    <  2:821281/91‑2 (MQ=18)
                      tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                   <  1:250578/91‑1 (MQ=255)
                        ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                 >  2:3454975/1‑91 (MQ=255)
                        ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                 >  1:1300989/1‑91 (MQ=255)
                           atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                              <  2:2771394/91‑1 (MQ=255)
                           atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                              <  1:2110483/91‑1 (MQ=255)
                           atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                              >  2:724974/1‑91 (MQ=255)
                           atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                             <  2:59580/91‑1 (MQ=255)
                            tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                             <  1:1438076/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                                                           <  2:1511584/91‑1 (MQ=255)
                                 gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                        >  1:539777/1‑91 (MQ=255)
                                 gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                        >  2:1123965/1‑91 (MQ=255)
                                 gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                        >  1:1002549/1‑91 (MQ=255)
                                  aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                       >  1:1979676/1‑91 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      <  2:1536748/91‑1 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      <  1:2556317/91‑1 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      <  1:1642124/91‑1 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      >  1:2258545/1‑91 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      <  1:2017924/91‑1 (MQ=255)
                                   aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                      >  2:2069944/1‑91 (MQ=255)
                                     gCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAg                                                    >  2:400696/1‑91 (MQ=255)
                                           aaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTa                                              <  2:1841127/91‑1 (MQ=255)
                                               tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                          <  2:2290477/91‑1 (MQ=255)
                                                atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                         >  1:174406/1‑91 (MQ=255)
                                                 tatCGCTAAATAATTATCATTTAAAAAGCTAACCTGTCTTAATATAGATTGGATAGTTTTTTCTTGATTTTCTAAAAAGGTGTTAAAGATa                                        >  2:3678730/1‑91 (MQ=255)
                                                 tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAATGATa                                        <  1:3126275/91‑1 (MQ=255)
                                                 aatCGCTAAATAATTATCATTTAAAAAGCTAACCCGTCGTTAGATAGATTGGTAAGTTTTTTCTTGCTTGGCTAACAAGGTGTTAAAGATa                                        <  2:46450/90‑1 (MQ=255)
                                                    cGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAat                                     <  1:2654916/91‑1 (MQ=255)
                                                      cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                   >  1:2339294/1‑91 (MQ=255)
                                                      cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                   >  2:3705745/1‑91 (MQ=255)
                                                       tAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattat                                  >  2:2003563/1‑91 (MQ=255)
                                                        aaaTAATTATCATTTAAAAAGCTATCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattatt                                 <  1:940290/91‑1 (MQ=255)
                                                        aaaTAATTATCATTTAAAAAGCTAACCAGTATTTATATAGATTGGTTAGTTTTTTAAAAAAAAAAAAAAAAGGTGTTAAAGATAAattatt                                 <  2:1767795/91‑1 (MQ=255)
                                                         aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                                <  1:3496103/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)
                                                             atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                            >  1:1135145/1‑91 (MQ=255)
                                                             atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                            <  2:1135886/91‑1 (MQ=255)
                                                              ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                           <  1:1925059/91‑1 (MQ=255)
                                                              ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                           >  2:3302691/1‑91 (MQ=255)
                                                                aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                         <  2:1248528/91‑1 (MQ=255)
                                                                aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                         <  1:3761969/91‑1 (MQ=255)
                                                                aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                         <  1:3463173/91‑1 (MQ=255)
                                                                 tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                        >  1:2006019/1‑91 (MQ=255)
                                                                 tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                        <  2:295816/91‑1 (MQ=255)
                                                                 tCATTTAAAAAGCTAACCAGTCTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                       <  2:869383/91‑1 (MQ=255)
                                                                 tCATGTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                        <  1:344061/91‑1 (MQ=255)
                                                                   aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                      <  2:721578/91‑1 (MQ=255)
                                                                   aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                      <  1:3313688/91‑1 (MQ=255)
                                                                      tAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAt                   >  2:2838866/1‑91 (MQ=255)
                                                                          aaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTg               >  2:3400653/1‑91 (MQ=255)
                                                                            gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag             <  1:674160/91‑1 (MQ=255)
                                                                            gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag             >  1:464450/1‑91 (MQ=255)
                                                                            gCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag             <  1:292693/91‑1 (MQ=255)
                                                                            acTAACCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgag             >  1:2618320/2‑91 (MQ=37)
                                                                             cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga            >  1:3520440/1‑91 (MQ=255)
                                                                             cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga            >  2:3171720/1‑91 (MQ=255)
                                                                             cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTGAAAGATAAATTATTTATAATGTTACCATTTTgaga            >  2:779575/1‑91 (MQ=255)
                                                                               aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg          >  1:707616/1‑91 (MQ=255)
                                                                               aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg          <  1:1586018/91‑1 (MQ=255)
                                                                                aCCAGTCTTTATCTAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGCTGTTAAAGATAAATTATTTATAATGTTACCTTTTTGAGATGa         >  2:392773/1‑91 (MQ=255)
                                                                                 ccAGTCTTTATATAGCTTGGCTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTCTCATGTTACCATTTTGAGCTGaa        >  2:1103051/1‑91 (MQ=255)
                                                                                 ccAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaa        <  1:1322520/91‑1 (MQ=255)
                                                                                  ccgTCTTCATATAGATTGGCTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGaaa       >  2:2492179/3‑91 (MQ=255)
                                                                                   aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg      <  2:1449490/91‑1 (MQ=255)
                                                                                    gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt     <  2:3830335/91‑1 (MQ=255)
                                                                                      cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa   <  1:3203848/91‑1 (MQ=255)
                                                                                      cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa   >  2:1490331/1‑91 (MQ=255)
                                                                                      cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa   >  2:12039/1‑91 (MQ=255)
                                                                                       tttATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATCATTTATAATGTTACCATTTTGAGATGAAAGTGaa  >  1:3944085/1‑91 (MQ=255)
                                                                                        |                                                                                         
GTGATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAA  >  USA300TCH1516_ALE/970593‑970770

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: