Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 109,958 G→A 7.4% intergenic (+37/‑184) 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,9580GA 7.4% 801.1 / 33.8 244intergenic (+37/‑184)plc/USA300HOU_RS005151‑phosphatidylinositol phosphodiesterase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (138/88);  new base A (8/10);  total (146/98)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 252 total aligned reads displayed.

ACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAA  >  USA300TCH1516_ALE/109868‑110048
                                                                                           |                                                                                          
aCGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAg                                                                                            <  2:590314/91‑1 (MQ=255)
 cGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGt                                                                                           >  2:1716265/1‑91 (MQ=255)
  gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                          >  2:3281183/1‑91 (MQ=255)
  gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                          >  1:1222115/1‑91 (MQ=255)
  gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                          >  1:2794383/1‑91 (MQ=255)
  gTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGtt                                                                                          <  1:1580453/91‑1 (MQ=255)
      ccTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGAt                                                                                      >  1:3421090/1‑91 (MQ=255)
      ccTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGAt                                                                                      <  1:1019218/91‑1 (MQ=255)
        tGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTc                                                                                    >  2:2657241/1‑91 (MQ=255)
        tGGATATGATGATATCGTAAGTGAAATCATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTc                                                                                    >  2:813613/1‑91 (MQ=255)
         ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                   >  1:889946/1‑91 (MQ=255)
         ggATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCt                                                                                   <  1:2240531/91‑1 (MQ=255)
          gATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTa                                                                                  >  2:1049376/1‑91 (MQ=255)
          gATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTa                                                                                  >  2:3485449/1‑91 (MQ=255)
           atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                 <  1:197668/91‑1 (MQ=255)
           atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                 <  2:2286862/91‑1 (MQ=255)
           atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                 >  2:1122688/1‑91 (MQ=255)
           atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                 <  1:82214/91‑1 (MQ=255)
           atatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTaa                                                                                 <  1:505676/91‑1 (MQ=255)
            tatGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAat                                                                                <  2:1354152/91‑1 (MQ=255)
             atgatgATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatt                                                                               <  2:3135759/91‑1 (MQ=255)
              tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAatta                                                                              >  2:3459880/1‑91 (MQ=255)
              tgatgaTATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAATTAGAATCTaacta                                                                              <  2:1788702/91‑4 (MQ=39)
               gatgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAattat                                                                             >  1:2251700/1‑91 (MQ=255)
                 tgatATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAATTAGAATCTAACTAttt                                                                           <  1:1741823/91‑1 (MQ=38)
                   atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGCTTCAATAATGATAATAAGACGAGTATGGAAATAGTTAGATTCTAATTATTTTc                                                                         >  1:3611363/1‑91 (MQ=255)
                   atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                         >  1:2185703/1‑91 (MQ=255)
                   atatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTc                                                                         >  1:2185604/1‑91 (MQ=255)
                    tatCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCa                                                                        >  2:1030744/1‑91 (MQ=255)
                     atCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCa                                                                        <  1:563921/91‑1 (MQ=25)
                      tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCat                                                                       <  2:981106/91‑2 (MQ=25)
                      tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                      >  1:1464661/1‑91 (MQ=255)
                      tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCact                                                                      <  2:3323073/91‑1 (MQ=255)
                      tCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAAGAGTTAGATTCTAATTATTTTCact                                                                      >  2:3600026/1‑91 (MQ=255)
                       cGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTGAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCacta                                                                     >  2:3417765/1‑91 (MQ=255)
                        gTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCactac                                                                    <  1:2086622/91‑1 (MQ=255)
                          aaGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTc                                                                  <  1:761959/91‑1 (MQ=255)
                           aGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCg                                                                 >  1:2430708/1‑91 (MQ=255)
                             tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGtc                                                               <  1:2886124/91‑2 (MQ=18)
                             tGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtt                                                               <  2:352106/91‑1 (MQ=255)
                              gAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGttt                                                              <  1:3483792/91‑1 (MQ=255)
                               aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                             <  1:3442096/91‑1 (MQ=255)
                               aaaTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGtttt                                                             >  2:3439882/1‑91 (MQ=255)
                                 aTTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTa                                                           >  2:400597/1‑91 (MQ=255)
                                    atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                        >  1:3687795/1‑91 (MQ=255)
                                    atagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAttt                                                        >  1:3170895/1‑91 (MQ=255)
                                     tagatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTAtttt                                                       <  1:3466911/91‑1 (MQ=255)
                                      agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTa                                                                                         <  2:2897098/56‑1 (MQ=255)
                                      agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTa                                                                                         >  1:2897098/1‑56 (MQ=255)
                                      agatagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTg                                                      >  1:1283020/1‑91 (MQ=255)
                                        atagTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGaa                                                    >  1:1366190/1‑91 (MQ=255)
                                           gTAATAAATAAGGATTCAATAATGATATTAAGACGGGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAAt                                                 >  2:1275124/1‑91 (MQ=255)
                                             aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                               <  2:2346258/91‑1 (MQ=255)
                                             aataaataAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATaa                                               <  1:2067066/91‑1 (MQ=255)
                                              ataaataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAg                                              <  2:2339483/91‑1 (MQ=255)
                                                aaataaGGATTCAATAATGCTATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGta                                            >  2:1228587/1‑91 (MQ=255)
                                                 aataaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaa                                           >  2:3128958/1‑91 (MQ=255)
                                                   taaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaata                                         >  1:2471739/1‑91 (MQ=255)
                                                    aaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaataa                                        <  2:1464661/91‑1 (MQ=255)
                                                    aaGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGtaataa                                        >  2:3069551/1‑91 (MQ=255)
                                                     aGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTaataat                                       >  2:2027631/1‑91 (MQ=255)
                                                      ggATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAAtt                                      <  1:1299028/91‑1 (MQ=255)
                                                         ttCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTcaa                                   >  1:944400/1‑91 (MQ=255)
                                                          tCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTcaac                                  >  2:2933603/1‑91 (MQ=255)
                                                             ataatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                               >  1:2325943/1‑91 (MQ=255)
                                                             ataatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAat                               <  2:2834566/91‑1 (MQ=255)
                                                              taatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAata                              <  1:190063/91‑1 (MQ=255)
                                                               aatGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAatat                             <  1:3246116/91‑1 (MQ=255)
                                                                atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                            >  2:1976804/1‑91 (MQ=255)
                                                                atGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatt                            >  1:2261654/1‑91 (MQ=255)
                                                                 tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                           >  2:2456411/1‑91 (MQ=255)
                                                                 tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                           <  1:10469/91‑1 (MQ=255)
                                                                 tGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAAtatta                           <  1:2091880/91‑1 (MQ=255)
                                                                     atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAATATAAGTAATAATTCAACAATATTATaaa                       >  2:3790383/1‑91 (MQ=255)
                                                                     atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                       >  2:2749666/1‑91 (MQ=255)
                                                                     atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                       >  2:2844137/1‑91 (MQ=255)
                                                                     atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                       >  1:3592554/1‑91 (MQ=255)
                                                                     atTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATaaa                       >  1:2952339/1‑91 (MQ=255)
                                                                       tAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAAtt                     >  2:2078086/1‑91 (MQ=255)
                                                                        aaaaCGAGTAGTGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTg                    >  2:2013453/4‑91 (MQ=17)
                                                                        aaGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTg                    >  2:2134357/1‑91 (MQ=255)
                                                                         aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                   >  1:1872787/1‑91 (MQ=255)
                                                                         aGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACACGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGa                   <  2:147880/91‑1 (MQ=255)
                                                                          gACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGaa                  <  2:565645/91‑1 (MQ=255)
                                                                           aCGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                 <  2:2505083/91‑1 (MQ=255)
                                                                            cGAGTAGTGAAAATAATTAGAATCTAACTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAAc                 >  1:649012/1‑91 (MQ=25)
                                                                            cGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACa                <  1:3750243/91‑1 (MQ=255)
                                                                              aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa              >  1:2657712/1‑91 (MQ=255)
                                                                              aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa              >  1:2907329/1‑91 (MQ=255)
                                                                              aGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGa              <  2:504737/91‑1 (MQ=255)
                                                                               gTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTc                                     >  1:957360/1‑67 (MQ=255)
                                                                                 agtGAAAATAATTAGAATCTAACTATTTTCATACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg           >  2:3789574/3‑91 (MQ=25)
                                                                                 aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg           >  1:2598789/1‑91 (MQ=255)
                                                                                 aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg           >  1:3047/1‑91 (MQ=255)
                                                                                 aTGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTCGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttg           >  1:3078476/1‑91 (MQ=255)
                                                                                    gAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtt         >  1:3796665/1‑91 (MQ=255)
                                                                                     aaaaTAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgttt        <  2:2179898/91‑1 (MQ=255)
                                                                                      aaaTAGTTAGATTCTAGTTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGAttgtttg       <  2:3706291/91‑1 (MQ=255)
                                                                                         tAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGa    >  1:1415994/1‑91 (MQ=255)
                                                                                           gTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGaaa  >  2:3008841/1‑91 (MQ=255)
                                                                                           |                                                                                          
ACGTGGCCTGGATATGATGATATCGTAAGTGAAATTATAGATAGTAATAAATAAGGATTCAATAATGATATTAAGACGAGTATGAAAATAGTTAGATTCTAATTATTTTCACTACTCGTTTTTATTTTGAAAATAAGTAATAATTCAACAATATTATAAATTGAACAGATTGTTTGTGAAA  >  USA300TCH1516_ALE/109868‑110048

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: