Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,139,620 T→A 9.3% intergenic (‑60/+851) USA300HOU_RS10910 ← / ← groL hypothetical protein/60 kDa chaperonin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,139,6200TA 9.3% 659.5 / 44.1 204intergenic (‑60/+851)USA300HOU_RS10910/groLhypothetical protein/60 kDa chaperonin
Reads supporting (aligned to +/- strand):  ref base T (102/82);  new base A (5/14);  total (107/97)
Fisher's exact test for biased strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 207 total aligned reads displayed.

AACGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCTG  >  USA300TCH1516_ALE/2139530‑2139710
                                                                                          |                                                                                          
aaCGTGCTTGCTTTGCAGTTAACTCTCTCATTTTATCACCCCCATAATTTTTTCTAATATGGTTTCATACCATAATATTACAGATTGTTCt                                                                                            <  2:3713056/91‑1 (MQ=255)
 aCGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATAGCATAATATTACAGATTGTTCTg                                                                                           >  2:3125477/1‑91 (MQ=255)
 aCGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTg                                                                                           >  1:3560381/1‑91 (MQ=255)
 aCGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTg                                                                                           >  1:2842711/1‑91 (MQ=255)
  cGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCTTACCATAATATTACAGATTGTTCAGa                                                                                          <  2:2347868/91‑1 (MQ=255)
  cGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGa                                                                                          <  1:1736276/91‑1 (MQ=255)
   gTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCAGaa                                                                                         <  2:1037269/91‑1 (MQ=255)
    tgcttgctTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTTTCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAAc                                                                                        <  2:2989959/91‑1 (MQ=255)
     gcttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACACTTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACa                                                                                       >  2:654274/1‑91 (MQ=255)
     gcttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACa                                                                                       >  2:1734650/1‑91 (MQ=255)
     gcttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACa                                                                                       >  1:2298417/1‑91 (MQ=255)
     gcttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACa                                                                                       >  2:2055451/1‑91 (MQ=255)
     gcttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACa                                                                                       >  1:1354223/1‑91 (MQ=255)
      cttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCAGAACaa                                                                                      <  1:1101732/91‑1 (MQ=255)
       ttgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAAt                                                                                     <  1:1258492/91‑1 (MQ=255)
        tgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATc                                                                                    <  2:2205707/91‑1 (MQ=255)
        tgcttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATc                                                                                    >  1:2821205/1‑91 (MQ=255)
          cttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTa                                                                                  <  1:1573775/91‑1 (MQ=255)
          cttTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTa                                                                                  >  1:689325/1‑91 (MQ=255)
           tttGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCGAATATGGTTTCATACCATAATATTACAGATTGTTCAGAACAATCtgt                                                                                 <  1:3883305/91‑3 (MQ=255)
              gCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGc                                                                              <  2:275821/91‑1 (MQ=255)
               cAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCa                                                                             >  2:1547523/1‑91 (MQ=255)
                aGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCac                                                                            <  1:53179/91‑1 (MQ=255)
                aGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCac                                                                            >  1:471751/1‑91 (MQ=255)
                  ttAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCacta                                                                          >  1:2784052/1‑91 (MQ=255)
                  ttAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCacta                                                                          <  2:413358/91‑1 (MQ=255)
                   tAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCactac                                                                         <  1:2562018/91‑1 (MQ=255)
                   tAACTCACGCATGTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCGGAACAATCTAAGGCactac                                                                         <  2:1646118/91‑1 (MQ=255)
                    aaCTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCactact                                                                        <  2:3214476/91‑1 (MQ=255)
                      ctcactcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTaa                                                                      <  2:3632319/91‑1 (MQ=255)
                       tcactcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAat                                                                     >  1:3833441/1‑91 (MQ=255)
                        cactcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAata                                                                    >  1:860753/1‑91 (MQ=255)
                         actcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAatat                                                                   <  1:1024997/91‑1 (MQ=255)
                          ctcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATc                                                                  >  2:1066218/1‑91 (MQ=255)
                          ctcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATc                                                                  <  1:765744/91‑1 (MQ=255)
                          ctcaTTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATc                                                                  >  1:3036969/1‑91 (MQ=255)
                               tttATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGAt                                                             >  1:1617774/1‑91 (MQ=255)
                               tttATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGAt                                                             <  1:1514212/91‑1 (MQ=255)
                                ttATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATa                                                            >  2:2498113/1‑91 (MQ=255)
                                  aTCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAAc                                                          >  2:966176/1‑91 (MQ=255)
                                  aTCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAAc                                                          <  1:2233044/91‑1 (MQ=255)
                                  aTCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAAc                                                          >  2:2703360/1‑91 (MQ=255)
                                  aTCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAAc                                                          >  1:2525293/1‑91 (MQ=255)
                                  aTCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAAc                                                          <  2:3472729/91‑1 (MQ=255)
                                   tCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACt                                                         <  1:1465125/91‑1 (MQ=255)
                                   tCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACt                                                         >  1:470545/1‑91 (MQ=255)
                                   tCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACt                                                         <  1:2960730/91‑1 (MQ=255)
                                     aCCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTaa                                                       >  1:916141/1‑91 (MQ=255)
                                     aCCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTaa                                                       >  1:591070/1‑91 (MQ=255)
                                     aCCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTaa                                                       <  1:1332432/91‑1 (MQ=255)
                                      cccccACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAg                                                      >  2:2460551/1‑91 (MQ=255)
                                       ccccACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGt                                                     >  1:3413297/1‑91 (MQ=255)
                                       ccccACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGt                                                     >  1:1627459/1‑91 (MQ=255)
                                          cacaATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCt                                                  <  1:3191571/91‑1 (MQ=255)
                                            caATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTg                                                >  2:3115075/1‑91 (MQ=255)
                                             aaTTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATATCTAAGTCCTTGt                                               >  2:2225156/1‑91 (MQ=255)
                                              aTTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGta                                              <  1:2763104/91‑1 (MQ=255)
                                              aTTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGta                                              <  1:1711260/91‑1 (MQ=255)
                                                 ttATCTGATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATaa                                           <  2:318579/91‑1 (MQ=255)
                                                  tATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAg                                          >  2:1883187/1‑91 (MQ=255)
                                                  tATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAg                                          >  2:1956544/1‑91 (MQ=255)
                                                    tCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAggg                                        <  1:2685538/91‑1 (MQ=255)
                                                          atGGTATGAAACCATAATATTACAGATTGTTCTGCACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCaa                                  >  2:1162622/1‑91 (MQ=255)
                                                            ggTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCaaaa                                >  1:972393/1‑91 (MQ=255)
                                                              tttCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAat                              >  2:465562/1‑91 (MQ=255)
                                                              tttCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGAGAACTAAGTCCTTGTATAAGGGAGTCAAAAat                              <  1:3213356/91‑1 (MQ=255)
                                                               ttCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAata                             <  1:1145575/91‑1 (MQ=255)
                                                                tCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAatat                            >  1:2044525/1‑91 (MQ=255)
                                                                tCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAatat                            <  2:3211732/91‑1 (MQ=255)
                                                                tCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAatat                            >  1:1810669/1‑91 (MQ=255)
                                                                tCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAatat                            >  2:182925/1‑91 (MQ=255)
                                                                 cATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAAtata                           >  2:864787/1‑91 (MQ=255)
                                                                    aCCATGATATTACAGATTGTTCTGAGCAATCTAAGGCACCACTAATATCTTGAAAACTAAGTCCTTGTATAAGGGAGTCAAAAATATaaaa                        <  1:3609738/91‑1 (MQ=255)
                                                                       ataataTTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAActc                     >  1:3455492/1‑91 (MQ=255)
                                                                         aataTTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCtt                   <  2:517393/91‑1 (MQ=255)
                                                                         aataTTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCtt                   <  1:141503/91‑1 (MQ=255)
                                                                          atatTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCttt                  >  1:3194016/1‑91 (MQ=255)
                                                                          atatTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCttt                  <  2:144377/91‑1 (MQ=255)
                                                                            atTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTAt                >  2:333492/1‑91 (MQ=255)
                                                                             ttACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATc               >  1:3175854/1‑91 (MQ=255)
                                                                               aCAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCgg             >  1:3251888/1‑91 (MQ=255)
                                                                               aCAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCgg             >  2:327536/1‑91 (MQ=255)
                                                                                gggATTGTTCTGAACAATCGAAGGCACTACGAATATCTGGATAACTAAGGCCTTGTGTAAGGGAGTGAGAAATATAAAACTCTTTATCGGt            <  2:3565271/89‑1 (MQ=255)
                                                                                 aGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTc           >  2:360973/1‑91 (MQ=255)
                                                                                  gATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCg          >  1:2939747/1‑91 (MQ=255)
                                                                                  gATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCg          <  1:1256328/91‑1 (MQ=255)
                                                                                  gATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCg          <  2:3036366/91‑1 (MQ=255)
                                                                                   aTTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGc         <  2:1342476/91‑1 (MQ=255)
                                                                                   aTTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGc         <  2:517853/91‑1 (MQ=255)
                                                                                    ttGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCt        >  2:2224200/1‑91 (MQ=255)
                                                                                     tGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTa       >  1:3693292/1‑91 (MQ=255)
                                                                                     tGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTa       <  1:3244755/91‑1 (MQ=255)
                                                                                      gTTCAGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTaa      >  1:2436777/1‑91 (MQ=255)
                                                                                      gTTCAGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTaa      >  2:1673955/1‑91 (MQ=255)
                                                                                      gTTCAGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTaa      >  2:693645/1‑91 (MQ=255)
                                                                                       ttCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAAt     <  2:2843279/91‑1 (MQ=255)
                                                                                         cTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCt   >  2:1210564/1‑91 (MQ=255)
                                                                                         cTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCt   >  1:3329541/1‑91 (MQ=255)
                                                                                          tGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCTg  >  2:2271659/1‑91 (MQ=255)
                                                                                          tGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCTg  >  1:3103497/1‑91 (MQ=255)
                                                                                          |                                                                                          
AACGTGCTTGCTTTGCAGTTAACTCACTCATTTTATCACCCCCACAATTTTATCTAATATGGTTTCATACCATAATATTACAGATTGTTCTGAACAATCTAAGGCACTACTAATATCTTGATAACTAAGTCCTTGTATAAGGGAGTCAAAAATATAAAACTCTTTATCGGTCGCTAATCTG  >  USA300TCH1516_ALE/2139530‑2139710

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: