Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,104,520 T→A 11.8% L309* (TTA→TAA)  pdhB → Pyruvate dehydrogenase E1 component subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,104,5200TA11.8% 401.3 / 37.5 127L309* (TTA→TAA) pdhBPyruvate dehydrogenase E1 component subunit beta
Reads supporting (aligned to +/- strand):  ref base T (68/44);  new base A (7/8);  total (75/52)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 130 total aligned reads displayed.

AACGTGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTATTTTAGTC  >  USA300TCH1516_ALE/1104430‑1104610
                                                                                          |                                                                                          
aaCGTGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGtt                                                                                            >  1:2851927/1‑91 (MQ=255)
  cGTGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTAc                                                                                          >  1:1006574/1‑91 (MQ=255)
   gTGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAAcc                                                                                         <  1:1656981/91‑1 (MQ=255)
    tGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACca                                                                                        <  1:3787459/91‑1 (MQ=255)
    tGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACca                                                                                        <  1:3786286/91‑1 (MQ=255)
    tGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAACca                                                                                        <  1:1928630/91‑1 (MQ=255)
     gCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaa                                                                                       >  1:1916966/1‑91 (MQ=255)
     gCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAACcaa                                                                                       <  1:3999939/91‑1 (MQ=255)
     gCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAACcaa                                                                                       <  1:3600933/91‑1 (MQ=255)
     gCAATCCCTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAACcaa                                                                                       <  1:2078323/91‑1 (MQ=255)
      cAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaa                                                                                      <  1:1491104/91‑1 (MQ=255)
      cAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaa                                                                                      <  1:1491210/91‑1 (MQ=255)
      cAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaa                                                                                      <  1:1492383/91‑1 (MQ=255)
      cAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaa                                                                                      <  1:2274211/91‑1 (MQ=255)
       aaTCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaac                                                                                     <  1:40976/91‑1 (MQ=255)
       aaTCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTAACcaaac                                                                                     <  1:1428977/91‑1 (MQ=255)
          ccTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACcaaacaaa                                                                                  <  1:1959784/91‑1 (MQ=255)
            tttCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACaaaaa                                                                                <  1:147859/91‑1 (MQ=255)
             ttCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAAt                                                                               >  1:2203588/1‑91 (MQ=255)
             ttCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAAt                                                                               >  1:2000912/1‑91 (MQ=255)
             ttCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAAt                                                                               >  1:107501/1‑91 (MQ=255)
              tCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATg                                                                              <  1:3824955/91‑1 (MQ=255)
              tCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATg                                                                              >  1:489396/1‑91 (MQ=255)
               cATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGa                                                                             <  1:766373/91‑1 (MQ=255)
               cATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGa                                                                             <  1:3573690/91‑1 (MQ=255)
                aTTAGAAGCACCTATTGGAAGCGTTGCAGCAGCAGATACAATTTCTCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGAc                                                                            >  1:2040296/1‑91 (MQ=255)
                aTTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGAc                                                                            >  1:2040371/1‑91 (MQ=255)
                    gAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGAcatca                                                                        <  1:1973570/91‑1 (MQ=255)
                    gAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGAcatca                                                                        <  1:1974005/91‑1 (MQ=255)
                     aaGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGAcatcat                                                                       >  1:1379315/1‑91 (MQ=255)
                       gCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCg                                                                     <  1:615882/91‑1 (MQ=255)
                        cACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGa                                                                    >  1:3505153/1‑91 (MQ=255)
                        cACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGa                                                                    >  1:1760958/1‑91 (MQ=255)
                        cACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGa                                                                    <  1:3377479/91‑1 (MQ=255)
                         aCCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGaa                                                                   >  1:1900201/1‑91 (MQ=255)
                             aTTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAg                                                               >  1:820398/1‑91 (MQ=255)
                              ttGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGc                                                              <  1:405367/91‑1 (MQ=255)
                              ttGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGc                                                              >  1:2098475/1‑91 (MQ=255)
                                ggAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCaa                                                            >  1:3530043/1‑91 (MQ=255)
                                ggAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCaa                                                            <  1:2395077/91‑1 (MQ=255)
                                  aaGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCaaaa                                                          >  1:2640654/1‑91 (MQ=255)
                                   agagTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAg                                                         <  1:757937/91‑1 (MQ=255)
                                    gagTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGa                                                        >  1:768372/1‑91 (MQ=255)
                                    gagTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGa                                                        >  1:195120/1‑91 (MQ=255)
                                    gagTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGa                                                        >  1:199212/1‑91 (MQ=255)
                                     agTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGaa                                                       >  1:3691515/1‑91 (MQ=255)
                                      gTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGaaa                                                      >  1:3922223/1‑91 (MQ=255)
                                         gcagcagcaGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACtt                                                   >  1:1004839/1‑91 (MQ=255)
                                         gcagcagcaGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACtt                                                   >  1:3701685/1‑91 (MQ=255)
                                          cagcagcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACttt                                                  >  1:1829612/1‑91 (MQ=255)
                                          cagcagcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACttt                                                  >  1:3063135/1‑91 (MQ=255)
                                            gcagcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAg                                                <  1:2827943/91‑1 (MQ=255)
                                             cagcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGa                                               >  1:180668/1‑91 (MQ=255)
                                               gcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAAt                                             >  1:1235277/1‑91 (MQ=255)
                                               gcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAAt                                             >  1:1852439/1‑91 (MQ=255)
                                               gcagATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAAt                                             >  1:3661140/1‑91 (MQ=255)
                                                  gATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAACCAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAAtttt                                          >  1:3292837/1‑91 (MQ=255)
                                                  gATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAAtttt                                          >  1:3373155/1‑91 (MQ=255)
                                                   aTACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTa                                         <  1:3429318/91‑1 (MQ=255)
                                                    tACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTaa                                        >  1:3429668/1‑91 (MQ=255)
                                                    tACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTaa                                        >  1:750909/1‑91 (MQ=255)
                                                      cAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATa                                      <  1:2686345/91‑1 (MQ=255)
                                                      cAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATa                                      >  1:2019507/1‑91 (MQ=255)
                                                       aaTTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATAc                                     <  1:2337084/91‑1 (MQ=255)
                                                       aaTTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATAc                                     <  1:1657313/91‑1 (MQ=255)
                                                       aaTTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATAc                                     <  1:2337507/91‑1 (MQ=255)
                                                       aaTTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATAc                                     <  1:2337631/91‑1 (MQ=255)
                                                        aTTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACa                                    >  1:3418557/1‑91 (MQ=255)
                                                             tCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTa                               <  1:1808155/91‑1 (MQ=255)
                                                              ccATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTaa                              >  1:564804/1‑91 (MQ=255)
                                                               cATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTaaa                             >  1:2025074/1‑91 (MQ=255)
                                                                 ttCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAg                           >  1:3524091/1‑91 (MQ=255)
                                                                      tCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAAc                      <  1:3795052/91‑1 (MQ=255)
                                                                        aaGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGa                    >  1:3166908/1‑91 (MQ=255)
                                                                         aGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGaa                   >  1:3201620/1‑91 (MQ=255)
                                                                         aGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGaa                   >  1:826075/1‑91 (MQ=255)
                                                                         aGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGaa                   >  1:3051077/1‑91 (MQ=255)
                                                                         aGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGaa                   >  1:3198374/1‑91 (MQ=255)
                                                                          gCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAg                  >  1:1996093/1‑91 (MQ=255)
                                                                              aaaaTGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAg              <  1:2587373/91‑1 (MQ=255)
                                                                              aaaaTGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAg              <  1:525716/91‑1 (MQ=255)
                                                                               aaaTGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGc             >  1:38130/1‑91 (MQ=255)
                                                                               aaaTGTTTGGTTACCAAACAAAAACGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGc             >  1:1774511/1‑91 (MQ=255)
                                                                                aaTGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCg            <  1:896579/91‑1 (MQ=255)
                                                                                 aTGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGt           >  1:2981147/1‑91 (MQ=255)
                                                                                  tGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTa          <  1:3041007/91‑1 (MQ=255)
                                                                                  tGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTa          <  1:1839341/91‑1 (MQ=255)
                                                                                   gTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAt         >  1:2265071/1‑91 (MQ=255)
                                                                                   gTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAt         >  1:2264837/1‑91 (MQ=255)
                                                                                    tttGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtt        >  1:2928840/1‑91 (MQ=255)
                                                                                    tttGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtt        <  1:1178697/91‑1 (MQ=255)
                                                                                     ttGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAttt       >  1:1095627/1‑91 (MQ=255)
                                                                                      tGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtttt      <  1:1422404/91‑1 (MQ=255)
                                                                                      tGGTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTCGAATTTTAATACATTTTAAAAGTTAACGACGTTAGCGTAtttt      >  1:4142947/1‑91 (MQ=255)
                                                                                      tGGTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtttt      >  1:3467355/1‑91 (MQ=255)
                                                                                      tGGTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtttt      >  1:2945820/1‑91 (MQ=255)
                                                                                      tGGTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtttt      >  1:711704/1‑91 (MQ=255)
                                                                                      tGGTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTAtttt      >  1:761138/1‑91 (MQ=255)
                                                                                       ggTAACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTATTTTa     >  1:1263976/1‑91 (MQ=255)
                                                                                          tACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTATTTTAGtc  <  1:2810900/91‑1 (MQ=255)
                                                                                          |                                                                                          
AACGTGCAATCCTTTCATTAGAAGCACCTATTGGAAGAGTTGCAGCAGCAGATACAATTTATCCATTCACTCAAGCTGAAAATGTTTGGTTACCAAACAAAAATGACATCATCGAAAAAGCAAAAGAAACTTTAGAATTTTAATACATTTTAAAAGTTAACGAAGTTAGCGTATTTTAGTC  >  USA300TCH1516_ALE/1104430‑1104610

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