Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 614,082 G→A 11.5% S798S (TCG→TCA sdrC → Serine‑aspartate repeat‑containing protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE614,0820GA11.5% 156.0 / 2.1 76S798S (TCG→TCAsdrCSerine‑aspartate repeat‑containing protein C
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base A (8/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ATAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAG  >  USA300TCH1516_ALE/613946‑614221
                                                                                                                                        |                                                                                                                                            
aTAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACt                                                                                                                                          <  1:433246/141‑1 (MQ=255)
aTAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCAGACt                                                                                                                                          >  1:510480/1‑141 (MQ=255)
 aagagatacagattcaGACAGCGACGAAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTc                                                                                                                                         <  2:379362/137‑1 (MQ=255)
  aGCGATTCATATTCAGACAGCGACTCAGACTCAGACAACGATACAGACTCGGACAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGTCTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCa                                                                                                                                        <  1:189258/141‑1 (MQ=17)
   gCGATTCAGATTCAGACAGCGACTCAGATTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAg                                                                                                                                       <  2:385010/141‑1 (MQ=255)
    tgattcagatttaTGTAGTGATTCAGGCTCAGACAGCGATTCGGCCTCTGATAGAGACTCAGACTCAGATAGCGACTCAGATTCGTATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGTTTCGGACTCAGa                                                                                                                                      <  2:15977/140‑1 (MQ=14)
    cgagacagatgcagaCAGCGAGATAGAGTCAGACAGCGATACAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGa                                                                                                                                      <  2:503718/136‑1 (MQ=14)
    cGATTCAGATTCAGACAGCGACTCAGACTCCGAACGCGATGCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGa                                                                                                                                      <  1:709656/141‑1 (MQ=255)
    cGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGa                                                                                                                                      <  1:199314/141‑1 (MQ=255)
     gattcagatgcaGAAAGCGAATAAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGAc                                                                                                                                     <  2:175312/141‑1 (MQ=255)
             ttcagaCAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATCGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGAgattc                                                                                                                             >  1:448498/1‑141 (MQ=255)
             ttcagaCAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGAAAGTgattc                                                                                                                             >  1:225024/1‑141 (MQ=255)
             ttcagaCAGCGACTCAGACGCAGACACCGATTCAGACGCGCATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTgattc                                                                                                                             <  1:337939/141‑1 (MQ=255)
              ttgaaCAGCGAGTCGGACTCATACAGCGATTCAGACTCGGATAACGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTgattca                                                                                                                            <  2:1365/137‑1 (MQ=17)
              tcagaCAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCTGCTTCGAACAACCACACAGATTCAGATATCCACGCAGATTCAGATAGCGATTCGGACTCAGACAGTgattca                                                                                                                            <  1:357893/141‑1 (MQ=14)
              tcagaCAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTgattca                                                                                                                            >  1:3621/1‑141 (MQ=255)
               cagaCAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCCGACTCAGACAGTgatacag                                                                                                                           >  2:434033/1‑141 (MQ=255)
                    aGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGAATCAGATTCAGATAGCGACTCGGAATCAGAAAGCaaatcagaacaa                                                                                                                      >  2:123082/1‑137 (MQ=14)
                     gCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGattcagattcag                                                                                                                     >  1:336068/1‑141 (MQ=255)
                      cGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGAttcagattcaga                                                                                                                    >  1:302342/1‑141 (MQ=255)
                         ctcagactcagaCAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGCTAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTc                                                                                                                 >  1:228260/1‑141 (MQ=255)
                           cagactcagacAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGACTCAGACTCAGACAGAGATCCAGATTCAGACCCAg                                                                                                               >  2:264136/1‑141 (MQ=17)
                            agactcagacAGCGATTCAGACTCGGATAGTGACGCAGACTCAGATAGTGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGa                                                                                                              <  2:107209/141‑1 (MQ=255)
                              actcagacAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCAGATAGCGACTCAGACTCAGACAGAGATTCAGATTCAGACACAGAtc                                                                                                            >  2:511792/1‑140 (MQ=14)
                                   gataGCGATTCAGGCTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGAc                                                                                                       <  2:586715/138‑1 (MQ=255)
                                   gacAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGAc                                                                                                       <  1:409994/141‑1 (MQ=255)
                                     cAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTc                                                                                                     <  1:14800/141‑1 (MQ=255)
                                      aGCGATTCAGACTCGGATAGCGACTCAGACTCAGATGGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCa                                                                                                    >  1:361518/1‑141 (MQ=255)
                                      aGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCa                                                                                                    <  1:232931/141‑1 (MQ=255)
                                       gCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAg                                                                                                   >  1:257185/1‑141 (MQ=255)
                                            tCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCAGATTCAGACACCGATACAGACTCAGacaca                                                                                              >  2:360042/1‑137 (MQ=14)
                                             cAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCAGATAGCGATTCGGACTCAGAAAGCGATTCACATACAGACACAGATCGCGACTCacagccag                                                                                             >  2:749256/1‑134 (MQ=14)
                                              aGACTCGGATAGAGACTTAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGAGATTCAGACTCAGACAGTGATTCAGATTCAGACAAAGATAGCGACTCAGCTTCTGa                                                                                            >  2:660963/1‑141 (MQ=18)
                                              aGACTCAGATAGTGCGTGAAACTCAGATAGCGACTCGGAGTGGGATAGCGACTCAGATTCAGATAGCGATTGAGATTCAGATAGCGATTTTGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGa                                                                                            <  2:140040/141‑1 (MQ=17)
                                                 cTCGGATAGCGACTCAGACTCAGATAGCGACTCAGCTTCGGTTAGTGGAACAGAATCAGCTCATTTCTCAAATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAg                                                                                        >  1:357993/1‑141 (MQ=18)
                                                   cGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCCGATAGCGACTCAGATTCTGACAGCg                                                                                      >  1:118145/1‑141 (MQ=255)
                                                    tgATAGCGCGTCAGACTCACACAGCGACTCTGATGCGGATAGCGACTTAGACTCAGATAGCGAGGCAGATTCAGATAGCGATCCTGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGa                                                                                     <  2:18243/140‑1 (MQ=21)
                                                     aaTAGCGACGCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGAt                                                                                    <  1:196957/140‑1 (MQ=255)
                                                           gactcagactcaGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATtcagac                                                                              <  1:491162/141‑1 (MQ=255)
                                                             ctcagactcagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATTGCAACTAAGATTTTGATAGAGAGacagactc                                                                            >  1:429274/1‑141 (MQ=17)
                                                                agactcagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATtcagactcaga                                                                         <  1:359373/141‑1 (MQ=255)
                                                                 gactcagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATtcagactcagac                                                                        <  1:69090/141‑1 (MQ=255)
                                                                  actcagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGACtccgactccgacc                                                                       >  1:46998/1‑140 (MQ=255)
                                                                  actcagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGAGATTCAGACTCAGACa                                                                       >  1:89828/1‑141 (MQ=255)
                                                                     cagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCg                                                                    <  1:524893/141‑1 (MQ=255)
                                                                     cagaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTCAGATTCAGACACCGACACAGACTCAGATTCTGATCAACACTCAGACTCAGACacca                                                                    >  2:96031/1‑138 (MQ=14)
                                                                      agaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCga                                                                   >  1:160037/1‑141 (MQ=255)
                                                                      agaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCga                                                                   >  1:411776/1‑141 (MQ=255)
                                                                       gaTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAATGAATCAAATTCAAACACAGATCGAGACTCAGATTCAGATAACGATTCACACTCAGACAACgaa                                                                  >  2:370427/1‑140 (MQ=18)
                                                                        aTAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCgact                                                                 >  1:503590/1‑141 (MQ=255)
                                                                           gCGACGCAGATTCGGATAGCGACTCGGACGCAGATAACAATAAAGATGCGGATAGCGCATCAGCCAATTTCACGGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCgactcag                                                              <  1:355301/141‑1 (MQ=12)
                                                                             gACTGAGATTCAGACAAGGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACTGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGATAGCGATTCAGACTCAGACAGCgactcagac                                                            <  2:271472/141‑1 (MQ=18)
                                                                             gACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGAGCCAGATAGCGACTCAGACACTGACAGCGATTCAGACTCAGACAGCgactcagca                                                            >  2:634738/1‑139 (MQ=18)
                                                                                 cAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCgactccgactcaa                                                        >  1:577792/1‑140 (MQ=255)
                                                                                     ttgggATAGCGACTCAGACTCAGATACCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAg                                                    <  1:641041/138‑1 (MQ=255)
                                                                                     ttCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAg                                                    <  1:721501/141‑1 (MQ=255)
                                                                                     atccgggaGCGTCTCAGTCGCAGAGACGGACTCAGATTCAGAGAGCGATTCGGACGCAGACAGTGATTCAGAGTCAGACTCAGATAGCGACTCAGACTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAg                                                    <  2:393141/134‑1 (MQ=18)
                                                                                     atcagatagataCTAAGACTCAGATTCAGATTCAGATACAGATACCGATTCAGACTCAGACAGAGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAg                                                    <  2:291477/137‑1 (MQ=18)
                                                                                      tCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCAGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATACTGACAGAGATTCAGACTCAGACCGCGACtccgactcagacgac                                                   >  2:716822/1‑138 (MQ=21)
                                                                                       cGGATGCCGACTACGACTCAGACAGCGATTCAGATTCAGACAGCGACTCGGACTCAGACTCTGATGCAGACTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTg                                                  <  2:197691/141‑1 (MQ=17)
                                                                                       cGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCAGATTCAGACAGAGATAGCGACTCAGATACCGACACCGATTCAGACTAAGACAGAGACTCAGACCTAGACaaag                                                  >  2:198179/1‑138 (MQ=18)
                                                                                           tAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCAGACTCAGACAGCGATTCAGATTCAGACACAGATGCAGACTCAGATCCCGACACGGCTTCAGACGCCGACAGCGCCTCAGCCTCAGACCCTgattc                                              >  2:169152/1‑141 (MQ=14)
                                                                                            aGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGCCAGTgattca                                             >  1:195204/1‑141 (MQ=255)
                                                                                                        tcagaTAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGc                                 >  1:609255/1‑141 (MQ=255)
                                                                                                          agaTAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCATACAGCGACTCAGAATCAGACAGTGATTCAGATTCAGACAGCGa                               >  2:303987/1‑141 (MQ=255)
                                                                                                           aacagagattcagattcaGACAGAGACTCGGACTCCGACAGTGATTCAGATTCAGAATCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGAc                              <  2:361227/135‑1 (MQ=18)
                                                                                                               gagatacagatttaGATAACGATACAGACACCGCCTGGGACTAAGATAGAGACTCAGATAGCGACTCAGATTCTGACAGCGATACAGACTCAGACAGCGACTCAGACTCAGACAGTGAATCAGATTCAGACAGCGACtcag                          <  2:50140/139‑1 (MQ=14)
                                                                                                               gCGATTCAGATTCAGATAGCGATTCTGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCTACTCAGATTCTGACAGCGATTCAGACCCAGACAGCGACTGAGAGTGAGACAGTGATTCAGATTCAGACAGCGACtcag                          <  1:678562/141‑1 (MQ=21)
                                                                                                                cGATTGATATTCACATGGAGATTCGAACTCAGACACTGATTCAGATACAGACTCAGATAGCGACTCAGACTCTGACAGCGACTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACtcaga                         <  2:291596/141‑1 (MQ=17)
                                                                                                                         ttcagatAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatag                <  1:411581/141‑1 (MQ=255)
                                                                                                                              atAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTGAGCTAGCGACTCAGATTCAGGCAGCGATTCAGAGTCAGAGAGCGACTCAGAGTCGGACAGGGATTAAGATCAAGACAGCGACTCAAATgcagacagcgcgt           >  2:547431/1‑138 (MQ=18)
                                                                                                                               caGCGACTCAGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACCCAGAACCACCTTTGCTTGCAGACAACGACTCAGACtcagatagcgactc          >  2:7546/2‑141 (MQ=17)
                                                                                                                                     ttCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactc    >  1:135263/1‑141 (MQ=255)
                                                                                                                                       cGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTGAGGCTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcag  >  1:276446/1‑141 (MQ=255)
                                                                                                                                       cGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcag  <  1:97020/141‑1 (MQ=255)
                                                                                                                                       cGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcacacagcgactcagactcac  >  2:443262/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                            
ATAGCGATTCAGATTCAGACAGCGACTCAGACTCAGACAGCGATTCAGACTCGGATAGCGACTCAGACTCAGATAGCGACTCAGATTCGGATAGCGACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAG  >  USA300TCH1516_ALE/613946‑614221

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