Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,313,113 A→T 10.7% intergenic (+17/+280) rny_1 → / ← USA300HOU_RS06485 Ribonuclease Y/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,313,1130AT10.7% 634.6 / 61.1 216intergenic (+17/+280)rny_1/USA300HOU_RS06485Ribonuclease Y/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (106/86);  new base T (3/20);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 1.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-01

Only 100 of 230 total aligned reads displayed.

ATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTG  >  USA300TCH1516_ALE/1313023‑1313202
                                                                                          |                                                                                         
aTGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTcaca                                                                                           <  2:1772636/91‑1 (MQ=255)
aTGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTcaca                                                                                           <  2:1772465/91‑1 (MQ=255)
 tGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACaa                                                                                          <  1:3660922/91‑1 (MQ=255)
  gAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACaaa                                                                                         <  1:3868046/91‑1 (MQ=255)
   aagtACAAGATCCTGGTCATATCAAGGGGACAGTTGTTCGAGAGACTAGAGCAGGAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAt                                                                                        <  2:3846359/88‑1 (MQ=255)
   aaTTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAt                                                                                        <  1:1170614/91‑1 (MQ=255)
    aTTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                       <  2:3491798/91‑1 (MQ=255)
    aTTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAAtt                                                                                       >  1:3356937/1‑91 (MQ=255)
    aTTACAATATCCTGGTCATATCAAGGGGACAGTTGTTCGGTAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtt                                                                                       <  2:1480570/91‑1 (MQ=255)
     ttACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAttt                                                                                      <  1:932673/91‑2 (MQ=255)
       aCAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACTAAtttgt                                                                                    <  1:1555446/91‑4 (MQ=255)
       aCAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGt                                                                                    >  2:1170482/1‑91 (MQ=255)
       aCAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGt                                                                                    <  2:3187379/91‑1 (MQ=255)
        cAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTg                                                                                   >  1:1328870/1‑91 (MQ=255)
           tatCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAgg                                                                                >  1:337048/1‑91 (MQ=255)
            atCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAggg                                                                               <  2:569873/91‑1 (MQ=255)
              ccTGGGCAGATCAAGGTGACAGTTGTTCGTGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAg                                                                             <  2:1040474/91‑1 (MQ=255)
                tGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCt                                                                           <  1:926806/91‑1 (MQ=255)
                 ggTCATTTCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtt                                                                          >  1:3076687/1‑91 (MQ=255)
                 ggTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtt                                                                          >  2:1726897/1‑91 (MQ=255)
                 ggTCATATCAAGGTGACAGTTGTTCGAGAGACCAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtt                                                                          >  2:3883303/1‑91 (MQ=255)
                  gTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCttt                                                                         >  2:3105296/1‑91 (MQ=255)
                   tCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtttt                                                                        >  1:741091/1‑91 (MQ=255)
                     atatCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCtttttt                                                                      <  1:2090927/91‑1 (MQ=255)
                       atCAAGGTGACAGTTGTTCGAGGGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTaa                                                                    <  2:3177645/91‑1 (MQ=255)
                        tCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAg                                                                   >  1:3251789/1‑91 (MQ=255)
                        tCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAg                                                                   >  2:1965965/1‑91 (MQ=255)
                           aGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTg                                                                >  1:1908488/1‑91 (MQ=255)
                            ggTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGt                                                               <  2:1890685/91‑1 (MQ=255)
                              tGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAg                                                             >  2:2539234/1‑91 (MQ=255)
                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  2:3388772/91‑1 (MQ=255)
                               gACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGt                                                            <  1:375335/91‑1 (MQ=255)
                                aCAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTc                                                           >  1:1382997/1‑91 (MQ=255)
                                 cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                          >  1:3606348/1‑91 (MQ=255)
                                 cAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCt                                                          <  2:2984850/91‑1 (MQ=255)
                                  aGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAACTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTCTAAGTTGTAGTCtt                                                         >  1:46870/1‑91 (MQ=255)
                                   gttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTa                                                        >  1:2211746/1‑91 (MQ=255)
                                   gttgttCGAGAGACTAGAGCAGCAGAATATGCGAAATAATTTTTGTCTCGCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTa                                                        <  2:3001752/91‑1 (MQ=255)
                                    ttgttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTaa                                                       <  1:1177066/91‑1 (MQ=255)
                                      gttCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATc                                                     >  2:592643/1‑91 (MQ=255)
                                      gttCGAGAGACGAGAGCAGTAGTCTATGCGAAATAATTTTTGTCTCCCTGACAAATTAGGGAGGGAGCTTTTTTAAGTTGTAGTCTTAATc                                                     <  2:3012508/91‑1 (MQ=255)
                                        tCGAGAGACTAGAGCAGTAGAATATGCGAAATCATTTCTGTCTCGCTCAAAAACTCGTGAGGGGGTTTTTTTAAGTTGTAGTCTTAATCtt                                                   >  2:3880137/1‑90 (MQ=255)
                                        tCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACACATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTa                                                   >  2:208207/1‑91 (MQ=255)
                                        tCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTa                                                   >  1:3454668/1‑91 (MQ=255)
                                        tCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTa                                                   >  1:1581906/1‑91 (MQ=255)
                                        tCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTa                                                   <  2:2092580/91‑1 (MQ=255)
                                         cGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAg                                                  >  1:1189961/1‑91 (MQ=255)
                                         cGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAg                                                  >  1:2409231/1‑91 (MQ=255)
                                             agaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAg                                              >  1:632526/1‑91 (MQ=255)
                                              gaCTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGa                                             <  1:1880186/91‑1 (MQ=255)
                                                cTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACa                                           >  2:2626669/1‑91 (MQ=255)
                                                  agagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGc                                         >  1:1799933/1‑91 (MQ=255)
                                                  agagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGc                                         <  2:2387494/91‑1 (MQ=255)
                                                   gagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACACATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCa                                        >  2:2812910/1‑91 (MQ=255)
                                                   gagCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCa                                        >  1:1594855/1‑91 (MQ=255)
                                                    agCCGTAGAATCTGCGAACTAATTCTTGTCTCACTCCCACATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTCATCTAGTTAGACAGCAc                                       >  2:1363290/1‑91 (MQ=255)
                                                    agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                       <  1:2336685/91‑1 (MQ=255)
                                                    agCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCAc                                       >  1:3005490/1‑91 (MQ=255)
                                                     gCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACt                                      <  2:647228/91‑1 (MQ=255)
                                                     gCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACt                                      >  1:533670/1‑91 (MQ=255)
                                                       aGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTCAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACttt                                    <  1:3963393/91‑1 (MQ=255)
                                                        gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                   <  2:3322471/91‑1 (MQ=255)
                                                        gTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTa                                   <  1:2036319/91‑1 (MQ=255)
                                                          aGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATc                                 >  2:3032176/1‑91 (MQ=255)
                                                           gAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCg                                <  1:3189128/91‑1 (MQ=255)
                                                           gAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCg                                >  1:988482/1‑91 (MQ=255)
                                                            aaTATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCgg                               >  1:1549544/1‑91 (MQ=255)
                                                             atatGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGt                              <  1:3857723/91‑1 (MQ=255)
                                                             atatGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGt                              <  2:2881267/91‑1 (MQ=255)
                                                             atatGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGt                              >  2:2669769/1‑91 (MQ=255)
                                                             atatGCGAAATAATTTTTGTATCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGt                              <  2:2705502/91‑1 (MQ=255)
                                                                   gagataatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  2:3709859/88‑1 (MQ=255)
                                                                   gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  2:1437013/91‑1 (MQ=255)
                                                                   gAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAAc                        <  1:1070455/91‑1 (MQ=255)
                                                                     aataatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACta                      >  1:1015758/1‑91 (MQ=255)
                                                                      ataatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACtat                     <  1:1019456/91‑1 (MQ=255)
                                                                       taatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACtata                    <  1:2427578/91‑1 (MQ=255)
                                                                        aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                   >  2:1815861/1‑91 (MQ=255)
                                                                        aatTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTatat                   <  2:1404029/91‑1 (MQ=255)
                                                                         atTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACtactgt                  <  1:706330/91‑5 (MQ=255)
                                                                         atTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATAtt                  >  2:335174/1‑91 (MQ=255)
                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  2:203945/1‑91 (MQ=255)
                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  1:1068471/1‑91 (MQ=255)
                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  2:1759787/1‑91 (MQ=255)
                                                                           ttttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  1:2819497/1‑91 (MQ=255)
                                                                           ttttGTCTCCCTCACAAATCAGCGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaa                >  2:459817/1‑91 (MQ=255)
                                                                            tttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaaa               >  2:3056145/1‑91 (MQ=255)
                                                                            tttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTaaa               >  1:261405/1‑91 (MQ=255)
                                                                             ttGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAAc              <  1:2498614/91‑1 (MQ=255)
                                                                               gTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACag            >  1:746822/1‑91 (MQ=255)
                                                                               gTCTCCCTCACAAAATAGTGAGGTAGCTTTTTTACGTTGTAGTCTTAATCTAGCTAGACAGCCCTTTACCGGTAATAACTATATTAACCag            >  2:85416/1‑91 (MQ=255)
                                                                               agctcCCTCACTAATTTGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACag            >  1:607216/3‑91 (MQ=39)
                                                                                 ctcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagta          <  2:90438/91‑1 (MQ=255)
                                                                                 ctcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagta          >  1:3893540/1‑91 (MQ=255)
                                                                                  tcCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtag         >  1:1373701/1‑91 (MQ=255)
                                                                                   cccTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACagtagt        <  2:1291296/91‑1 (MQ=255)
                                                                                     cTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATTTTAAACAGTAGTTa      <  2:3370490/91‑1 (MQ=255)
                                                                                      tCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTAt     <  1:443399/91‑1 (MQ=255)
                                                                                         caAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTg  >  1:2075443/1‑91 (MQ=255)
                                                                                         caAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTg  <  1:826809/91‑1 (MQ=255)
                                                                                          |                                                                                         
ATGAATTACAATATCCTGGTCATATCAAGGTGACAGTTGTTCGAGAGACTAGAGCAGTAGAATATGCGAAATAATTTTTGTCTCCCTCACAAATTAGTGAGGGAGCTTTTTTAAGTTGTAGTCTTAATCTAGTTAGACAGCACTTTATCGGTAATAACTATATTAAACAGTAGTTATTTG  >  USA300TCH1516_ALE/1313023‑1313202

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