Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 970,686 T→A 10.5% intergenic (+32/‑297) 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,6860TA10.5% 743.0 / 55.3 239intergenic (+32/‑297)USA300HOU_RS04805/USA300HOU_RS04810Ureidoglycolate lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (107/105);  new base A (4/21);  total (111/128)
Fisher's exact test for biased strand distribution p-value = 1.15e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 252 total aligned reads displayed.

ATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTG  >  USA300TCH1516_ALE/970596‑970776
                                                                                          |                                                                                          
atgaAGTTAAAGTGACTATTGATAATAATGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTtatat                                                                                            >  2:1874061/1‑91 (MQ=255)
  gaAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAg                                                                                          >  1:53901/1‑91 (MQ=255)
   aaGTTAAACTGCCTCTTGCGAATATTAGACCGCTGACAACTCATACCCCTAAATAATTGTCATTTAAAAAGCTAACCAATCTATATAAAGa                                                                                         <  2:3884991/91‑1 (MQ=255)
      ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTg                                                                                      <  2:273359/91‑1 (MQ=255)
      ttAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGactg                                                                                      <  2:882885/91‑4 (MQ=37)
       tAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTgg                                                                                     <  2:1080964/91‑1 (MQ=37)
        aaaGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGt                                                                                    <  2:2422600/91‑1 (MQ=38)
          aGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTa                                                                                  <  1:864309/91‑1 (MQ=255)
           gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                                 >  1:3092351/1‑91 (MQ=255)
           gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAg                                                                                 >  1:335328/1‑91 (MQ=255)
           gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAATCTATATAAAGACTGGTTAg                                                                                 <  1:351578/91‑1 (MQ=37)
           gTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAg                                                                                 <  2:387889/91‑1 (MQ=25)
            tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                                <  2:1453459/91‑1 (MQ=255)
            tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGt                                                                                <  2:2778534/91‑1 (MQ=255)
            tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAGAGACTGGTTAgc                                                                                <  2:2193292/91‑2 (MQ=37)
            tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgc                                                                                <  1:3279213/91‑2 (MQ=25)
            tGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAgc                                                                                <  1:3709693/91‑2 (MQ=25)
            cgACTATTAAAAATTTTGAAACGCTGACACCTTATGTCGCTAAATAGTTATCATTTAAAAAGCTAACCGATCTATATAGAGACTGGTTAgc                                                                                <  2:1807975/90‑2 (MQ=255)
              aCTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttt                                                                              >  2:1322520/1‑91 (MQ=255)
               cTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttt                                                                             >  1:966611/1‑91 (MQ=255)
                tATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGttttt                                                                            <  1:2206747/91‑1 (MQ=255)
                 aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                           >  2:3237152/1‑91 (MQ=255)
                 aTTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGtttttt                                                                           >  2:1527875/1‑91 (MQ=255)
                  ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                          >  1:3407630/1‑91 (MQ=255)
                  ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                          <  1:1919495/91‑1 (MQ=255)
                  ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTc                                                                          >  2:2107215/1‑91 (MQ=255)
                  ttGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAACTAACCAATCTATATAAAGACTGGTTAGcttttta                                                                          <  2:821260/91‑2 (MQ=21)
                   tGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCt                                                                         <  1:2546247/91‑1 (MQ=255)
                    gATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCtt                                                                        <  1:2954204/91‑1 (MQ=255)
                     ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                       >  1:937154/1‑91 (MQ=255)
                     ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                       >  1:3900477/1‑91 (MQ=255)
                     ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                       >  1:1300989/1‑91 (MQ=255)
                     ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                       >  2:3454975/1‑91 (MQ=255)
                     ataataTTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTg                                                                       <  1:1642237/91‑1 (MQ=255)
                        atatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCtt                                                                    <  1:2110483/91‑1 (MQ=255)
                         tatTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCttt                                                                   <  1:1438076/91‑1 (MQ=255)
                           ttGGAACGCTGACAACTTATATTACTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                                 >  2:3573968/1‑91 (MQ=255)
                           ttGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTc                                                                 <  2:1436448/91‑1 (MQ=255)
                              gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                              >  1:1002549/1‑91 (MQ=255)
                              gAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaa                                                              >  2:1123965/1‑91 (MQ=255)
                               aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                             >  1:2383654/1‑91 (MQ=255)
                               aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                             >  1:1979676/1‑91 (MQ=255)
                               aaCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaa                                                             <  1:3487768/91‑1 (MQ=255)
                                aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                            <  1:505257/91‑1 (MQ=255)
                                aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                            >  1:2258841/1‑91 (MQ=255)
                                aCGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaa                                                            <  2:2650256/91‑1 (MQ=255)
                                 cGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTaaaaa                                                           >  2:401054/1‑91 (MQ=255)
                                   cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                         >  2:636550/1‑91 (MQ=255)
                                   cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                         >  2:94599/1‑91 (MQ=255)
                                   cTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAgg                                                         <  1:3699551/91‑1 (MQ=255)
                                    tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                        <  1:2799770/91‑1 (MQ=255)
                                    tGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgt                                                        >  1:2667628/1‑91 (MQ=255)
                                      aCAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGgtgt                                                      <  1:2905732/91‑1 (MQ=255)
                                       aaaCTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGACTGGTTAGTTTTTCCTAGCTTTTCTACAAAGGTGtt                                                     >  2:3378102/2‑91 (MQ=255)
                                            tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGa                                                <  2:2290477/91‑1 (MQ=255)
                                            tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAGGTGTTAAAGAt                                               <  1:1576199/91‑1 (MQ=255)
                                            tataTCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGATTTAATAAAAAGGTGTTAAAGa                                                <  2:2788100/91‑1 (MQ=255)
                                             atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                               >  2:834661/1‑91 (MQ=255)
                                             atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                               >  1:174406/1‑91 (MQ=255)
                                             atatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGAt                                               <  2:1698476/91‑1 (MQ=255)
                                              tatCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATa                                              >  1:1375185/1‑91 (MQ=255)
                                                 cGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAat                                           >  2:2938014/1‑91 (MQ=255)
                                                   cTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAatta                                         >  2:3705745/1‑91 (MQ=255)
                                                    tAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattat                                        >  2:2003563/1‑91 (MQ=255)
                                                     aaaTAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattatt                                       <  1:3528849/91‑1 (MQ=255)
                                                     aaaTAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAattatt                                       >  2:3806918/1‑91 (MQ=255)
                                                      aataatTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAAttattt                                      <  1:1273176/91‑1 (MQ=255)
                                                          atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                                  <  2:1135886/91‑1 (MQ=255)
                                                          atTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATaa                                  >  1:1135145/1‑91 (MQ=255)
                                                           ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                                 >  2:3302691/1‑91 (MQ=255)
                                                           ttATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAAt                                 <  1:1925059/91‑1 (MQ=255)
                                                             aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                               <  1:3463173/91‑1 (MQ=255)
                                                             aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                               >  1:2065057/1‑91 (MQ=255)
                                                             aTCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGt                               <  2:1248528/91‑1 (MQ=255)
                                                              tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                              >  1:2006019/1‑91 (MQ=255)
                                                              tCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGtt                              <  2:295816/91‑1 (MQ=255)
                                                               cATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTa                             <  1:41483/91‑1 (MQ=255)
                                                                aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                            <  2:722183/91‑1 (MQ=255)
                                                                aTTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTAc                            >  2:2812915/1‑91 (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)
                                                                    aaaaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtt                        >  1:3145054/1‑91 (MQ=255)
                                                                      aaaGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCAtttt                      >  2:485642/1‑91 (MQ=255)
                                                                       aaGCTAACCAGTCTTTATATAGATTGGTTAGTATTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTg                     >  2:1571328/1‑91 (MQ=255)
                                                                        aGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTga                    <  1:2640596/91‑1 (MQ=255)
                                                                          cTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTgaga                  >  2:3171720/1‑91 (MQ=255)
                                                                            aaCCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg                <  1:1586018/91‑1 (MQ=255)
                                                                            aaCCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg                >  1:2828213/1‑91 (MQ=37)
                                                                            aaCCAATCTATATAAAGACTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATg                >  1:2609351/1‑91 (MQ=38)
                                                                             aCCAGTCTTTATCTAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGCTGTTAAAGATAAATTATTTATAATGTTACCTTTTTGAGATGa               >  2:392773/1‑91 (MQ=255)
                                                                              ccAGTCTTTATATAGCTTGGCTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTCTCATGTTACCATTTTGAGCTGaa              >  2:1103051/1‑91 (MQ=255)
                                                                               cAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATTaaa             >  2:673090/1‑91 (MQ=255)
                                                                                aGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAg            <  1:376174/91‑1 (MQ=255)
                                                                                 gTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGt           <  2:917727/91‑1 (MQ=255)
                                                                                   cTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGa         >  2:12039/1‑91 (MQ=255)
                                                                                   cTTTATATAGATTGGTTAGTTTTTTCTTGCTCTTCTAAAAAGGTGTTAAAGATCAATTATTTATAATGTTACCATTTTGAGATGAACGTGa         >  2:1491527/1‑91 (MQ=255)
                                                                                       atataGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAata     >  2:3963410/1‑91 (MQ=255)
                                                                                          taGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTg  <  1:3255682/91‑1 (MQ=255)
                                                                                          |                                                                                          
ATGAAGTTAAAGTGACTATTGATAATATTGGAACGCTGACAACTTATATCGCTAAATAATTATCATTTAAAAAGCTAACCAGTCTTTATATAGATTGGTTAGTTTTTTCTTGCTTTTCTAAAAAGGTGTTAAAGATAAATTATTTATAATGTTACCATTTTGAGATGAAAGTGAAATATTG  >  USA300TCH1516_ALE/970596‑970776

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: