Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,123,779 T→A 5.8% intergenic (‑39/‑94) USA300HOU_RS10825 ← / → lexA_2 hypothetical protein/LexA repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,123,7790TA 5.8% 911.0 / 29.6 269intergenic (‑39/‑94)USA300HOU_RS10825/lexA_2hypothetical protein/LexA repressor
Reads supporting (aligned to +/- strand):  ref base T (131/122);  new base A (8/8);  total (139/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 287 total aligned reads displayed.

ACTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123689‑2123869
                                                                                          |                                                                                          
aCTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAt                                                                                            <  2:3285604/91‑1 (MQ=255)
aCTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAt                                                                                            >  1:3589713/1‑91 (MQ=255)
  tctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAttt                                                                                          <  2:191240/91‑1 (MQ=255)
  tctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAttt                                                                                          >  2:1992248/1‑91 (MQ=255)
  tctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAttt                                                                                          >  2:2512214/1‑91 (MQ=255)
  tctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCAGGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAttt                                                                                          >  2:3416745/1‑91 (MQ=255)
   ctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAtttt                                                                                         <  2:2894820/91‑1 (MQ=255)
    tttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTc                                                                                        >  2:1557600/1‑91 (MQ=255)
     ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       <  2:814091/91‑1 (MQ=255)
        tGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTAt                                                                                    >  1:3489633/1‑91 (MQ=255)
         gattCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTAtt                                                                                   <  1:342801/88‑1 (MQ=255)
           gTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGt                                                                                 >  1:2001701/1‑91 (MQ=255)
            tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                >  1:1484831/1‑91 (MQ=255)
            tCAATCCGTAATCAGTCCTCAACTCATTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                <  1:2565388/91‑1 (MQ=255)
               aTCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAAc                                                                             >  2:2980571/1‑91 (MQ=255)
                 ccGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAt                                                                           <  2:1502801/91‑1 (MQ=255)
                  cGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAtt                                                                          >  1:528422/1‑91 (MQ=255)
                   gTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTa                                                                         >  1:557334/1‑91 (MQ=255)
                    tAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaa                                                                        >  1:1668907/1‑91 (MQ=255)
                     aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       >  2:1785108/1‑91 (MQ=255)
                     aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       >  1:1739336/1‑91 (MQ=255)
                     aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       >  1:294809/1‑91 (MQ=255)
                      aTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAAt                                                                      <  2:2803889/91‑1 (MQ=255)
                      aTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAAt                                                                      <  1:3381584/91‑1 (MQ=255)
                      aTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAAt                                                                      <  1:1841232/91‑1 (MQ=255)
                       tCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATa                                                                     >  2:1146939/1‑91 (MQ=255)
                        cAGTCCTCAACTCTTTTGATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATAc                                                                    <  1:3322256/91‑1 (MQ=255)
                        cAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATAc                                                                    <  1:2122650/91‑1 (MQ=255)
                          gTCATCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAACCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGt                                                                  >  1:3083018/1‑91 (MQ=255)
                           tCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtt                                                                 <  1:2432426/91‑1 (MQ=255)
                            ccTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttt                                                                >  1:3210379/1‑91 (MQ=255)
                             cTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtttt                                                               >  2:2726445/1‑91 (MQ=255)
                               cAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtttttt                                                             >  1:2712904/1‑91 (MQ=255)
                               cAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtttttt                                                             >  2:3182574/1‑91 (MQ=255)
                                aaCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttttt                                                            >  1:69237/1‑91 (MQ=255)
                                   tctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACga                                                         <  2:1611638/91‑2 (MQ=255)
                                    ctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTa                                                        >  1:3796440/1‑91 (MQ=255)
                                    ctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTa                                                        >  2:526633/1‑91 (MQ=255)
                                      tttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaa                                                      <  1:1189607/91‑1 (MQ=255)
                                      tttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACGaaaa                                                      <  2:2746561/91‑1 (MQ=255)
                                        tAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaa                                                    >  2:1897306/1‑91 (MQ=255)
                                        tAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaa                                                    <  1:1614404/91‑1 (MQ=255)
                                         aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   >  1:3673425/1‑91 (MQ=255)
                                         aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   >  2:1836346/1‑91 (MQ=255)
                                         aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACGAaaaaaa                                                   <  2:2668702/91‑1 (MQ=255)
                                         aaTGTTTTCATGTTGTCACCGCCTTTAGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   >  2:2210472/1‑91 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  2:891978/91‑1 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  2:1571152/91‑1 (MQ=255)
                                             ttttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACttt                                               <  1:1557600/91‑1 (MQ=255)
                                             ttttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACGAAAAAAAcgta                                               <  1:1604191/91‑4 (MQ=255)
                                              tttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTa                                              <  2:431066/91‑1 (MQ=255)
                                              tttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTa                                              >  2:549965/1‑91 (MQ=255)
                                               ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  2:1799889/1‑91 (MQ=255)
                                                 cATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACta                                           >  2:2616715/1‑91 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          >  1:764255/1‑91 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          <  1:2649526/91‑1 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          >  2:3775527/1‑91 (MQ=255)
                                                   tgttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatg                                         >  2:1320968/1‑91 (MQ=255)
                                                       gtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     <  1:843345/91‑1 (MQ=255)
                                                       gtCACCGCCTTTCGTAAACCTAATATAATACGAAATCTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     <  2:3841171/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    >  1:2786424/1‑91 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:2619026/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:18854/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:2054957/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:448459/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTAATTAACAAAAAAAACTTTACtatgatatg                                    <  2:138688/91‑1 (MQ=255)
                                                         cACCGCCTTTCGTAAACCTAATATAAGACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatga                                   <  2:910725/91‑1 (MQ=255)
                                                          aCCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaa                                  <  1:593881/91‑1 (MQ=255)
                                                            cGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaaaa                                >  1:731218/1‑91 (MQ=255)
                                                             gCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAt                               <  2:2871547/91‑1 (MQ=255)
                                                              ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              <  1:322540/91‑1 (MQ=255)
                                                                tttCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTc                            >  1:1357031/1‑91 (MQ=255)
                                                                 ttCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCg                           >  2:30812/1‑91 (MQ=255)
                                                                    gTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGtat                        >  2:3039982/1‑91 (MQ=255)
                                                                     tAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGtata                       >  2:609255/1‑91 (MQ=255)
                                                                     tAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGtata                       >  1:1717524/1‑91 (MQ=255)
                                                                      aaaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataa                      <  2:2180976/91‑1 (MQ=255)
                                                                       aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:2215726/91‑1 (MQ=255)
                                                                        aCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataata                    >  1:1748163/1‑91 (MQ=255)
                                                                         ccTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTAtaataa                   <  1:55492/91‑1 (MQ=255)
                                                                         ccTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTAtaataa                   >  1:2126351/1‑91 (MQ=255)
                                                                         ccTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTAtaataa                   <  2:232104/91‑1 (MQ=255)
                                                                         ccTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTAtaataa                   <  1:91842/91‑1 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATaaaaaa               <  2:1481455/91‑1 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               >  1:670623/1‑91 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               <  1:941694/91‑1 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               >  1:305011/1‑91 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               <  2:69237/91‑1 (MQ=255)
                                                                              tataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTTGTAAAAAACTTTACTATGATATGAAAATTTCGAATAATAaaaaaa              >  2:3527652/1‑91 (MQ=255)
                                                                              tataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaa              >  2:1646501/1‑91 (MQ=255)
                                                                               ataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaaa             >  1:1033168/1‑91 (MQ=255)
                                                                                taataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaaaa            >  2:2256860/1‑91 (MQ=255)
                                                                                 aataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAag           <  2:3187119/91‑1 (MQ=255)
                                                                                 aataCGAAATTTGCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAag           <  2:3228242/91‑1 (MQ=255)
                                                                                   taCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         <  1:2458701/91‑1 (MQ=255)
                                                                                   taCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         <  1:2683116/91‑1 (MQ=255)
                                                                                   taCGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         >  2:1269355/1‑91 (MQ=255)
                                                                                     cGAAATTTTCGTATTTTCAACATTAAATACGTTTTTTTCGTAAAAAACATTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       <  2:1748634/91‑1 (MQ=255)
                                                                                     cGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       >  2:2094799/1‑91 (MQ=255)
                                                                                        aaTTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaa    <  1:2502430/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123689‑2123869

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: