Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 109,949:1 +G 5.2% intergenic (+28/‑193) plc → / → USA300HOU_RS00515 1‑phosphatidylinositol phosphodiesterase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE109,9491.G 5.2% 695.7 / 28.6 155intergenic (+28/‑193)plc/USA300HOU_RS005151‑phosphatidylinositol phosphodiesterase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (101/46);  new base G (3/5);  total (104/51)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 163 total aligned reads displayed.

GCAGGATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTG  >  USA300TCH1516_ALE/109859‑110039
                                                                                           |                                                                                          
gCAGGATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTa                                                                                             >  1:2467759/1‑91 (MQ=255)
gCAGGATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTa                                                                                             <  1:372212/91‑1 (MQ=255)
 cAGGATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAt                                                                                           >  1:196719/1‑91 (MQ=255)
   ggATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGa                                                                                         >  1:2960060/1‑91 (MQ=255)
   ggATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGa                                                                                         >  1:2960397/1‑91 (MQ=255)
      tataCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGaaaa                                                                                      >  1:2673902/1‑91 (MQ=255)
      tataCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGaaaa                                                                                      >  1:2728241/1‑91 (MQ=255)
      tataCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGaaaa                                                                                      >  1:53539/1‑91 (MQ=255)
       ataCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAAt                                                                                     <  1:2353993/91‑1 (MQ=255)
        taCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATa                                                                                    <  1:2289914/91‑1 (MQ=255)
          cGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGt                                                                                  >  1:3080960/1‑91 (MQ=255)
          cGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGt                                                                                  >  1:4021603/1‑91 (MQ=255)
          cGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGt                                                                                  >  1:3882024/1‑91 (MQ=255)
           gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                 <  1:1814171/91‑1 (MQ=255)
           gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                 <  1:1813975/91‑1 (MQ=255)
            tGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTa                                                                                >  1:1825370/1‑91 (MQ=255)
             ggCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGCGTATGAAAATAGTTAg                                                                               >  1:1484914/1‑91 (MQ=255)
             ggCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAg                                                                               >  1:3377201/1‑91 (MQ=255)
             ggCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAg                                                                               >  1:1977571/1‑91 (MQ=255)
               ccTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGAt                                                                             >  1:466731/1‑91 (MQ=255)
                 tGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTc                                                                           >  1:419902/1‑91 (MQ=255)
                   aatatGATGATATCGTAAGAGAAATTATAGATAGTAATAAAGAAGGATTCAATAATGATATTTAGACGAGTATGAAAATAGTGAGATTCTa                                                                         <  1:704669/90‑1 (MQ=255)
                    atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                        >  1:786243/1‑91 (MQ=255)
                    atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                        >  1:646007/1‑91 (MQ=255)
                    atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                        <  1:376933/91‑1 (MQ=255)
                     tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                       >  1:2382464/1‑91 (MQ=255)
                     tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                       >  1:2383496/1‑91 (MQ=255)
                       tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                     >  1:329016/1‑91 (MQ=255)
                       tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                     >  1:330357/1‑91 (MQ=255)
                         atgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAattatt                                                                   >  1:337989/1‑91 (MQ=255)
                         atgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAattatt                                                                   <  1:1658171/91‑1 (MQ=255)
                          tgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTAttt                                                                  >  1:1140706/1‑91 (MQ=255)
                            atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                >  1:1002141/1‑91 (MQ=255)
                             tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTc                                                                <  1:3071598/91‑1 (MQ=25)
                             tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTc                                                                <  1:1010776/91‑1 (MQ=25)
                               tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                             >  1:2732255/1‑91 (MQ=255)
                                 gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCactac                                                           >  1:367500/1‑91 (MQ=255)
                                  tAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCatact                                                           <  1:4082427/91‑5 (MQ=21)
                                  tAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCatact                                                           <  1:3559161/91‑5 (MQ=21)
                                   aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTCAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTc                                                         >  1:3291041/1‑91 (MQ=255)
                                   aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTc                                                         <  1:696217/91‑1 (MQ=21)
                                   aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTc                                                         <  1:3630619/91‑1 (MQ=255)
                                    aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                        >  1:3835912/1‑91 (MQ=255)
                                    aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                        >  1:3563735/1‑91 (MQ=255)
                                     gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                       <  1:3955203/91‑1 (MQ=255)
                                     gTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGt                                                       <  1:3955932/91‑1 (MQ=255)
                                      tGAAATTATAGATAGTAATAAATAAGGATTCAGTAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTGTTTTCACTACTCGtt                                                      <  1:3538768/91‑1 (MQ=255)
                                      tGAAATTATAGATAGTAATAAATAAGGATTCAGTAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTGTTTTCACTACTCGtt                                                      <  1:3538023/91‑1 (MQ=255)
                                      tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                      >  1:3073510/1‑91 (MQ=255)
                                       gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                     >  1:1230103/1‑91 (MQ=255)
                                       gAAATTATAGATAGGAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                     <  1:840379/91‑1 (MQ=255)
                                        aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                    <  1:2990780/91‑1 (MQ=255)
                                        aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                    >  1:3752511/1‑91 (MQ=255)
                                         aaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttttt                                                   >  1:870984/1‑91 (MQ=255)
                                          aTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTa                                                  >  1:3915478/1‑91 (MQ=255)
                                             atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                               >  1:1415545/1‑91 (MQ=255)
                                             atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                               >  1:2832641/1‑91 (MQ=255)
                                             atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                               >  1:1415356/1‑91 (MQ=255)
                                               agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTg                                             >  1:4118879/1‑91 (MQ=255)
                                               agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTg                                             >  1:149341/1‑91 (MQ=255)
                                                gatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGa                                            >  1:3485677/1‑91 (MQ=255)
                                                gatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGa                                            >  1:63267/1‑91 (MQ=255)
                                                  tagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaaa                                          >  1:287191/1‑91 (MQ=255)
                                                      aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                      <  1:1751836/91‑1 (MQ=255)
                                                      aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                      <  1:36723/91‑1 (MQ=255)
                                                      aataaataAGGAGTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                      >  1:3193895/1‑91 (MQ=255)
                                                        taaataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGt                                    >  1:605947/1‑91 (MQ=255)
                                                        taaataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGt                                    >  1:605821/1‑91 (MQ=255)
                                                          aataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaa                                  >  1:2857133/1‑91 (MQ=255)
                                                            taaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaata                                <  1:502889/91‑1 (MQ=255)
                                                                gATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                            >  1:1646860/1‑91 (MQ=255)
                                                                 aTTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTca                           >  1:930397/1‑91 (MQ=255)
                                                                    cAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATACTTcaaca                        >  1:2073293/1‑91 (MQ=255)
                                                                      ataatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                      >  1:1540504/1‑91 (MQ=255)
                                                                       taatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAata                     >  1:730621/1‑91 (MQ=255)
                                                                         atGATATTAAGACGAGTATGAAAATAGTTAGATTCTACTTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                   >  1:3509405/1‑91 (MQ=255)
                                                                         atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                   >  1:1249339/1‑91 (MQ=255)
                                                                          tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                  >  1:1897144/1‑91 (MQ=255)
                                                                           gATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtattat                 >  1:2492306/1‑91 (MQ=255)
                                                                           gATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtattat                 <  1:60882/91‑1 (MQ=255)
                                                                             tatTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaa               >  1:790959/1‑91 (MQ=255)
                                                                              atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa              <  1:2627001/91‑1 (MQ=255)
                                                                              atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa              >  1:2428699/1‑91 (MQ=255)
                                                                              atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa              >  1:2136340/1‑91 (MQ=255)
                                                                                 aaGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTg           >  1:1588404/1‑91 (MQ=255)
                                                                                  aGACGAGTATGAAAATAGTTAGGTTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa          <  1:80979/91‑1 (MQ=255)
                                                                                   gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa         >  1:2806287/1‑91 (MQ=255)
                                                                                   gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa         >  1:2371560/1‑91 (MQ=255)
                                                                                   gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa         <  1:3300696/91‑1 (MQ=255)
                                                                                   gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa         >  1:1847722/1‑91 (MQ=255)
                                                                                   aaCGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa         >  1:3818921/2‑91 (MQ=18)
                                                                                     cGAGTAGTGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc        >  1:1168100/1‑91 (MQ=25)
                                                                                     cGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACa       <  1:2332434/91‑1 (MQ=255)
                                                                                      gAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg      <  1:2184619/91‑1 (MQ=255)
                                                                                      gAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAg      <  1:3621462/91‑1 (MQ=255)
                                                                                       aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa     >  1:2885651/1‑91 (MQ=255)
                                                                                        gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAt    >  1:1653732/1‑91 (MQ=255)
                                                                                         tATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAtt   >  1:2881991/1‑91 (MQ=255)
                                                                                         tATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAtt   >  1:415776/1‑91 (MQ=255)
                                                                                          aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg  <  1:3681186/91‑1 (MQ=255)
                                                                                           |                                                                                          
GCAGGATATACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTG  >  USA300TCH1516_ALE/109859‑110039

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