Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 690,101 A→T 15.8% intergenic (‑30/+5) USA300HOU_RS03320 ← / ← USA300HOU_RS03325 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE690,1010AT15.8% 329.0 / 51.9 120intergenic (‑30/+5)USA300HOU_RS03320/USA300HOU_RS03325hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (48/53);  new base T (5/14);  total (53/67)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 123 total aligned reads displayed.

CGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAA  >  USA300TCH1516_ALE/690011‑690191
                                                                                          |                                                                                          
cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                            >  1:4066488/1‑91 (MQ=255)
cGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATa                                                                                            >  1:4065887/1‑91 (MQ=255)
 gCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAg                                                                                           >  1:4004640/1‑91 (MQ=255)
   cAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAggg                                                                                         <  1:1523916/91‑1 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                       <  1:2225884/91‑2 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                       <  1:2226215/91‑2 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                       <  1:3713049/91‑2 (MQ=255)
     gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                       <  1:860001/91‑2 (MQ=255)
      ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                      <  1:800152/91‑3 (MQ=255)
      ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                      <  1:800368/91‑3 (MQ=255)
       tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                     <  1:731506/91‑4 (MQ=255)
       tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                     <  1:423671/91‑4 (MQ=255)
       tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                     <  1:2340500/91‑4 (MQ=255)
       tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGTGCCTa                                                                                     >  1:1814818/1‑91 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAACATAACCCCAATAGGGGCTAt                                                                                    <  1:3334394/91‑1 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                    <  1:1189600/91‑1 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                    <  1:1102125/91‑1 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                    <  1:3438701/91‑1 (MQ=255)
        ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAt                                                                                    <  1:2341015/91‑1 (MQ=255)
         tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                                   <  1:2313393/91‑1 (MQ=255)
         tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                                   <  1:3651481/91‑1 (MQ=255)
           gTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGCCTATTTa                                                                                 >  1:2085402/1‑91 (MQ=255)
              aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                              >  1:3496904/1‑91 (MQ=255)
              aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                              <  1:2098044/91‑1 (MQ=255)
                 aaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGtt                                                                           >  1:679291/1‑91 (MQ=255)
                   cAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAt                                                                         <  1:3429254/91‑1 (MQ=255)
                    aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                        <  1:2426814/91‑1 (MQ=255)
                     aaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttt                                                                       <  1:2095022/91‑1 (MQ=255)
                       ttATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttt                                                                     <  1:2339649/91‑1 (MQ=255)
                        tATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtttttt                                                                    >  1:432966/1‑91 (MQ=255)
                         aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                   >  1:2635182/1‑91 (MQ=255)
                         aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                   >  1:2635155/1‑91 (MQ=255)
                         aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                                   >  1:1368435/1‑91 (MQ=255)
                           ctctTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTAt                                                                 >  1:4011673/1‑91 (MQ=255)
                                aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                            <  1:560177/91‑1 (MQ=255)
                                aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                            <  1:1463278/91‑1 (MQ=255)
                                 aaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGat                                                           <  1:267220/91‑1 (MQ=255)
                                 aaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGat                                                           <  1:267071/91‑1 (MQ=255)
                                   ccTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGatat                                                         >  1:4135983/1‑91 (MQ=255)
                                    cTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATAtt                                                        >  1:537640/1‑91 (MQ=255)
                                    cTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATAtt                                                        >  1:537625/1‑91 (MQ=255)
                                      gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                      >  1:878465/1‑91 (MQ=255)
                                      gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                      >  1:197550/1‑91 (MQ=255)
                                       cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                     <  1:3768329/91‑1 (MQ=255)
                                         tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                   >  1:363728/1‑91 (MQ=255)
                                         tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                                   >  1:362317/1‑91 (MQ=255)
                                           cGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTActtc                                                 >  1:211427/1‑91 (MQ=255)
                                                 ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTAATTTTTTTATCATGATATTCTACTTCTtcatc                                           >  1:3528615/1‑91 (MQ=255)
                                                    tAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATa                                        >  1:2973441/1‑91 (MQ=255)
                                                      aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                      <  1:3897859/91‑1 (MQ=255)
                                                      aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                      >  1:2770504/1‑91 (MQ=255)
                                                      aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                      <  1:3637444/91‑1 (MQ=255)
                                                       aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                     >  1:2482337/1‑91 (MQ=255)
                                                       aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                     >  1:3740872/1‑91 (MQ=255)
                                                       aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                     <  1:2343374/91‑1 (MQ=255)
                                                         tCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTa                                   >  1:518776/1‑91 (MQ=255)
                                                          cATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTaa                                  >  1:1424503/1‑91 (MQ=255)
                                                           aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                 <  1:2407799/91‑1 (MQ=255)
                                                           aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                                 <  1:778392/91‑1 (MQ=255)
                                                            tGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAtt                                >  1:1368084/1‑91 (MQ=255)
                                                              tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                              >  1:1674958/1‑91 (MQ=255)
                                                                aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                            <  1:946273/91‑1 (MQ=255)
                                                                aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                            >  1:2811094/1‑91 (MQ=255)
                                                                aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                            <  1:3265505/91‑1 (MQ=255)
                                                                 tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                           >  1:1875606/1‑91 (MQ=255)
                                                                 tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                           >  1:1826875/1‑91 (MQ=255)
                                                                    cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                        <  1:2944272/91‑1 (MQ=255)
                                                                     tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                       <  1:1348628/91‑1 (MQ=255)
                                                                      ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                      <  1:1638754/91‑1 (MQ=255)
                                                                      ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                      <  1:239018/91‑1 (MQ=255)
                                                                      ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                      <  1:239131/91‑1 (MQ=255)
                                                                      ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                      <  1:239098/91‑1 (MQ=255)
                                                                       cAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCtt                     <  1:2119100/91‑1 (MQ=255)
                                                                         atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                   >  1:1245771/5‑91 (MQ=37)
                                                                         atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa                   >  1:1245626/5‑91 (MQ=37)
                                                                         aTACAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTCTCTTaa                   <  1:3220580/91‑1 (MQ=255)
                                                                          tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATGATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc                  <  1:3287492/91‑1 (MQ=255)
                                                                          tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc                  <  1:3559417/91‑1 (MQ=255)
                                                                           aaaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACt                 <  1:2595097/91‑1 (MQ=255)
                                                                             aaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGc               >  1:2260855/1‑91 (MQ=255)
                                                                              aaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca              >  1:3110697/1‑91 (MQ=255)
                                                                               aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTGATCATAACTTAATTTAACTTTTCTTAACTGcac             >  1:3404598/1‑91 (MQ=255)
                                                                               aTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGcac             >  1:3404578/1‑91 (MQ=255)
                                                                                 agccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAt           >  1:933621/3‑91 (MQ=255)
                                                                                  gccccTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          >  1:1675409/2‑91 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          >  1:356908/1‑91 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:1287514/91‑1 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:2465576/91‑1 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:1289164/91‑1 (MQ=255)
                                                                                  aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt          <  1:1287547/91‑1 (MQ=255)
                                                                                   ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa         >  1:1078421/1‑91 (MQ=255)
                                                                                   ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa         <  1:2599744/91‑1 (MQ=255)
                                                                                   ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa         >  1:3781043/1‑91 (MQ=255)
                                                                                    cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa        >  1:3208031/1‑91 (MQ=255)
                                                                                    cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa        >  1:3207891/1‑91 (MQ=255)
                                                                                     ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa       <  1:82643/91‑1 (MQ=255)
                                                                                     ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa       >  1:2825228/1‑91 (MQ=255)
                                                                                      cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa      >  1:1695619/1‑91 (MQ=255)
                                                                                      cAATAGGGGCTATTTAAATGTTATTTTTTTAGCATGATATTCTAGTTCTTCATTATAGCTTAATTTAACTTTTCTTAACTGCACATTaaaa      <  1:3311904/91‑1 (MQ=255)
                                                                                          aGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACaa  >  1:3943371/1‑91 (MQ=255)
                                                                                          |                                                                                          
CGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAAACAA  >  USA300TCH1516_ALE/690011‑690191

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