Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,934,780 A→T 6.8% N156K (AAT→AAA USA300HOU_RS09605 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,934,7800AT 6.8% 939.9 / 36.4 284N156K (AAT→AAAUSA300HOU_RS09605hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (135/128);  new base T (19/2);  total (154/130)
Fisher's exact test for biased strand distribution p-value = 4.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAATA  >  USA300TCH1516_ALE/1934690‑1934870
                                                                                          |                                                                                          
aCTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGaaa                                                                                            >  2:1825412/1‑91 (MQ=255)
    aaGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATACTTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGCAATTCa                                                                                        >  1:746368/1‑91 (MQ=255)
    aaGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCa                                                                                        >  1:1075790/1‑91 (MQ=255)
     agtagtTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaa                                                                                       <  1:3228059/91‑1 (MQ=255)
     aGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaa                                                                                       <  1:3593620/91‑1 (MQ=255)
     aGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTATTTGAATTTCaa                                                                                       <  1:477725/91‑1 (MQ=255)
       tATCTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaat                                                                                     <  1:3006656/91‑1 (MQ=255)
        aTTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaata                                                                                    <  1:678415/91‑1 (MQ=255)
         ttttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataa                                                                                   <  1:2817534/91‑1 (MQ=255)
         ttttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataa                                                                                   <  1:2064976/91‑1 (MQ=255)
          tttttCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaa                                                                                  <  2:3649985/91‑1 (MQ=255)
           ttttcttttcTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaat                                                                                 <  2:2409563/91‑1 (MQ=255)
           ttttcttttcTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCaaataaat                                                                                 <  1:984286/91‑1 (MQ=255)
             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                               >  2:692594/1‑91 (MQ=255)
             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                               <  2:1255511/91‑1 (MQ=255)
             ttcttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAatat                                                                               >  1:1747306/1‑91 (MQ=255)
               cttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTc                                                                             >  1:1721009/1‑91 (MQ=255)
                ttttcttCAGACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCt                                                                            <  2:1148486/91‑1 (MQ=255)
                ttttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCt                                                                            <  1:3501559/91‑1 (MQ=255)
                ttttctgCGACCGCTTCTTCCAGTTCCCAATAATTTTTGACTGATTCTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCt                                                                            <  2:584728/91‑1 (MQ=255)
                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCtt                                                                           >  1:1029937/1‑91 (MQ=255)
                 tttcttCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCtt                                                                           >  1:1028861/1‑91 (MQ=255)
                  ttcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCttt                                                                          <  2:371853/91‑1 (MQ=255)
                  ttcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCttt                                                                          <  1:1226735/91‑1 (MQ=255)
                  ttcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCttt                                                                          <  1:2804135/91‑1 (MQ=255)
                   tcttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTg                                                                         <  1:3204886/91‑1 (MQ=255)
                     ttcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCg                                                                       <  1:3693120/91‑1 (MQ=255)
                      tcAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTTAATTTCAAATAAATATTCTTTTCGt                                                                      >  2:2796385/1‑91 (MQ=255)
                       cAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTa                                                                     >  2:3370926/1‑91 (MQ=255)
                       cAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTa                                                                     >  2:2007327/1‑91 (MQ=255)
                       cAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTa                                                                     <  2:819131/91‑1 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   >  2:2618396/1‑91 (MQ=255)
                         aaCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtt                                                                   >  1:1860189/1‑91 (MQ=255)
                          aCTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAttt                                                                  >  2:3213419/1‑91 (MQ=255)
                           ctctTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTAtttt                                                                 <  2:1280297/91‑1 (MQ=255)
                            tcttctTCGAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTg                                                                <  2:3363466/91‑1 (MQ=255)
                             cttcttCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGa                                                               <  2:1594566/91‑1 (MQ=255)
                             cttcttCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGa                                                               <  2:1545622/91‑1 (MQ=255)
                              ttcttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAt                                                              >  2:1990973/1‑91 (MQ=255)
                                cttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGAttt                                                            >  2:3028739/1‑91 (MQ=255)
                                 ttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTc                                                           >  2:2416755/1‑91 (MQ=255)
                                 ttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTc                                                           >  2:1321569/1‑91 (MQ=255)
                                 ttcCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTc                                                           >  1:3741158/1‑91 (MQ=255)
                                   ccAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCAt                                                         <  1:2604751/91‑1 (MQ=255)
                                     aGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCt                                                       <  1:2800918/91‑1 (MQ=255)
                                            aataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaa                                                >  2:517225/1‑91 (MQ=255)
                                             ataatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaaa                                               >  2:522996/1‑91 (MQ=255)
                                              taatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaaaa                                              <  2:3810461/91‑1 (MQ=255)
                                              taatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTaaaa                                              <  2:2561067/91‑1 (MQ=255)
                                               aatTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAat                                             <  2:1624600/91‑1 (MQ=255)
                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTTATCTTATCTAAAATATc                                          >  2:2904598/1‑91 (MQ=255)
                                                  tattGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATc                                          <  2:1354937/91‑1 (MQ=255)
                                                   attGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCa                                         >  1:2195471/1‑91 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        >  1:2055608/1‑91 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAt                                        >  1:1677502/1‑91 (MQ=255)
                                                    ttGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGACTTCATCTTATCTAAAATATCAt                                        >  2:2949483/1‑91 (MQ=255)
                                                      gACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAttt                                      >  1:1157117/1‑91 (MQ=255)
                                                       aCTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCAtttt                                     <  2:2374763/91‑1 (MQ=255)
                                                        cTGATTGTTCTAAGTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTc                                    <  2:3201155/91‑1 (MQ=255)
                                                        cTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTc                                    >  1:1737555/1‑91 (MQ=255)
                                                          gATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGa                                  <  1:3325972/91‑1 (MQ=255)
                                                          gATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGa                                  <  2:2877974/91‑1 (MQ=255)
                                                            ttGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGata                                <  1:1945324/91‑1 (MQ=255)
                                                            ttGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGata                                <  1:1371871/91‑1 (MQ=255)
                                                             tGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGatat                               >  1:1825761/1‑91 (MQ=255)
                                                              gTTCTAAATTAGAGAATATTTTTTTGAAAGTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGAtata                              <  1:3113311/91‑1 (MQ=255)
                                                               ttCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGCTATaa                             >  2:983458/1‑91 (MQ=255)
                                                                tCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAc                            <  1:3371741/91‑1 (MQ=255)
                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                          <  2:3530301/91‑1 (MQ=255)
                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                          >  2:1112904/1‑91 (MQ=255)
                                                                  tAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCt                          >  2:275402/1‑91 (MQ=255)
                                                                   aaaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtt                         >  1:3037838/1‑91 (MQ=255)
                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                        <  1:1290790/91‑1 (MQ=255)
                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                        <  2:448287/91‑1 (MQ=255)
                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                        <  2:1677502/91‑1 (MQ=255)
                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                        >  1:2657861/1‑91 (MQ=255)
                                                                    aaTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCttt                        <  2:3063339/91‑1 (MQ=255)
                                                                     aTTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCtttt                       <  1:1917695/91‑1 (MQ=255)
                                                                      ttAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTa                      <  1:2309169/91‑1 (MQ=255)
                                                                       tAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTaa                     >  1:1154432/1‑91 (MQ=255)
                                                                           gaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCt                 <  1:3640160/91‑1 (MQ=255)
                                                                           gaATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATAGAACCTTTTAATTCt                 <  2:3292236/91‑1 (MQ=255)
                                                                            aaTATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCtt                >  1:345281/1‑91 (MQ=255)
                                                                               atttatttGAATTTCCAAAAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTTATTCTTTTa             >  1:1570548/1‑91 (MQ=38)
                                                                               atttatttGAATTTCAAAAAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa             >  1:842703/1‑91 (MQ=38)
                                                                               atTTTTTTGAAATTCAAATAAATAATCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTa             <  2:3610391/91‑1 (MQ=255)
                                                                                tttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAg            >  2:1217656/1‑91 (MQ=255)
                                                                                 ttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGa           <  2:2180973/91‑1 (MQ=255)
                                                                                 ttttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGa           >  2:2130428/1‑91 (MQ=255)
                                                                                 ttatttGAATTTCAAAAAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGa           >  1:2218762/4‑91 (MQ=39)
                                                                                  tttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGat          <  1:1263750/91‑1 (MQ=255)
                                                                                   ttttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGata         >  2:1131054/1‑91 (MQ=255)
                                                                                    tttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGatat        >  1:556026/1‑91 (MQ=255)
                                                                                     ttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATg       >  1:1458708/1‑91 (MQ=255)
                                                                                     ttGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATg       <  1:142052/91‑1 (MQ=255)
                                                                                      tGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCCTCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGa      <  1:2759510/91‑1 (MQ=255)
                                                                                      tGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAAAATTTTAATTCTTTTAGGTTTGa      <  2:189118/91‑1 (MQ=255)
                                                                                       gAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGaa     >  1:971578/1‑91 (MQ=255)
                                                                                       gAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGaa     <  2:2228061/91‑1 (MQ=255)
                                                                                        aaaTTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGaaa    <  1:2464788/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACTTAAGTATTTTTTCTTTTCTTCAAACTCTTCTTCCAGTTCCAAATAATTATTGACTGATTGTTCTAAATTAGAGAATATTTTTTTGAAATTCAAATAAATATTCTTTGCGTATTTTGATTTCATCTTATCTAAAATATCATTTTCGATATAACCTTTTAATTCTTTTAGATATGAAATA  >  USA300TCH1516_ALE/1934690‑1934870

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: