Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,211,664 Δ1 bp 6.8% intergenic (+657/+156) USA300HOU_RS11275 → / ← yidC hypothetical protein/Membrane protein insertase YidC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,211,6640C. 6.8% 725.8 / 42.4 206intergenic (+657/+156)USA300HOU_RS11275/yidChypothetical protein/Membrane protein insertase YidC
Reads supporting (aligned to +/- strand):  ref base C (98/94);  new base . (7/7);  total (105/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 207 total aligned reads displayed.

GCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA  >  USA300TCH1516_ALE/2211574‑2211752
                                                                                          |                                                                                         
gCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGAc                                                                                           >  2:2172772/1‑91 (MQ=255)
  gATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACtt                                                                                         <  2:115225/91‑1 (MQ=255)
   aTGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACttt                                                                                        <  2:723893/91‑1 (MQ=255)
   aTGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACttt                                                                                        >  1:996287/1‑91 (MQ=255)
     gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                      >  2:1990859/1‑91 (MQ=255)
     gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                      >  1:2752809/1‑91 (MQ=255)
     gaagaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTca                                                                                      >  2:3172936/1‑91 (MQ=255)
       agaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTcaca                                                                                    <  2:205871/91‑1 (MQ=255)
        gaaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAt                                                                                   <  1:3779285/91‑1 (MQ=255)
         aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                                  <  1:1573813/91‑1 (MQ=255)
         aaTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAtt                                                                                  >  2:2864583/1‑91 (MQ=255)
          aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                                 <  2:2618816/91‑1 (MQ=255)
          aTGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACAttt                                                                                 <  2:2392137/91‑1 (MQ=255)
           tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                                >  2:1921324/1‑91 (MQ=255)
           tGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTg                                                                                >  1:293928/1‑91 (MQ=255)
            gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                               >  2:1310114/1‑91 (MQ=255)
            gACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGa                                                                               >  2:3355236/1‑91 (MQ=255)
             aCAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaa                                                                              <  2:1786425/91‑1 (MQ=255)
              cAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaaa                                                                             >  1:3491773/1‑91 (MQ=255)
              cAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGaaa                                                                             >  1:3491724/1‑91 (MQ=255)
                 ctctGTTTTTATCCACTTTTTACATACTGAAAACGCGTTAATATTATGTCATATAGAAACAGCAGTAAGATGATTTCAAATGTGaaagtca                                                                         <  2:1000027/91‑5 (MQ=12)
                 ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                         >  2:2112277/1‑91 (MQ=255)
                 ctctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatca                                                                         >  2:518902/1‑91 (MQ=255)
                  tctGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAatcat                                                                        <  2:1691108/91‑1 (MQ=255)
                    tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCt                                                                      <  2:1384482/91‑1 (MQ=255)
                    tGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGtcatcc                                                                      <  1:229421/91‑2 (MQ=21)
                        attaACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTAct                                                                  >  1:161390/1‑91 (MQ=255)
                          taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                <  2:2781689/91‑1 (MQ=255)
                          taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                <  2:1429446/91‑1 (MQ=255)
                          taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                <  1:3910473/91‑1 (MQ=255)
                          taACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgc                                                                >  1:3886122/1‑91 (MQ=255)
                           aaCCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgct                                                               >  1:1326570/1‑91 (MQ=255)
                            aCCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTActgctg                                                              <  2:3758978/91‑1 (MQ=255)
                             ccACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGt                                                             >  2:1211834/1‑91 (MQ=255)
                              cACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGtt                                                            <  1:1727057/91‑1 (MQ=255)
                              cACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGtt                                                            <  1:2138649/91‑1 (MQ=17)
                               aCTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttt                                                           <  1:304257/91‑1 (MQ=255)
                               aCTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGttt                                                           <  1:337939/91‑1 (MQ=255)
                                cTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTc                                                          <  2:1189645/91‑1 (MQ=255)
                                   tttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATTATCTTACTGCTGTTTCTAt                                                       >  2:69657/1‑91 (MQ=255)
                                    ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                      <  2:158450/91‑1 (MQ=255)
                                    ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                      >  2:1679845/1‑91 (MQ=255)
                                    ttACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAtt                                                      >  2:1858365/1‑91 (MQ=255)
                                     tACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAttt                                                     >  2:1845697/1‑91 (MQ=255)
                                     tACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTAttt                                                     >  1:124952/1‑91 (MQ=255)
                                       cATACTGAAAGCACGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTAt                                                   <  2:1237509/91‑1 (MQ=255)
                                        aTACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATc                                                  <  2:2407825/91‑1 (MQ=255)
                                         tACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATCTGAAATCATCTTACTGCTGTTTCTATTTATca                                                 >  1:3364436/1‑91 (MQ=255)
                                           cTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcaca                                               <  2:2563653/91‑1 (MQ=255)
                                           cTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATcaca                                               >  1:1808503/1‑91 (MQ=255)
                                                aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGGCTTTCACATTTGAAATCATCTTACTGCTGTTTCTATATATCACATATTg                                          <  2:3462046/91‑1 (MQ=255)
                                                aaCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTg                                          <  2:2730562/91‑1 (MQ=255)
                                                 aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                         <  1:544812/91‑1 (MQ=255)
                                                 aCTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGt                                         <  1:370999/91‑1 (MQ=255)
                                                   tCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtat                                       <  1:1864282/91‑1 (MQ=255)
                                                    cGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGtata                                      >  2:1614908/1‑91 (MQ=255)
                                                     gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                     >  2:2270791/1‑91 (MQ=255)
                                                     gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                     >  2:3490508/1‑91 (MQ=255)
                                                     gTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATaa                                     >  2:146981/1‑91 (MQ=255)
                                                      ttAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAt                                    >  2:310355/1‑91 (MQ=255)
                                                       tAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtg                                   >  2:468721/1‑91 (MQ=255)
                                                        aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                  <  2:864767/91‑1 (MQ=255)
                                                        aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                  <  1:3259578/91‑1 (MQ=255)
                                                        aaTATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAAtgt                                  <  2:2831971/91‑1 (MQ=255)
                                                          tattatTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGa                                <  1:959374/91‑1 (MQ=255)
                                                          tattatTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGa                                <  2:2533879/91‑1 (MQ=255)
                                                           attattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGac                               >  2:3905992/1‑91 (MQ=255)
                                                            ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                              >  1:1053891/1‑91 (MQ=255)
                                                            ttattTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGaca                              >  1:74612/1‑91 (MQ=255)
                                                              attTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACt                            >  1:3580974/1‑91 (MQ=255)
                                                               tttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTa                           <  2:2354220/91‑1 (MQ=255)
                                                                ttCAGATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                          <  2:2689392/91‑1 (MQ=255)
                                                                ttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                          <  1:2209828/91‑1 (MQ=255)
                                                                ttCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTaa                          <  1:1597400/91‑1 (MQ=255)
                                                                  cAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGt                        >  1:568899/1‑91 (MQ=255)
                                                                    aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                      >  2:802669/1‑91 (MQ=255)
                                                                    aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                      >  2:1179118/1‑91 (MQ=255)
                                                                    aaTAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGttt                      <  2:60826/91‑1 (MQ=255)
                                                                      tAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCg                    >  1:2465573/1‑91 (MQ=255)
                                                                       agaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGc                   >  1:2339208/3‑91 (MQ=11)
                                                                        gaaaCAGCAGTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTTTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt                  >  1:1158714/2‑91 (MQ=11)
                                                                        aaaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCt                  >  1:73374/1‑91 (MQ=255)
                                                                         aaaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTa                 >  1:2633553/1‑91 (MQ=255)
                                                                          aaCAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCGCATATTGTATAATGTGACACTAAGTTTCGCTAt                >  1:432940/1‑91 (MQ=255)
                                                                            cagcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTg              >  1:283664/1‑91 (MQ=255)
                                                                             agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa             >  1:3109753/1‑91 (MQ=255)
                                                                             agcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa             >  1:2558705/1‑91 (MQ=255)
                                                                             agcagTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGa             >  1:1491866/1‑91 (MQ=17)
                                                                              gcagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa            >  1:2892630/1‑91 (MQ=255)
                                                                              gcagTAAGATGATTTCAAATGTGAAAGTCATCCTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGaa            >  1:1712336/1‑91 (MQ=17)
                                                                               cagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAg           >  1:1237846/1‑91 (MQ=255)
                                                                               cagTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAg           >  2:2149931/1‑91 (MQ=255)
                                                                                agTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGc          <  1:1040281/91‑1 (MQ=255)
                                                                                 gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg         >  1:2046946/1‑91 (MQ=255)
                                                                                 gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg         <  2:329944/91‑1 (MQ=255)
                                                                                 gTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCg         <  1:3572804/91‑1 (MQ=255)
                                                                                  tAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGa        >  2:2928951/1‑91 (MQ=255)
                                                                                   aGGATGACTTTCACAGTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaa       <  2:532464/91‑1 (MQ=255)
                                                                                    ggATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGaaa      >  1:1871280/1‑91 (MQ=255)
                                                                                       tGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTAATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAaat   <  1:931672/91‑1 (MQ=255)
                                                                                          |                                                                                         
GCGATGAAGAATGACAACTCTGTTATTAACCACTTTTTACATACTGAAAACTCGTTAATATTATTTCAAATAAAAACAGCAGTAGGATGACTTTCACATTTGAAATCATCTTACTGCTGTTTCTATTTATCACATATTGTATAATGTGACACTAAGTTTCGCTATTGAAGCGAAAAATA  >  USA300TCH1516_ALE/2211574‑2211752

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: