Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 937,158 T→A 5.5% intergenic (+65/‑66) 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,1580TA 5.5% 1090.6 / 40.6 310intergenic (+65/‑66)USA300HOU_RS04665/pepA_1NADH dehydrogenase‑like protein/Cytosol aminopeptidase
Reads supporting (aligned to +/- strand):  ref base T (142/151);  new base A (8/9);  total (150/160)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 320 total aligned reads displayed.

TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248
                                                                                          |                                                                                          
tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAgtgt                                                                                            >  1:4414841/1‑91 (MQ=255)
  gCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                          >  1:5091304/1‑91 (MQ=255)
  gCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                          >  2:4797043/1‑91 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCAc                                                                                         >  2:3141150/1‑91 (MQ=255)
   cATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCAc                                                                                         <  1:1821482/91‑1 (MQ=255)
     ttCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTg                                                                                       >  2:290268/1‑91 (MQ=255)
     atCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTg                                                                                       <  2:1893602/90‑1 (MQ=255)
      tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGa                                                                                      >  2:3894171/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     >  2:1977018/1‑91 (MQ=255)
       cAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGaa                                                                                     <  1:4585659/91‑1 (MQ=255)
        aaaaaaGGTAAATTCTAGTCTTTAAATAAAGTTTCAGCTAAACACCTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                                    <  2:5668659/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:4456534/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:338576/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAgg                                                                                   <  1:1007992/91‑1 (MQ=255)
         aaaaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgg                                                                                   <  1:1653364/91‑1 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 >  2:5254438/1‑91 (MQ=255)
           aaaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 >  1:1929664/1‑91 (MQ=255)
            aaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCt                                                                                <  1:946158/91‑1 (MQ=255)
            aaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCt                                                                                >  2:4543937/1‑91 (MQ=255)
            aaGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCt                                                                                >  1:346628/1‑91 (MQ=255)
             aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                               >  2:5831079/1‑91 (MQ=255)
             aGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCtt                                                                               >  1:2977770/1‑91 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAGTTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              <  2:3939429/91‑1 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              >  2:2761159/1‑91 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCttt                                                                              >  1:3088935/1‑91 (MQ=255)
              ggTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAGGACttt                                                                              <  2:1031391/91‑1 (MQ=38)
               gTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTa                                                                             >  1:1729968/1‑91 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                 >  2:1770595/1‑85 (MQ=255)
                 aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAtt                                                                           >  1:706642/1‑91 (MQ=255)
                   aTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                         <  2:2191103/91‑1 (MQ=255)
                   aTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgt                                                                         <  1:601180/91‑1 (MQ=255)
                    ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                        >  2:2286738/1‑91 (MQ=255)
                    ttCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGacac                                                                                         >  2:501311/1‑74 (MQ=255)
                     tCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgttg                                                                       <  2:1729715/91‑1 (MQ=255)
                     tCTAGGCATTAAATAAAGTATCGGCTAAACACATAAAATTAAATTATTCATTTGTAAAGACCTTCGGTGTCACTGAAGGGCTTTAttgttg                                                                       <  2:3606399/91‑1 (MQ=255)
                       tAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTa                                                                     <  1:3701957/91‑1 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    <  1:5033195/91‑1 (MQ=255)
                        aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTAc                                                                    >  2:2021451/1‑91 (MQ=255)
                         ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                   >  1:3417235/1‑91 (MQ=255)
                         ggCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                   >  2:3328213/1‑91 (MQ=255)
                           cTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                 >  1:558433/1‑91 (MQ=255)
                                aataaaGTTTAAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGaa                                                            >  2:4861703/1‑91 (MQ=255)
                                 ataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGATTTATTGTTGTACAATTTGAAg                                                           <  2:4466844/91‑1 (MQ=255)
                                   aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                         >  1:248413/1‑91 (MQ=255)
                                   aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATATGAAGTg                                                         <  2:299973/91‑1 (MQ=255)
                                    aaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGa                                                        <  2:1588124/91‑1 (MQ=255)
                                     aGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTAACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                       >  2:4239558/1‑91 (MQ=255)
                                      gTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGt                                                      >  1:5741651/1‑91 (MQ=255)
                                       tttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTa                                                     >  1:2134273/1‑91 (MQ=255)
                                        ttCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTaa                                                    <  1:4182241/91‑1 (MQ=255)
                                         tCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                   >  1:2320771/1‑91 (MQ=255)
                                          cAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                  >  2:345669/1‑91 (MQ=255)
                                           aGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACat                                                 >  1:1140649/1‑91 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  2:1096703/91‑1 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  1:5340114/91‑1 (MQ=255)
                                             cTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACatat                                               <  2:2382026/91‑1 (MQ=255)
                                              tAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                              <  1:1843089/91‑1 (MQ=255)
                                               aaaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCa                                             <  2:770111/91‑1 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            <  1:1041389/91‑1 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            <  1:2930909/91‑1 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            >  1:4034014/1‑91 (MQ=255)
                                                aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaa                                            >  2:194976/1‑91 (MQ=255)
                                                    caTAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                        <  1:4930259/91‑1 (MQ=255)
                                                    caTAAAATTAAATCATTCATTTATAAACTCATTCAGTGTCACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCt                                        >  1:4662047/1‑91 (MQ=255)
                                                      tAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                      >  1:1088565/1‑91 (MQ=255)
                                                      tAAAATTAAATTATTCATTTATAAAGTCCTTCAGAGTCACTGAAGAGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTaa                                      >  2:4112989/1‑91 (MQ=255)
                                                       aaaaTTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAag                                     <  1:303758/91‑1 (MQ=255)
                                                            taaattaTTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGt                                >  1:3283581/1‑91 (MQ=255)
                                                               attattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACa                             >  1:2030867/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTTACTGAAGGGCTTTATTGATTTACAACTCAACGCGAGTACCATATCAACTCTAACGATGTACAt                            >  2:4338719/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            >  2:910287/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            >  2:2632277/1‑91 (MQ=255)
                                                                ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAt                            <  1:4797043/91‑1 (MQ=255)
                                                                   ttCATTTATAAAGTCCTTCAGTGTCACTGACGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTAt                         >  2:5389262/1‑91 (MQ=255)
                                                                    tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                        <  2:3290755/91‑1 (MQ=255)
                                                                    tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                        <  2:618653/91‑1 (MQ=255)
                                                                    tCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                        <  2:520431/91‑1 (MQ=255)
                                                                     cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGa                       <  2:665724/91‑1 (MQ=255)
                                                                      aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                      <  2:2544316/91‑1 (MQ=255)
                                                                      aTTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGaa                      <  1:84163/91‑1 (MQ=255)
                                                                       tttATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAt                     >  2:3687149/1‑91 (MQ=255)
                                                                         tataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   <  1:538526/91‑1 (MQ=255)
                                                                         aataAAGCCCTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAttt                   >  1:1238473/2‑91 (MQ=37)
                                                                          ataAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtttt                  >  1:3495845/1‑91 (MQ=255)
                                                                             aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa               <  1:3141150/91‑1 (MQ=255)
                                                                             aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaa               <  1:1545410/91‑1 (MQ=255)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              <  2:421409/91‑1 (MQ=255)
                                                                              aGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaat              <  1:2567511/91‑1 (MQ=255)
                                                                               gTCCTTCAGTGTCACTGAAGGGCATTATTGTTGTAAAATTTGAAGAGAGTAACATATCAAATATAAGGAGGTAAATTATGAATTttaaatt             >  2:3199284/1‑91 (MQ=255)
                                                                                  cTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaaa          <  2:1925088/91‑1 (MQ=255)
                                                                                  cTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaattaaa          >  1:2464032/1‑91 (MQ=255)
                                                                                   ttCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTtaaattaaat         >  1:230618/1‑91 (MQ=255)
                                                                                    tCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa        <  2:4602279/91‑1 (MQ=255)
                                                                                    tCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa        >  1:2879004/1‑91 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCACATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       <  1:4493748/91‑1 (MQ=255)
                                                                                     cAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa       <  1:2721351/91‑1 (MQ=255)
                                                                                      aGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      <  1:3835719/91‑1 (MQ=255)
                                                                                      aGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAac      >  2:3750288/1‑91 (MQ=255)
                                                                                         gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca   >  1:5843697/1‑91 (MQ=255)
                                                                                         gtCACTGAAGGGCTATATTGTTGTACAATTTGAAGTGAGTAACATATAAAAACTAAGGAGGTACATTATGAATTTTAAATTAAATAacaca   >  2:3790850/1‑91 (MQ=255)
                                                                                          |                                                                                          
TAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACAC  >  USA300TCH1516_ALE/937068‑937248

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: