Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,107,417 G→A 5.8% intergenic (+52/‑116) pdhD → / → USA300HOU_RS05475 Dihydrolipoyl dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,107,4170GA 5.8% 977.8 / 33.3 290intergenic (+52/‑116)pdhD/USA300HOU_RS05475Dihydrolipoyl dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (156/117);  new base A (14/3);  total (170/120)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 302 total aligned reads displayed.

GCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACAA  >  USA300TCH1516_ALE/1107327‑1107507
                                                                                          |                                                                                          
gcagAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACg                                                                                            <  2:4361444/91‑1 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          >  2:3696011/1‑91 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          >  1:1651139/1‑91 (MQ=255)
  agAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGc                                                                                          <  1:4457866/91‑1 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         >  2:1805255/1‑91 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCt                                                                                         >  1:5091632/1‑91 (MQ=255)
   gAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGACTCAGTCATTAAAAGCCGTAGCATATGCTACAGCt                                                                                         <  2:3485802/91‑1 (MQ=39)
    aaaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCCGTAGCATATGCTACAGCtt                                                                                        <  2:2061647/91‑1 (MQ=255)
     aaaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCttt                                                                                       <  2:5099409/91‑1 (MQ=255)
      aaaGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttt                                                                                      <  2:4570666/91‑1 (MQ=255)
        aGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttttt                                                                                    >  2:1818813/1‑91 (MQ=255)
        aGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttttt                                                                                    >  1:2898250/1‑91 (MQ=255)
        aGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCtttttt                                                                                    >  2:5634327/1‑91 (MQ=255)
         gCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTg                                                                                   <  2:176889/91‑1 (MQ=255)
         gCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTg                                                                                   >  1:3057734/1‑91 (MQ=255)
         gCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGGCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTg                                                                                   <  1:4148529/91‑1 (MQ=255)
          cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                                  >  1:5376947/1‑91 (MQ=255)
          cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                                  >  2:4037420/1‑91 (MQ=255)
          cTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGt                                                                                  >  2:2780774/1‑91 (MQ=255)
            aTCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGttt                                                                                >  1:1094666/1‑91 (MQ=255)
             tCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGtttt                                                                               >  2:1636161/1‑91 (MQ=255)
               ggATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAg                                                                             >  2:3900286/1‑91 (MQ=255)
                    cccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTaaa                                                                        >  2:3107239/1‑91 (MQ=255)
                    cccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTaaa                                                                        <  2:2925653/91‑1 (MQ=255)
                     ccAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAg                                                                       >  1:5531079/1‑91 (MQ=255)
                      cAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGt                                                                      >  2:2434769/1‑91 (MQ=255)
                         tCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAt                                                                   <  1:3675667/91‑1 (MQ=255)
                         tCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAt                                                                   <  1:3863523/91‑1 (MQ=255)
                         tCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAt                                                                   >  2:103356/1‑91 (MQ=255)
                         tCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAAt                                                                   >  2:4442525/1‑91 (MQ=255)
                          ccATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATg                                                                  >  1:1892279/1‑91 (MQ=255)
                           cATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGt                                                                 <  1:1805255/91‑1 (MQ=255)
                              aaaaTGTAATAACTGATTATCTATAAAGATTCAGACATTAAAAGCTGTAGCATATTCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                                              >  2:4914732/3‑91 (MQ=255)
                              aCAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAg                                                              >  1:1381308/1‑91 (MQ=255)
                               cAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAgg                                                             >  1:3398986/1‑91 (MQ=255)
                                aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                                            <  1:2638582/91‑1 (MQ=255)
                                aaTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGa                                                            <  2:3588286/91‑1 (MQ=255)
                                 aTGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaa                                                           <  1:814489/91‑1 (MQ=255)
                                  tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  1:1726012/91‑1 (MQ=255)
                                  tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  1:3724287/91‑1 (MQ=255)
                                  tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  1:3519861/91‑1 (MQ=255)
                                  tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  2:3433794/91‑1 (MQ=255)
                                  tGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaaa                                                          <  2:1138167/91‑1 (MQ=255)
                                   gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTANGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAt                                                         <  1:2421749/91‑1 (MQ=255)
                                   gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGaagt                                                         >  1:5196183/1‑89 (MQ=255)
                                   gTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAAt                                                         <  1:227175/91‑1 (MQ=255)
                                      ataaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGa                                                      >  1:4536194/1‑91 (MQ=255)
                                       taaCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAt                                                     <  1:2670643/91‑1 (MQ=255)
                                         aCTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAttt                                                   >  1:3835428/1‑91 (MQ=255)
                                          cTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTg                                                  <  2:897951/91‑1 (MQ=255)
                                             aTTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTgaga                                               <  1:1206148/91‑1 (MQ=255)
                                              ttATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGat                                              <  2:2415472/91‑1 (MQ=255)
                                              ttATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGat                                              >  1:1054340/1‑91 (MQ=255)
                                              ttATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGat                                              >  1:4386019/1‑91 (MQ=255)
                                               tATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGAttt                                                   >  2:1544097/1‑85 (MQ=255)
                                                 tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                                           >  1:3174450/1‑91 (MQ=255)
                                                 tCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATc                                           <  2:4317795/91‑1 (MQ=255)
                                                   tatcAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGt                                         <  2:380062/91‑1 (MQ=255)
                                                    ataAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGtt                                        <  2:4972838/91‑1 (MQ=255)
                                                         gATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTataaa                                   >  2:5499383/1‑87 (MQ=255)
                                                            tCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACAtgtg                                <  1:2470400/91‑1 (MQ=255)
                                                              aGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAc                              >  2:759271/1‑91 (MQ=255)
                                                                tCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcat                            >  1:5445120/1‑91 (MQ=255)
                                                                 cATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatg                           <  2:541978/91‑1 (MQ=255)
                                                                  tttAAAAGCTGTAGCATATGCTACGGTTTTTTTGTTTTAGGTAAAGTAATGTATGGAAATTGATTTGAAATATCGTTAACATGTGAcatgc                          <  2:3635227/90‑1 (MQ=255)
                                                                    tAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcat                        <  2:3235052/91‑1 (MQ=255)
                                                                     aaaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg                       <  1:904744/91‑1 (MQ=255)
                                                                     aaaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGAcatgcatg                       <  1:3365445/91‑1 (MQ=255)
                                                                      aaaGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGt                      <  1:1268943/91‑1 (MQ=255)
                                                                      aaaGCCGTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGt                      >  2:122682/1‑91 (MQ=255)
                                                                        aGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTta                    <  1:358339/91‑1 (MQ=255)
                                                                         gCTGTAGCATATGCTACGGCTTTTTTTTTTTATGTTTAGTAATTTAAGGAAATTGATTTGATATATCGTTCACATGTGACATGCATGTtat                   >  2:5335643/1‑91 (MQ=255)
                                                                         gCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTTACATGCAAGTtat                   >  1:1485621/1‑91 (MQ=255)
                                                                         gCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtat                   >  1:2323474/1‑91 (MQ=255)
                                                                          ccgTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  >  1:4189800/3‑91 (MQ=255)
                                                                          cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  <  2:4633732/91‑1 (MQ=255)
                                                                          cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  <  2:984414/91‑1 (MQ=255)
                                                                          cTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtata                  >  2:3311785/1‑91 (MQ=255)
                                                                           tGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTtataa                 >  2:1883992/1‑90 (MQ=255)
                                                                            gTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGCTTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  2:4182557/1‑91 (MQ=255)
                                                                            gTAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACt                >  1:4298222/1‑91 (MQ=255)
                                                                             tAGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTa               >  2:642688/1‑91 (MQ=255)
                                                                              aGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAg              >  2:2882681/1‑91 (MQ=255)
                                                                              aGCATATGCTACAGCTTTTTTTTTTTTTTTTAAGTAATGTAAGGAAATTGATTTTAGATATCTTTAACATGTGTCATGCTTGTTTTACTAg              >  2:4239818/1‑91 (MQ=255)
                                                                              aGCATATGCTACAGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAg              >  1:2048370/1‑91 (MQ=255)
                                                                              aGCATATACTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAg              <  1:4652753/91‑1 (MQ=255)
                                                                               gCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGc             >  1:668732/1‑91 (MQ=255)
                                                                               gCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGc             >  1:4082576/1‑91 (MQ=255)
                                                                               gCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGc             >  2:2349042/1‑91 (MQ=255)
                                                                                cATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCg            >  1:5314512/1‑91 (MQ=255)
                                                                                cATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCg            >  2:5045051/1‑91 (MQ=255)
                                                                                  tatGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGAt          >  2:5237098/1‑91 (MQ=255)
                                                                                  tatGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGAt          <  1:435018/91‑1 (MQ=255)
                                                                                   atGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATg         <  2:2200080/91‑1 (MQ=255)
                                                                                     gCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCt       <  2:2286701/91‑1 (MQ=255)
                                                                                     gCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCt       <  2:250881/91‑1 (MQ=255)
                                                                                      cTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTa      <  1:3945054/91‑1 (MQ=255)
                                                                                      cTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTa      <  2:2665481/91‑1 (MQ=255)
                                                                                       tACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTaa     <  2:5196183/91‑1 (MQ=255)
                                                                                         cGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACa   >  1:5650593/1‑91 (MQ=255)
                                                                                          |                                                                                          
GCAGAAAAAGCTATCGGATACCCAATCCATACAATGTAATAACTGATTATCTATAAAGATTCAGTCATTAAAAGCTGTAGCATATGCTACGGCTTTTTTGTTTTAGGTAAAGTAATGTAAGGAAATTGATTTGAGATATCGTTAACATGTGACATGCATGTTATACTAGCGATGCTAACAA  >  USA300TCH1516_ALE/1107327‑1107507

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: