Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,123,779 T→A 9.9% 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 9.9% 540.2 / 40.0 162intergenic (‑39/‑94)USA300HOU_RS10825/lexA_2hypothetical protein/LexA repressor
Reads supporting (aligned to +/- strand):  ref base T (79/67);  new base A (5/11);  total (84/78)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 171 total aligned reads displayed.

ACTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123689‑2123869
                                                                                          |                                                                                          
  tctTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAttt                                                                                          <  1:1504413/91‑1 (MQ=255)
     ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       >  1:3810888/1‑91 (MQ=255)
     ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       <  1:3952646/91‑1 (MQ=255)
     ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       <  1:4042260/91‑1 (MQ=255)
     ttttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCg                                                                                       <  1:721435/91‑1 (MQ=255)
      tttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGt                                                                                      >  1:1025475/1‑91 (MQ=255)
      tttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGt                                                                                      <  1:716794/91‑1 (MQ=255)
      tttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGt                                                                                      <  1:569671/91‑1 (MQ=255)
      tttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTAGTAAACCTAATATAATACGAAAATTTCGt                                                                                      <  1:186901/91‑1 (MQ=255)
       ttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTa                                                                                     >  1:3125990/1‑91 (MQ=255)
       ttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTa                                                                                     >  1:1175337/1‑91 (MQ=255)
       ttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTa                                                                                     <  1:3088225/91‑1 (MQ=255)
       ttGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTa                                                                                     <  1:879477/91‑1 (MQ=255)
        tGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTAt                                                                                    <  1:954254/91‑1 (MQ=255)
        tGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTAt                                                                                    <  1:953058/91‑1 (MQ=255)
         gAGTCAATCCGTAATCAGTCCTCAACTCTGTAAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTAtt                                                                                   <  1:3319714/91‑1 (MQ=255)
            tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                <  1:655573/91‑1 (MQ=255)
            tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTc                                                                                >  1:2953101/1‑91 (MQ=255)
            tCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAAATTTCGTAttata                                                                                <  1:1256288/91‑4 (MQ=255)
             cAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCa                                                                               >  1:853663/1‑91 (MQ=255)
              aaTCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTCTTGTCaa                                                                              >  1:2070905/1‑91 (MQ=255)
                 ccGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAt                                                                           >  1:1123209/1‑91 (MQ=255)
                 ccGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAt                                                                           >  1:1123943/1‑91 (MQ=255)
                 ccGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAt                                                                           <  1:3792197/91‑1 (MQ=255)
                  cGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAtt                                                                          >  1:1712566/1‑91 (MQ=255)
                  cGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACAtt                                                                          >  1:468967/1‑91 (MQ=255)
                   gTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTa                                                                         >  1:17393/1‑91 (MQ=255)
                    tAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaa                                                                        >  1:229418/1‑91 (MQ=255)
                     aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       >  1:3954871/1‑91 (MQ=255)
                     aaTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTaaa                                                                       >  1:1427358/1‑91 (MQ=255)
                      aTCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAAt                                                                      <  1:2258837/91‑1 (MQ=255)
                       tCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATa                                                                     <  1:1874844/91‑1 (MQ=255)
                        cAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATAc                                                                    >  1:2490255/1‑91 (MQ=255)
                        cAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATAc                                                                    >  1:2331354/1‑91 (MQ=255)
                         aGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACg                                                                   <  1:315444/91‑1 (MQ=255)
                           tCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtc                                                                 <  1:1918675/91‑2 (MQ=255)
                            ccTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttt                                                                <  1:919675/91‑1 (MQ=255)
                              tCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttt                                                              <  1:3926858/91‑1 (MQ=255)
                               cAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGtttttt                                                             >  1:1310162/1‑91 (MQ=255)
                                aaCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttttt                                                            <  1:133358/91‑1 (MQ=255)
                                aaCTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGttttttt                                                            >  1:2883120/1‑91 (MQ=255)
                                   tctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGt                                                         <  1:3202686/91‑1 (MQ=255)
                                   tctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGt                                                         <  1:3202181/91‑1 (MQ=255)
                                   tctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGt                                                         <  1:1827712/91‑1 (MQ=255)
                                    ctTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTa                                                        <  1:1544760/91‑1 (MQ=255)
                                     ttttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaa                                                       <  1:1607207/91‑1 (MQ=255)
                                       ttAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaa                                                     >  1:3613538/1‑91 (MQ=255)
                                         aaTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTaaaaaa                                                   <  1:603468/91‑1 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  1:254209/91‑1 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  >  1:9240/1‑91 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  1:1769627/91‑1 (MQ=255)
                                          aTGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAAc                                                  <  1:2665718/91‑1 (MQ=255)
                                            gTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACtt                                                <  1:3652018/91‑1 (MQ=255)
                                             ttttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACttt                                               >  1:2639524/1‑91 (MQ=255)
                                             ttttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTACGAAAAAAAcgta                                               <  1:2130661/91‑4 (MQ=255)
                                               ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  1:242083/1‑91 (MQ=255)
                                               ttCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTAc                                             >  1:2175801/1‑91 (MQ=255)
                                                 cATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTATAAAAATTTACta                                           >  1:1131702/1‑91 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          >  1:3951655/1‑91 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          >  1:770356/1‑91 (MQ=255)
                                                  aTGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtat                                          >  1:2968798/1‑91 (MQ=255)
                                                   tgttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatg                                         >  1:3693690/1‑91 (MQ=255)
                                                    gttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatga                                        >  1:1606487/1‑91 (MQ=255)
                                                     ttgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgat                                       <  1:962954/91‑1 (MQ=255)
                                                      tgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgata                                      >  1:3257804/1‑91 (MQ=255)
                                                      tgtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgata                                      <  1:3139413/91‑1 (MQ=255)
                                                       gtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     <  1:3674051/91‑1 (MQ=255)
                                                       gtCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatat                                     <  1:1224302/91‑1 (MQ=255)
                                                        tCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACtatgatatg                                    <  1:81740/91‑1 (MQ=255)
                                                          aCCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGaa                                  <  1:3445797/91‑1 (MQ=255)
                                                             gCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAt                               >  1:2190589/1‑91 (MQ=255)
                                                             gCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAt                               >  1:3627211/1‑91 (MQ=255)
                                                              ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              >  1:1413578/1‑91 (MQ=255)
                                                              ccTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAAtt                              <  1:3776946/91‑1 (MQ=255)
                                                                tttCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTc                            >  1:3367404/1‑91 (MQ=255)
                                                                    gTAAACCTAATATAATACGAAATTTTCGTATTGTCAATATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGtat                        >  1:3509165/1‑91 (MQ=255)
                                                                       aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGtattat                     >  1:996503/1‑88 (MQ=255)
                                                                       aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:467922/91‑1 (MQ=255)
                                                                       aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:1304635/91‑1 (MQ=255)
                                                                       aaCCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:3483789/91‑1 (MQ=255)
                                                                       aaCCTAATAGAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTataat                     <  1:3963359/91‑1 (MQ=255)
                                                                         ccTAATATAATGCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTAtaataa                   <  1:1382756/91‑1 (MQ=255)
                                                                          cTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAg                  <  1:2470906/91‑1 (MQ=255)
                                                                          cTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAg                  <  1:2803015/91‑1 (MQ=255)
                                                                           taatataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAagaaga                 <  1:3278631/91‑1 (MQ=255)
                                                                            aatataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaa                >  1:1787789/1‑91 (MQ=255)
                                                                             atataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaa               >  1:1015451/1‑91 (MQ=255)
                                                                              tataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGTAAATTTCGTATAATAAGaaaa              >  1:726111/1‑91 (MQ=255)
                                                                              tataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaa              >  1:2444427/1‑91 (MQ=255)
                                                                              tataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGAGATGAAAATTTCGTATAATAAGaaaa              >  1:3864511/1‑91 (MQ=255)
                                                                               ataataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGaaaaa             <  1:1351/91‑1 (MQ=255)
                                                                                  ataCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagg          <  1:1819052/91‑1 (MQ=255)
                                                                                   taCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         >  1:394688/1‑91 (MQ=255)
                                                                                   taCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         >  1:1452001/1‑91 (MQ=255)
                                                                                   taCGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCATAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAagga         >  1:3487366/1‑91 (MQ=255)
                                                                                     cGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       <  1:3895395/91‑1 (MQ=255)
                                                                                     cGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       >  1:944643/1‑91 (MQ=255)
                                                                                     cGAAAATTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAaggagg       >  1:3638756/1‑91 (MQ=255)
                                                                                        aaTTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaa    >  1:550491/1‑91 (MQ=255)
                                                                                          ttttCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGgtaagt  <  1:48865/91‑1 (MQ=255)
                                                                                          |                                                                                          
ACTCTTTTTGAGTCAATCCGTAATCAGTCCTCAACTCTTTTAATGTTTTCATGTTGTCACCGCCTTTCGTAAACCTAATATAATACGAAATTTTCGTATTGTCAACATTAAATACGTTTTTTTCGTAAAAAACTTTACTATGATATGAAAATTTCGTATAATAAGAAAAAAGGAGGTAAGT  >  USA300TCH1516_ALE/2123689‑2123869

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