Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 77,984 T→C 8.7% 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 8.7% 666.6 / 39.9 219intergenic (+93/‑334)USA300HOU_RS00350/USA300HOU_RS00355Monoacylglycerol lipase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (102/98);  new base C (7/12);  total (109/110)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 243 total aligned reads displayed.

CAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTACTC  >  USA300TCH1516_ALE/77894‑78074
                                                                                           |                                                                                          
cAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttt                                                                                            >  2:2805059/1‑91 (MQ=255)
cAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGAtttt                                                                                            >  1:2108949/1‑91 (MQ=255)
    ataGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCa                                                                                        >  2:32792/1‑91 (MQ=255)
      aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATa                                                                                      >  2:921458/1‑91 (MQ=255)
      aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATa                                                                                      <  1:3359322/91‑1 (MQ=255)
      aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATa                                                                                      <  2:3406879/91‑1 (MQ=255)
      aGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATa                                                                                      <  1:3036191/91‑1 (MQ=21)
      aGTACCGTGTATTTGATTAGAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATTACGATAGATTTTTTCATa                                                                                      <  2:3720479/91‑1 (MQ=255)
       gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                                     <  1:3379085/91‑1 (MQ=21)
       gTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAg                                                                                     <  2:2448723/91‑1 (MQ=21)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                                    <  1:1752749/91‑1 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                                    <  1:2965726/91‑1 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                                    <  1:1835152/91‑1 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                                    <  2:3788390/91‑1 (MQ=255)
        ttCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATAATAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGc                                                                                    <  2:3484656/91‑1 (MQ=255)
         tCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCt                                                                                   >  1:3312864/1‑91 (MQ=255)
          ccGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTa                                                                                  >  2:1390731/1‑91 (MQ=255)
           cGTGTATTTGATTATAAGACATTATGAGGATATTAAAGATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAt                                                                                 >  2:3520393/1‑91 (MQ=255)
            gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAg                                                                                     <  2:439466/86‑1 (MQ=255)
            gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAg                                                                                     >  1:439466/1‑86 (MQ=255)
            gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAtt                                                                                >  1:238495/1‑91 (MQ=255)
            gtgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaat                                                                                <  2:2819344/91‑3 (MQ=18)
             tgtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttt                                                                               >  1:765055/1‑91 (MQ=255)
              gtATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTaatct                                                                              <  1:15982/91‑5 (MQ=17)
               tATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTAttttt                                                                             <  1:3686622/91‑1 (MQ=255)
                 tttGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTtta                                                                           >  2:3588342/1‑91 (MQ=255)
                   tGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatt                                                                         >  2:426604/1‑90 (MQ=255)
                     aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                                       <  1:3757961/91‑1 (MQ=255)
                     aTTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagt                                                                       <  2:1149553/91‑1 (MQ=255)
                       tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                                     <  2:2727877/91‑1 (MQ=255)
                       tataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagtta                                                                     <  2:100541/91‑1 (MQ=255)
                        ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                                    <  2:23568/91‑1 (MQ=255)
                        ataAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttat                                                                    <  2:3151965/91‑1 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCAt                                                                                       <  1:1190687/71‑1 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                                   <  2:1953272/91‑1 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                                   <  2:1832985/91‑1 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                                   <  1:2427367/91‑1 (MQ=255)
                         taAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTttatagttata                                                                   >  1:1753012/1‑91 (MQ=255)
                           aGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATaga                                                                 <  2:3206048/91‑1 (MQ=255)
                             aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                               >  2:1955472/1‑91 (MQ=255)
                             aaaTTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATagaga                                                               >  1:3510646/1‑91 (MQ=255)
                                ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGTCGGTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAgagga                                                            <  2:2518665/91‑1 (MQ=255)
                                ttATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATTGCTATTTTTTATAGTTATAGAgagga                                                            <  2:2217044/91‑1 (MQ=255)
                                    gAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAg                                                        >  1:2292555/1‑91 (MQ=255)
                                     aGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTTATAGCTATTTTTTATAGTTATAGAGAGGAGTAGa                                                       >  2:2678982/1‑91 (MQ=255)
                                      ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                                      >  2:120217/1‑91 (MQ=255)
                                      ggATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGAc                                                      >  2:833680/1‑91 (MQ=255)
                                         tatTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGt                                                   <  1:3476304/91‑1 (MQ=255)
                                          atTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTc                                                  <  2:751199/91‑1 (MQ=255)
                                           ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                                 >  1:624948/1‑91 (MQ=255)
                                           ttAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTcc                                                 >  2:1291217/1‑91 (MQ=255)
                                                cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                            >  1:2607304/1‑91 (MQ=255)
                                                cATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGAct                                            >  2:1856627/1‑91 (MQ=255)
                                                 aTACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActc                                           <  1:3297641/91‑1 (MQ=255)
                                                 aTACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActc                                           <  2:1972868/91‑1 (MQ=255)
                                                  tACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGActct                                          >  1:2623302/1‑91 (MQ=255)
                                                    cTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTg                                        <  1:1427435/91‑1 (MQ=255)
                                                     tAAGAGTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTgg                                       <  2:3569031/91‑1 (MQ=255)
                                                      aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                                      >  2:1367004/1‑91 (MQ=255)
                                                      aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                                      >  1:1158487/1‑91 (MQ=255)
                                                      aaGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGa                                      >  1:3085584/1‑91 (MQ=255)
                                                       aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                                     <  2:2138056/91‑1 (MQ=255)
                                                       aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                                     <  1:2126401/91‑1 (MQ=255)
                                                       aGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAt                                     <  2:3543353/91‑1 (MQ=255)
                                                        aaatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtt                                    >  2:563238/4‑91 (MQ=17)
                                                         aatAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   >  2:2673880/3‑91 (MQ=17)
                                                         aTTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAttt                                   <  1:2812970/91‑1 (MQ=255)
                                                          ttAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                                  <  2:2300456/91‑1 (MQ=255)
                                                          atAGCTATGAAAAAATCTATCGTCATAGATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGAtttt                                  >  1:536936/2‑91 (MQ=17)
                                                              cTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAAt                              >  1:3652431/1‑91 (MQ=255)
                                                               tATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATc                             <  1:1072397/91‑1 (MQ=255)
                                                                aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                            >  2:903309/1‑91 (MQ=255)
                                                                aTGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATcc                            >  1:264917/1‑91 (MQ=255)
                                                                 tGAAGAAATCTATGACGNTAGATNTTTTCATAGCTATTTTTTANAGTTATAGAGAGGAGNAGACTGTCCAGACTCTTGGATTTTAAATNCg                           <  2:1020852/91‑1 (MQ=17)
                                                                 tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                           <  1:838269/91‑1 (MQ=255)
                                                                 tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                           <  1:2493558/91‑1 (MQ=255)
                                                                 tGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCg                           <  1:838508/91‑1 (MQ=255)
                                                                  gaagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGt                          >  1:100186/1‑91 (MQ=255)
                                                                   aagaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATACGTa                         <  2:1099938/91‑1 (MQ=255)
                                                                    agaaATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa                        <  2:3629105/91‑1 (MQ=255)
                                                                    agaaATCTATGACGATAGATTTTCTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaa                        <  2:3613280/91‑1 (MQ=255)
                                                                         tCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaa                   >  1:3484791/1‑91 (MQ=255)
                                                                          cTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaaa                  >  1:3397792/1‑91 (MQ=255)
                                                                           tATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTaaaaaaaaa                 <  1:932701/91‑2 (MQ=255)
                                                                            aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa                >  2:3050412/1‑91 (MQ=255)
                                                                            aTGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACa                >  2:3881373/1‑91 (MQ=255)
                                                                             tGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACaa               >  1:3729901/1‑91 (MQ=255)
                                                                               gACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAg              <  2:1644069/91‑1 (MQ=255)
                                                                               gACAATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAATAAACAAg              >  2:3441415/1‑91 (MQ=255)
                                                                                   aTAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAg          >  1:1102760/1‑91 (MQ=255)
                                                                                      gATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCtt       <  1:1869136/91‑1 (MQ=255)
                                                                                      gATTTCTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCtt       >  1:2281050/1‑91 (MQ=255)
                                                                                         tttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTAc    <  1:2285997/91‑1 (MQ=255)
                                                                                          ttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGTCTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTAct   <  2:2996764/91‑1 (MQ=255)
                                                                                          ttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTAct   <  1:1371029/91‑1 (MQ=255)
                                                                                          ttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTAct   >  1:2300844/1‑91 (MQ=255)
                                                                                          ttttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTAct   >  2:1814459/1‑91 (MQ=255)
                                                                                           tttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTActc  >  2:52751/1‑91 (MQ=255)
                                                                                           tttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTActc  <  2:3320425/91‑1 (MQ=255)
                                                                                           tttCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGCCCAGACTCTTGGATTTTAAATCCGTACAAAAAACAAGTCAGCTTTActc  >  2:2944941/1‑91 (MQ=255)
                                                                                           |                                                                                          
CAGTATAGTTCCGTGTATTTGATTATAAGAAATTATGAGGATATTAAACATACTAAGATTAGCTATGAAGAAATCTATGACGATAGATTTTTTCATAGCTATTTTTTATAGTTATAGAGAGGAGTAGACTGTCCAGACTCTTGGATTTTAAATCCGTAAAAAAAACAAGTCAGCTTTACTC  >  USA300TCH1516_ALE/77894‑78074

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