Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 77,984 T→C 7.8% intergenic (+93/‑334) USA300HOU_RS00350 → / → USA300HOU_RS00355 Monoacylglycerol lipase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE77,9840TC 7.8% 862.9 / 42.4 289intergenic (+93/‑334)USA300HOU_RS00350/USA300HOU_RS00355Monoacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (130/135);  new base C (17/6);  total (147/142)
Fisher's exact test for biased strand distribution p-value = 2.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 305 total aligned reads displayed.

CAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTACTC  >  USA300TCH1516_ALE/77894‑78074
                                                                                           |                                                                                          
 aGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTT                                                                                           <  2:3240491‑M1/91‑1 (MQ=255)
 aGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTT                                                                                           >  2:3455486‑M1/1‑91 (MQ=255)
  gTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTT                                                                                          >  1:2283166‑M1/1‑91 (MQ=255)
   tataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTC                                                                                         <  1:1145107‑M1/91‑1 (MQ=255)
    ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCA                                                                                        >  1:1252595‑M1/1‑91 (MQ=255)
    ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCA                                                                                        <  1:283636‑M1/91‑1 (MQ=255)
    ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCA                                                                                        <  2:2568794‑M1/91‑1 (MQ=255)
     taGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCAT                                                                                       >  1:853498‑M1/1‑91 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATTGATTTTTTCATAGC                                                                                    <  2:3358707‑M1/91‑1 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGC                                                                                    <  2:1615603‑M1/91‑1 (MQ=255)
        ttCCGTGTATTGGGTTATAAGAAATTATGGGGATATGAAACGTACTAAGCTTAGCTATGAAGAAATCTAGGACGATAGATTTTTTCATAGc                                                                                    <  2:170550/91‑1 (MQ=21)
          ccGTGTATTTGATTCTAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTA                                                                                  >  2:2586526‑M1/1‑91 (MQ=255)
          ccGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTA                                                                                  <  1:1518‑M1/91‑1 (MQ=255)
            gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtt                                                                                <  1:2171757/91‑1 (MQ=255)
             tgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTAatc                                                                               <  2:1180021‑M1/91‑4 (MQ=255)
               tATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttttt                                                                             >  2:307702/1‑91 (MQ=255)
               tATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttttt                                                                             >  2:1841754/1‑91 (MQ=255)
               tATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTActttt                                                                             <  1:3125797/91‑1 (MQ=255)
                 tttGATTATAAGAAATTATGAGGATATTAAATATATTAAGATTAGCTATGAAAAAATCTATCGTTATAGATTTCTTCATAGCTAatcttag                                                                           <  2:1766633‑M1/91‑8 (MQ=255)
                 tttGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATCGTCATAGATTTCTTCATAGCTAatcttag                                                                           >  1:559112‑M1/1‑84 (MQ=255)
                    gATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatag                                                                        <  2:2255472/91‑1 (MQ=255)
                     aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                                       <  2:1083897/91‑1 (MQ=255)
                       tataAGATATTATGAGGATATTAAACATACTATGTTTAGCTATGAAGAAATCTATGACGATAGATTGTTTCATAGCTATTTTttatagtta                                                                     <  2:2172547/91‑1 (MQ=18)
                       tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                                     <  1:340966/91‑1 (MQ=255)
                       tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                                     <  2:801755/91‑1 (MQ=255)
                       tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAGTTTttatagtta                                                                     <  1:3048356/91‑1 (MQ=255)
                        ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                                    >  2:2200912/1‑91 (MQ=255)
                        ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                                    >  2:682227/1‑91 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                                   <  1:2228738/91‑1 (MQ=255)
                          aaGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttatag                                                                  >  2:388206/1‑91 (MQ=255)
                            gAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagag                                                                >  2:437582/1‑91 (MQ=255)
                             aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                               <  2:2897238/91‑1 (MQ=255)
                                  aTGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGCGAGGAGt                                                          >  1:1956885/1‑91 (MQ=255)
                                   tGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTa                                                         <  1:1103667/91‑1 (MQ=255)
                                     aGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGa                                                       >  1:583516/1‑91 (MQ=255)
                                      ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                                      >  2:2784387/1‑91 (MQ=255)
                                        atatTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTg                                                    <  2:3169879/91‑1 (MQ=255)
                                        atatTAAACATACTAAGAGTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTg                                                    <  1:2426637/91‑1 (MQ=255)
                                             aaaCATACTAAGATTAGCTCTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTTCaa                                               >  2:2086803/1‑90 (MQ=255)
                                             aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                               >  1:1563685/1‑91 (MQ=255)
                                             aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                               <  1:2923511/91‑1 (MQ=255)
                                             aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                               <  1:1645639/91‑1 (MQ=255)
                                             aaaCATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAg                                               >  2:3848833/1‑91 (MQ=255)
                                                cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                            >  1:483009/1‑91 (MQ=255)
                                                cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                            >  2:872756/1‑91 (MQ=255)
                                                  tACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActct                                          >  2:2468044/1‑91 (MQ=255)
                                                     tAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTgg                                       <  2:471180/91‑1 (MQ=255)
                                                      aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGCTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                                      >  2:2412628/1‑91 (MQ=255)
                                                       aaaatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                                     >  2:1748504/5‑91 (MQ=17)
                                                       aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGCCTCTTGGAt                                     >  1:3784346/1‑91 (MQ=255)
                                                       aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                                     <  1:428575/91‑1 (MQ=255)
                                                         aatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   >  2:1029627/3‑91 (MQ=17)
                                                         aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   >  2:1924693/1‑91 (MQ=255)
                                                         aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   <  1:1559035/91‑1 (MQ=255)
                                                         aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   <  1:3640573/91‑1 (MQ=255)
                                                          ttAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                                  >  2:2497896/1‑91 (MQ=255)
                                                          ttAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                                  >  2:1511807/1‑91 (MQ=255)
                                                          atAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                                  >  1:2385044/2‑91 (MQ=17)
                                                            aGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTaa                                >  2:2309844/1‑91 (MQ=255)
                                                            aGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTaa                                >  1:3286485/1‑91 (MQ=17)
                                                             gCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTaaa                               >  2:2985784/1‑91 (MQ=255)
                                                             gCTATGAAAAACTCTATCGTCATAGCTTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGCGTAGACTGTCCAGACTCTTGGATTTTcaa                               >  1:2933711/1‑91 (MQ=17)
                                                             gCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTaaa                               >  1:1026248/1‑91 (MQ=17)
                                                              cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                              <  1:220371/91‑1 (MQ=255)
                                                              cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                              <  2:2527232/91‑1 (MQ=255)
                                                              cTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                              >  2:3773655/1‑91 (MQ=17)
                                                                aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                            >  2:1111310/1‑91 (MQ=255)
                                                                 tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                           <  1:3934888/91‑1 (MQ=255)
                                                                  gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                          <  1:802009/91‑1 (MQ=255)
                                                                  gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                          >  1:701206/1‑91 (MQ=255)
                                                                  gaagaaATCTATGACGATAGATTTTGTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                          <  1:3690948/91‑1 (MQ=255)
                                                                    agaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa                        >  2:9173/1‑91 (MQ=255)
                                                                    agaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa                        <  1:2717772/91‑1 (MQ=255)
                                                                    agaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTACAGACTCTTGGATTTTAAATCCGTaa                        <  2:3594398/91‑1 (MQ=255)
                                                                     gaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaa                       >  1:1843114/1‑91 (MQ=255)
                                                                       aaTCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGACAGGAGTAGACTGTCCAGACACTTGGATTTTAAATTTGTaaaaa                     >  2:3684438/1‑91 (MQ=17)
                                                                        aTCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaa                    >  1:598805/1‑91 (MQ=255)
                                                                        aTCGATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAACCGTaaaaaa                    <  2:2025845/91‑1 (MQ=255)
                                                                         tCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaa                   >  2:2852213/1‑91 (MQ=255)
                                                                           tATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAAc                 <  1:2880051/91‑1 (MQ=255)
                                                                            aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa                >  1:417249/1‑91 (MQ=255)
                                                                            aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa                >  1:1716293/1‑91 (MQ=255)
                                                                             aggACGATAGGTTTTTTCATAGCTATTTTTTATAGTTGGAGAGAGGAGTAGACTGTCCTGACTCTTGGATTTTAAATACGAAAAAAAAACa                <  2:2657547/89‑1 (MQ=14)
                                                                             aggACGAGAGATTGTGTCAAAGCGAGTTTTTATAGTTAGAGAGAGGAGTAGACTGTCCAGACGCTTGGATTTTAAATCAGAAAAAAAAACa                <  2:3260881/89‑1 (MQ=255)
                                                                               gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg              <  1:2624129/91‑1 (MQ=255)
                                                                                aCGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGt             >  2:1816430/1‑91 (MQ=255)
                                                                                aCGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGt             >  1:2187334/1‑91 (MQ=255)
                                                                                 cGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCATACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTc            >  2:348172/1‑91 (MQ=255)
                                                                                 cGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTc            >  2:1474814/1‑91 (MQ=255)
                                                                                 cGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTc            <  1:3007939/91‑1 (MQ=255)
                                                                                 cGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTc            <  1:3507670/91‑1 (MQ=255)
                                                                                  gATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCa           >  1:1199956/1‑91 (MQ=255)
                                                                                   aTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAg          <  2:3277159/91‑1 (MQ=255)
                                                                                   aTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAg          <  1:1394314/91‑1 (MQ=255)
                                                                                       aTTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCttt      <  2:563631/91‑1 (MQ=255)
                                                                                        ttttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTa     >  2:1406351/1‑91 (MQ=255)
                                                                                        ttttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTa     >  2:858795/1‑91 (MQ=255)
                                                                                         tttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCATAAAAAAAACAAGTCAGCTTTAc    <  1:2093346/91‑1 (MQ=255)
                                                                                           tttCATATCTTTTTTTTATCGTTATAGAGAGGAGTAGCCTGTCCTGACTGTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTActt  >  2:178757/1‑90 (MQ=255)
                                                                                           tttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTActc  >  1:2256908/1‑91 (MQ=255)
                                                                                           |                                                                                          
CAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTACTC  >  USA300TCH1516_ALE/77894‑78074

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