Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 690,098 A→T 15.3% intergenic (‑27/+8) 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,0980AT15.3% 348.7 / 49.2 123intergenic (‑27/+8)USA300HOU_RS03320/USA300HOU_RS03325hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (51/53);  new base T (5/14);  total (56/67)
Fisher's exact test for biased strand distribution p-value = 8.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 126 total aligned reads displayed.

CTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAA  >  USA300TCH1516_ALE/690008‑690187
                                                                                          |                                                                                         
cTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCa                                                                                           >  1:2463333/1‑91 (MQ=255)
cTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCa                                                                                           <  1:1492249/91‑1 (MQ=255)
 ttCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCaa                                                                                          <  1:1256086/91‑1 (MQ=255)
  tCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAAt                                                                                         >  1:1878806/1‑91 (MQ=255)
   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:3713049/91‑2 (MQ=255)
        gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                   <  1:3713070/91‑2 (MQ=255)
        gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                   <  1:2226215/91‑2 (MQ=255)
        gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggt                                                                                   <  1:2225884/91‑2 (MQ=255)
        gTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGc                                                                                   <  1:17967/91‑1 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                  <  1:800368/91‑3 (MQ=255)
         ttttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtt                                                                                  <  1:800152/91‑3 (MQ=255)
          tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                 <  1:2340500/91‑4 (MQ=255)
          tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTggggtta                                                                                 <  1:731506/91‑4 (MQ=255)
          tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTa                                                                                 >  1:2216619/1‑91 (MQ=255)
          tttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTa                                                                                 >  1:3225963/1‑91 (MQ=255)
           ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                <  1:1189600/91‑1 (MQ=255)
           ttCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAt                                                                                <  1:1102125/91‑1 (MQ=255)
            tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAGCCCCTATTGGGGTTAtt                                                                               <  1:3340876/91‑1 (MQ=255)
            tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                               <  1:2313710/91‑1 (MQ=255)
            tCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAtt                                                                               <  1:2313393/91‑1 (MQ=255)
             cGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTAttt                                                                              >  1:837680/1‑91 (MQ=255)
              gTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGCCTATTTa                                                                             >  1:2085402/1‑91 (MQ=255)
                 aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                          <  1:2098044/91‑1 (MQ=255)
                 aataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAAt                                                                          >  1:3496904/1‑91 (MQ=255)
                  ataaCAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATg                                                                         <  1:3972797/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)
                       aaaTTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAtt                                                                    <  1:2426900/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:1368435/1‑91 (MQ=255)
                            aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                               >  1:2635182/1‑91 (MQ=255)
                            aTCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTAttttttt                                                               >  1:2634939/1‑91 (MQ=255)
                                   aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                        <  1:1463278/91‑1 (MQ=255)
                                   aaaCCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGa                                                        <  1:560177/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:537625/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                  >  1:878465/1‑91 (MQ=255)
                                         gCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCt                                                  >  1:197550/1‑91 (MQ=255)
                                          cATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTa                                                 <  1:1278711/91‑1 (MQ=255)
                                            tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                               >  1:363728/1‑91 (MQ=255)
                                            tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                               >  1:362329/1‑91 (MQ=255)
                                            tACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTAct                                               >  1:362317/1‑91 (MQ=255)
                                                    ttttAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTAATTTTTTTATCATGATATTCTACTTCTtcatc                                       >  1:3528615/1‑91 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                  <  1:3637444/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                  <  1:3897859/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                  <  1:74485/91‑1 (MQ=255)
                                                         aaTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAAc                                  >  1:2770504/1‑91 (MQ=255)
                                                          aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                 >  1:3740872/1‑91 (MQ=255)
                                                          aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                 >  1:2482337/1‑91 (MQ=255)
                                                          aTTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACt                                 <  1:2343374/91‑1 (MQ=255)
                                                            tCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTa                               >  1:518776/1‑91 (MQ=255)
                                                            tCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTa                               >  1:519310/1‑91 (MQ=255)
                                                              aTGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAt                             <  1:2407799/91‑1 (MQ=255)
                                                               tGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAAtt                            >  1:1368084/1‑91 (MQ=255)
                                                                 tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                          >  1:1674958/1‑91 (MQ=255)
                                                                 tCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTa                          <  1:4080863/91‑1 (MQ=255)
                                                                   aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                        >  1:2811094/1‑91 (MQ=255)
                                                                   aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                        <  1:3265505/91‑1 (MQ=255)
                                                                   aTCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAAc                        <  1:946273/91‑1 (MQ=255)
                                                                    tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                       >  1:1826875/1‑91 (MQ=255)
                                                                    tCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACt                       >  1:1875606/1‑91 (MQ=255)
                                                                       cTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACtttt                    <  1:2944272/91‑1 (MQ=255)
                                                                        tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                   <  1:1602656/91‑1 (MQ=255)
                                                                        tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                   <  1:3083575/91‑1 (MQ=255)
                                                                        tCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTc                   <  1:1348628/91‑1 (MQ=255)
                                                                         ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                  <  1:239018/91‑1 (MQ=255)
                                                                         ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                  <  1:238979/91‑1 (MQ=255)
                                                                         ccAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCt                  <  1:239098/91‑1 (MQ=255)
                                                                          cAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCtt                 <  1:2119100/91‑1 (MQ=255)
                                                                            atttaaaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTaa               >  1:1245626/5‑91 (MQ=37)
                                                                             tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATGATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc              <  1:3287492/91‑1 (MQ=255)
                                                                             tAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAAc              <  1:3559417/91‑1 (MQ=255)
                                                                              aaaaaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACt             <  1:2595097/91‑1 (MQ=255)
                                                                                 aaTAGCCCCTATTGGGGTTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca          >  1:3110697/1‑91 (MQ=255)
                                                                                 aaTAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGca          <  1:3501840/91‑1 (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:1287514/91‑1 (MQ=255)
                                                                                     aCCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACAtt      <  1:1287547/91‑1 (MQ=255)
                                                                                      ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa     >  1:3781043/1‑91 (MQ=255)
                                                                                      ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa     >  1:1365776/1‑91 (MQ=255)
                                                                                      ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa     >  1:1078421/1‑91 (MQ=255)
                                                                                      ccccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTa     <  1:2599744/91‑1 (MQ=255)
                                                                                       cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa    >  1:3207891/1‑91 (MQ=255)
                                                                                       cccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaa    >  1:3208031/1‑91 (MQ=255)
                                                                                        ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa   >  1:2825228/1‑91 (MQ=255)
                                                                                        ccAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaa   <  1:82643/91‑1 (MQ=255)
                                                                                         cAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTaaaa  >  1:1695619/1‑91 (MQ=255)
                                                                                         cAATAGGGGCTATTTAAATGTTATTTTTTTAGCATGATATTCTAGTTCTTCATTATAGCTTAATTTAACTTTTCTTAACTGCACATTaaaa  <  1:3311904/91‑1 (MQ=255)
                                                                                          |                                                                                         
CTTCGCCAGTTTTCGTTAATAACAAATTATCTCTTAAACCTGCATACGGTGCTTTTAAATTCATGTCATCCCTCCAATAAAAATAACCCCAATAGGGGCTATTTAAATGTTATTTTTTTATCATGATATTCTACTTCTTCATCATAACTTAATTTAACTTTTCTTAACTGCACATTAAAA  >  USA300TCH1516_ALE/690008‑690187

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: