Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,148 T→C 11.7% intergenic (+55/‑76) USA300HOU_RS04665 → / → pepA_1 NADH dehydrogenase‑like protein/Cytosol aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE937,1480TC11.7% 681.1 / 70.3 256intergenic (+55/‑76)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (109/117);  new base C (4/26);  total (113/143)
Fisher's exact test for biased strand distribution p-value = 2.91e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

Only 100 of 262 total aligned reads displayed.

GGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTA  >  USA300TCH1516_ALE/937058‑937238
                                                                                          |                                                                                          
ggtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTACATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGt                                                                                            >  1:24686/1‑91 (MQ=255)
ggtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGt                                                                                            >  1:2059328/1‑91 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        <  1:910828/91‑1 (MQ=255)
    tCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCtt                                                                                        >  1:688289/1‑91 (MQ=255)
        tttAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgt                                                                                    <  1:3789447/91‑1 (MQ=255)
        tttAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgt                                                                                    <  2:160170/91‑1 (MQ=255)
          tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                  >  1:1263901/1‑91 (MQ=255)
          tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                  <  2:57843/91‑1 (MQ=255)
           ggCAGTCAAAAAAGGTAAATTCTAGGCTGTAAATAAAGTTTCAGCTAAGCACATAAAATTAAATTATTCATTTATAAAGCCCTTCAgtgac                                                                                 <  2:2471826/90‑3 (MQ=255)
           aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATGAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                 >  2:3902961/1‑91 (MQ=255)
           aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgac                                                                                 <  1:2474405/91‑3 (MQ=255)
           aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                 <  1:3619686/91‑1 (MQ=255)
           aGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                 <  2:3459499/91‑1 (MQ=255)
            gCATTCAAAAAAGGTAAATTCTAGTCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                <  2:1762865/91‑1 (MQ=255)
            gCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAgtgaca                                                                                <  1:456432/91‑4 (MQ=255)
             cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGacac                                                                               <  2:2430015/91‑1 (MQ=255)
              aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                              <  2:2658237/91‑1 (MQ=255)
              aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                              <  2:1467843/91‑1 (MQ=255)
              aTTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                              <  2:472449/91‑1 (MQ=255)
               ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTg                                                                             <  1:1949292/91‑1 (MQ=255)
                tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                            >  1:2403669/1‑91 (MQ=255)
                tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGa                                                                            <  1:2671698/91‑1 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                           >  1:2717264/1‑91 (MQ=255)
                 cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGaa                                                                           <  1:3927584/91‑1 (MQ=255)
                  aaaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                          >  2:1918051/1‑91 (MQ=255)
                   aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                         <  2:1196304/91‑1 (MQ=255)
                   aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                         <  1:1508815/91‑1 (MQ=255)
                   aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgg                                                                         <  2:308214/91‑1 (MQ=255)
                    aaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAggg                                                                        <  1:2579496/91‑1 (MQ=255)
                     aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                       <  1:2106721/91‑1 (MQ=255)
                     aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                       <  2:1858865/91‑1 (MQ=255)
                     aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                       >  1:2186674/1‑91 (MQ=255)
                     aaaGGTAAATTCAAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                       >  1:3049281/1‑91 (MQ=255)
                       aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACtt                                                                     <  1:155035/91‑1 (MQ=38)
                        ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAATTTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                    >  2:3475606/1‑91 (MQ=255)
                         gTAAAGTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAGAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTa                                                                   <  2:2643814/91‑1 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:3890396/1‑91 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:2985756/1‑91 (MQ=255)
                           aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                 >  1:1206838/1‑91 (MQ=255)
                            aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttg                                                                <  2:3232409/91‑1 (MQ=255)
                              ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                              <  1:1677493/91‑1 (MQ=255)
                               tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                             <  2:1034111/91‑1 (MQ=255)
                                 tAGGCTTTAAATAAAGTTTCAGCTAAGCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCGCTGAAGGGCTTTATTGTTGTa                                                           <  2:1044917/91‑1 (MQ=255)
                                  aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATATAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                          <  2:3632641/91‑1 (MQ=255)
                                  aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                          >  2:1501414/1‑91 (MQ=255)
                                  aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                          <  2:1638898/91‑1 (MQ=255)
                                  aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                          >  2:1871941/1‑91 (MQ=255)
                                   ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                         <  1:1730255/91‑1 (MQ=255)
                                     cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                       >  1:1737014/1‑91 (MQ=255)
                                     cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                       >  1:1154746/1‑91 (MQ=255)
                                      tttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAtt                                                      >  1:2901081/1‑91 (MQ=255)
                                       ttAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAttt                                                     <  1:3828666/91‑1 (MQ=255)
                                          aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                  >  1:3510991/1‑91 (MQ=255)
                                          aataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                  >  1:2555464/1‑91 (MQ=255)
                                               aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                             >  2:3921189/1‑91 (MQ=255)
                                               aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                             <  1:1091313/91‑1 (MQ=255)
                                               aGTTTCAGCTAAACACATAAAATTAAATTATACATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAGTTTGAAGTGAg                                             >  2:3024800/1‑91 (MQ=255)
                                               aGTTTCAGCTAAACACATAAAATTAAATTATACATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTAATTGTTGTACAGTTTGAAGTGAg                                             >  2:3024596/1‑91 (MQ=255)
                                                gTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGt                                            >  1:1121745/1‑91 (MQ=255)
                                                gTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGt                                            <  1:2133068/91‑1 (MQ=255)
                                                 tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                           <  1:2651219/91‑1 (MQ=255)
                                                  ttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTTTTGTTGTACAATTTGAAGTGAGTaa                                          >  2:1722832/1‑91 (MQ=255)
                                                   tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                         <  2:2161350/91‑1 (MQ=255)
                                                     aGCTAAACACATAACATTAAATTATTCATTTCTAAAGTCCTTCAGTGTCACTCAACGGCTTTATTGTTGTACAATTTGCAGTGAGTAACat                                       >  2:3355668/1‑91 (MQ=255)
                                                     aGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACat                                       <  1:2354701/91‑1 (MQ=255)
                                                        tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                    <  2:2784920/91‑1 (MQ=255)
                                                          aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTCACATATCaa                                  >  1:3788938/1‑91 (MQ=255)
                                                           acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                 <  1:2263980/91‑1 (MQ=255)
                                                           acacATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                 <  2:297277/91‑1 (MQ=255)
                                                              caTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                              <  1:1139211/91‑1 (MQ=255)
                                                               aTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTa                             >  1:412395/1‑91 (MQ=255)
                                                                tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                            >  1:563877/1‑91 (MQ=255)
                                                                tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                            >  2:245861/1‑91 (MQ=255)
                                                                 aaaaTTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGATGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAag                           >  1:3571491/1‑91 (MQ=255)
                                                                    attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAACTAaggag                        <  1:329081/91‑1 (MQ=255)
                                                                    attaaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATAGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAaggag                        <  1:3419759/91‑1 (MQ=255)
                                                                      taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                      >  2:3764023/1‑91 (MQ=255)
                                                                      taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                      >  1:1802090/1‑91 (MQ=255)
                                                                       aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                     <  1:1856815/91‑1 (MQ=255)
                                                                       aaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTa                     >  1:1326830/1‑91 (MQ=255)
                                                                        aattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTAc                    >  2:810161/1‑91 (MQ=255)
                                                                         attattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACa                   >  1:840827/1‑91 (MQ=255)
                                                                         attattCATTTATAAAGTCCTAAAGTGCCACTGAAGGGTTTTATTGTTGTACAATTTGAAGTGAGTAACATATGAAATCTAGGGAGGAAcg                   >  2:3524273/1‑90 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                <  2:525529/91‑1 (MQ=255)
                                                                            attCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                >  1:2823086/1‑91 (MQ=255)
                                                                              tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg              >  2:1792846/1‑91 (MQ=255)
                                                                              tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg              >  1:1471733/1‑91 (MQ=255)
                                                                               cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa             <  2:238878/91‑1 (MQ=255)
                                                                               cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa             <  1:1258590/91‑1 (MQ=255)
                                                                                aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTCATGTTGTACAATCTGAAGTGAGTAACATATAAAATCTAAGGCGGTACAGTATGaa            >  2:1143619/1‑91 (MQ=255)
                                                                                aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAATTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa            <  1:3343707/91‑1 (MQ=255)
                                                                                aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa            <  2:1195055/91‑1 (MQ=255)
                                                                                 tttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAt           >  2:1289431/1‑91 (MQ=255)
                                                                                    ataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtttt        <  1:2265345/91‑1 (MQ=255)
                                                                                     taAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta       >  2:2047581/1‑91 (MQ=255)
                                                                                     taAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTtta       >  1:1370591/1‑91 (MQ=38)
                                                                                      aaaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa      <  1:2260043/91‑1 (MQ=255)
                                                                                      aaaGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaa      >  1:3742382/1‑91 (MQ=38)
                                                                                         gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatt   <  1:403989/91‑1 (MQ=255)
                                                                                          tCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta  <  2:2778323/91‑1 (MQ=255)
                                                                                          |                                                                                          
GGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTA  >  USA300TCH1516_ALE/937058‑937238

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: