Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,987,374 A→C 6.7% intergenic (‑437/+340) USA300HOU_RS09895 ← / ← USA300HOU_RS09905 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,987,3740AC 6.7% 993.7 / 6.2 295intergenic (‑437/+340)USA300HOU_RS09895/USA300HOU_RS09905hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (125/143);  new base C (22/5);  total (147/148)
Fisher's exact test for biased strand distribution p-value = 5.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 308 total aligned reads displayed.

ACAGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCTTAAA  >  USA300TCH1516_ALE/1987284‑1987463
                                                                                          |                                                                                         
aCAGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCaaa                                                                                           >  2:1921276/1‑91 (MQ=255)
aCAGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCaaa                                                                                           >  1:4548457/1‑91 (MQ=255)
 cAGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCaaaa                                                                                          >  1:4155677/1‑91 (MQ=255)
  aGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACAt                                                                                         >  2:1725006/1‑91 (MQ=255)
  aGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAAt                                                                                         <  1:4821212/91‑1 (MQ=255)
      aTTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTg                                                                                     <  1:517870/91‑1 (MQ=255)
      aTTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTg                                                                                     <  1:2122340/91‑1 (MQ=255)
       ttGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTga                                                                                    <  2:1049133/91‑1 (MQ=255)
        tGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTgag                                                                                   <  2:2699148/91‑1 (MQ=255)
         gCAAATTGGGGAAACGGGCCACAAAGTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTgaga                                                                                  <  2:5221844/91‑1 (MQ=255)
         gCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAgagaaa                                                                                  >  2:2582600/1‑89 (MQ=255)
         gCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTgaga                                                                                  >  1:4080860/1‑91 (MQ=255)
           aaaTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGGGGCCCCAAAATAGAGAAATc                                                                                >  2:4306438/1‑91 (MQ=255)
           aaaTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGAGAAATc                                                                                >  2:2255829/1‑91 (MQ=255)
            aaTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCg                                                                               <  1:4529309/91‑1 (MQ=255)
            aaTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCg                                                                               <  1:1327445/91‑1 (MQ=255)
             aTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCgg                                                                              >  1:5412557/1‑91 (MQ=255)
             aTTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGAAAAATCgg                                                                              >  2:5572640/1‑91 (MQ=255)
               tGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGat                                                                            <  1:2314857/91‑1 (MQ=255)
                ggggAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatt                                                                           >  1:2397615/1‑91 (MQ=255)
                ggggAAACGCGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatt                                                                           <  2:5353059/91‑1 (MQ=255)
                 gggAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGattt                                                                          <  1:3628924/91‑1 (MQ=255)
                 gggAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGattt                                                                          >  1:761103/1‑91 (MQ=255)
                  ggAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttc                                                                         >  1:1712887/1‑91 (MQ=255)
                  ggAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttc                                                                         <  2:1960083/91‑1 (MQ=255)
                    aaaCGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACATAGTGAGATCGGatctcta                                                                       >  2:235420/1‑91 (MQ=255)
                    aaaCGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttcta                                                                       <  1:1392992/91‑1 (MQ=255)
                     aaCGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatgtctaa                                                                      >  1:3895726/1‑91 (MQ=255)
                      aCGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttctaat                                                                     <  2:569310/91‑1 (MQ=255)
                      aCGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttctaat                                                                     <  2:3182213/91‑1 (MQ=255)
                       cGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGAGAGATCGGatctctaatt                                                                    >  2:2641783/1‑91 (MQ=255)
                          gCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGatttctaatttct                                                                 <  2:3315398/91‑1 (MQ=255)
                            cacaAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACATAGAGAGATCGGATCTCTAATTTCTAc                                                               >  2:3224237/1‑91 (MQ=255)
                             acaAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACa                                                              >  1:3080066/1‑91 (MQ=255)
                                aaCTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGAc                                                           <  2:1067841/91‑1 (MQ=255)
                                 aCTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACa                                                          >  1:570471/1‑91 (MQ=255)
                                  cTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAt                                                         <  2:2727839/91‑1 (MQ=255)
                                  cTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAt                                                         <  2:3072649/91‑1 (MQ=255)
                                  cTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAt                                                         <  1:3985371/91‑1 (MQ=255)
                                   tCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAtt                                                        <  1:3868131/91‑1 (MQ=255)
                                   tCAGAAACTGGTGGAAAGTCAGCTGAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAtt                                                        <  1:5350966/91‑1 (MQ=255)
                                    cAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTg                                                       >  1:4271257/1‑91 (MQ=255)
                                    cAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTg                                                       <  2:738317/91‑1 (MQ=255)
                                    cAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTAGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTg                                                       >  1:5582089/1‑91 (MQ=255)
                                     aGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACATAGAGAGATCGGATTTCTAATTTCTACAGACATTgt                                                      >  2:399123/1‑90 (MQ=255)
                                     aGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGc                                                      >  1:2831474/1‑91 (MQ=255)
                                      gAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCa                                                     <  1:1877526/91‑1 (MQ=255)
                                      gAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCa                                                     <  2:1034217/91‑1 (MQ=255)
                                      gAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGAGAGATCGGATTTCTAATTTCTACAGACATTGCa                                                     >  2:1617465/1‑91 (MQ=255)
                                       aaaCTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCaa                                                    >  1:1018131/1‑91 (MQ=255)
                                        aaCTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCaaa                                                   >  1:379267/1‑91 (MQ=255)
                                        aaCTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCaaa                                                   <  1:2536815/91‑1 (MQ=255)
                                            ggtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACATAGAGAGATCGGATCTCTAATTTCTACAGACAATGTAAATTAg                                               >  2:3691733/1‑91 (MQ=17)
                                            ggtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAg                                               <  1:459576/91‑1 (MQ=255)
                                            ggtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAg                                               <  1:1200677/91‑1 (MQ=255)
                                            ggtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAg                                               >  1:3579621/1‑91 (MQ=255)
                                            ggtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATCTCTAATTTCTACAGACATTGCAAATTAg                                               >  2:2999206/1‑91 (MQ=255)
                                             gtggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAgg                                              <  1:826305/91‑1 (MQ=255)
                                             gaggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAACTAGAGAGACCGGATTTCTAATTTCTACAGACATTGTAAGTTCgg                                              >  2:4716796/3‑91 (MQ=21)
                                              tggGAAGTCAGGTTAAAATAACATGCAAGTTGGTGAGGCTCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAggg                                             <  2:348036/91‑1 (MQ=25)
                                              tggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAACATAGAGAAATCGGATCTCTAATTTCTACAGACAATGAAAATTAggg                                             >  2:4079321/1‑91 (MQ=17)
                                              tggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAggg                                             >  1:4244796/1‑91 (MQ=255)
                                               ggAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGa                                            >  1:510406/1‑91 (MQ=255)
                                               ggAAAGGCAGCATAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGa                                            <  2:3785240/91‑1 (MQ=255)
                                                gAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGaa                                           <  2:1267055/91‑1 (MQ=255)
                                                     tCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACggg                                      <  1:3951079/91‑1 (MQ=255)
                                                      cAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGc                                     >  1:1652612/1‑91 (MQ=255)
                                                      cAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGc                                     >  1:368803/1‑91 (MQ=255)
                                                       aGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGcc                                    <  1:3506870/91‑1 (MQ=255)
                                                        gCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGCAACGGGCca                                   >  2:598497/1‑91 (MQ=255)
                                                        gCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCca                                   <  1:1366588/91‑1 (MQ=255)
                                                         cTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCcac                                  >  2:612032/1‑91 (MQ=255)
                                                          tttaaaTAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCcaca                                 <  2:3841978/88‑1 (MQ=255)
                                                          ttAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCcaca                                 <  1:4960833/91‑1 (MQ=255)
                                                              aaTAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACt                             <  2:2834467/91‑1 (MQ=255)
                                                                tAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCa                           <  2:352453/91‑1 (MQ=255)
                                                                tAACATGCAAGGTGGCGAGGCCCAAAAATAGTGAGAGCGGATTTCAGATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACGCa                           <  2:5467353/91‑1 (MQ=21)
                                                                  aCATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGa                         <  1:5454234/91‑1 (MQ=255)
                                                                   cATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGaa                        >  1:3136600/1‑91 (MQ=255)
                                                                   cATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGACCACAAACTCAGaa                        >  1:4683647/1‑91 (MQ=255)
                                                                    aTGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAg                       >  1:3824959/1‑91 (MQ=255)
                                                                     tGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGt                      <  1:2682278/91‑1 (MQ=255)
                                                                       cAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTtg                    >  1:1123217/1‑91 (MQ=255)
                                                                         aGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTtggt                  >  1:1954858/1‑91 (MQ=255)
                                                                         aGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTtggt                  >  1:4308627/1‑91 (MQ=255)
                                                                          gTTGGCGAGGCCACAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTtggtg                 <  1:2255039/91‑1 (MQ=255)
                                                                           ttGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTtggtgg                <  1:3454019/91‑1 (MQ=255)
                                                                              gCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGaaa             <  1:1921276/91‑1 (MQ=255)
                                                                               cGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAg            >  2:1508220/1‑91 (MQ=255)
                                                                               cGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACAATGCAAAATAGGGAAACGAACCACAAACTCAGAAGTTGGTGGAAAg            >  2:5099703/1‑91 (MQ=38)
                                                                                ggagCCCCAAAATAGTGAAATTGGATTTCAAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGt           <  2:851125/88‑1 (MQ=14)
                                                                                    ccccAAAATAGTGGGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGc       <  2:4472459/91‑1 (MQ=255)
                                                                                    ccccAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGc       <  1:3239364/91‑1 (MQ=255)
                                                                                     cccAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCt      <  1:853956/91‑1 (MQ=255)
                                                                                     cccAAAATAGTGAGATCGGATTTCCAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCt      <  2:5251538/91‑1 (MQ=255)
                                                                                     cccAAAATAGTGAAATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCt      <  2:1412560/91‑1 (MQ=255)
                                                                                      ccAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCtt     <  1:4402803/91‑1 (MQ=255)
                                                                                       cAACATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCTTa    <  2:1365138/91‑1 (MQ=255)
                                                                                       cAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCTTa    >  2:4131714/1‑91 (MQ=255)
                                                                                        aaaaTAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCTTaa   >  1:3362587/1‑91 (MQ=255)
                                                                                          |                                                                                         
ACAGACATTGCAAATTGGGGAAACGGGCCACAAACTCAGAAACTGGTGGAAAGTCAGCTTAAAATAACATGCAAGTTGGCGAGGCCCCAAAATAGTGAGATCGGATTTCTAATTTCTACAGACATTGCAAATTAGGGAAACGGGCCACAAACTCAGAAGTTGGTGGAAAGTCAGCTTAAA  >  USA300TCH1516_ALE/1987284‑1987463

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: