Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 770,888 C→T 6.3% intergenic (+74/‑176) ybaK → / → glcR Cys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE770,8880CT 6.3% 781.9 / 24.3 235intergenic (+74/‑176)ybaK/glcRCys‑tRNA(Pro)/Cys‑tRNA(Cys) deacylase YbaK/HTH‑type transcriptional repressor GlcR
Reads supporting (aligned to +/- strand):  ref base C (105/115);  new base T (10/5);  total (115/120)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 238 total aligned reads displayed.

AGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTCA  >  USA300TCH1516_ALE/770799‑770978
                                                                                           |                                                                                          
aGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGcc                                                                                           <  2:2740997/91‑1 (MQ=255)
 gTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCca                                                                                          <  2:398888/91‑1 (MQ=255)
  ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                                         <  2:3872570/91‑1 (MQ=255)
  ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                                         >  2:1612294/1‑91 (MQ=255)
  ttATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCcac                                                                                         <  1:242860/91‑1 (MQ=255)
   tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTcaca                                                                                        <  1:3571723/91‑1 (MQ=14)
   tATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTcaca                                                                                        <  1:2057318/91‑1 (MQ=17)
    aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCacac                                                                                       <  2:2062024/91‑1 (MQ=255)
    aTCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCacac                                                                                       <  1:2306590/91‑1 (MQ=255)
     tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACt                                                                                      <  2:3746455/91‑1 (MQ=255)
     tCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACt                                                                                      <  1:3113137/91‑1 (MQ=255)
      ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTATGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTc                                                                                     <  2:1425979/91‑1 (MQ=255)
      ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTc                                                                                     <  2:3450600/91‑1 (MQ=14)
      ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTc                                                                                     <  1:1824883/91‑1 (MQ=255)
      ccATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTc                                                                                     <  1:836085/91‑1 (MQ=255)
        atgaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCt                                                                                   <  2:2117609/91‑1 (MQ=255)
          gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                                                 >  2:2363890/1‑91 (MQ=255)
          gaatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCttt                                                                                 >  1:936538/1‑91 (MQ=255)
           aatgaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCtttt                                                                                >  1:3803276/1‑91 (MQ=255)
              gaTTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttga                                                                             <  1:1065641/91‑1 (MQ=255)
                 taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgattg                                                                          >  2:2668342/1‑91 (MQ=255)
                 taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgattg                                                                          >  2:19544/1‑91 (MQ=255)
                 taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttgattg                                                                          >  1:3186472/1‑91 (MQ=255)
                 taataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTttg                                                                              <  2:2355851/87‑1 (MQ=255)
                   ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAGAT‑‑TTCATTTCTGCCACACTCCTTTTTGATTGaa                                                                        <  2:2876691/91‑1 (MQ=255)
                   ataaCAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGaa                                                                        <  1:1576720/91‑1 (MQ=255)
                     aacaacaacAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAAtt                                                                      >  1:1634004/1‑91 (MQ=255)
                       caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                                                    <  1:620742/91‑1 (MQ=255)
                       caacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAg                                                                    <  1:3312877/91‑1 (MQ=255)
                        aacaaAAAGTATGGGGCAAGATTAGGAGTGTGGCAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                                                   <  2:2826343/91‑1 (MQ=255)
                        aacaaAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGc                                                                   <  1:1954382/91‑1 (MQ=255)
                            aaaaGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCAttt                                                               <  1:2850817/91‑1 (MQ=255)
                               aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTAc                                                            >  2:3069897/1‑91 (MQ=255)
                               aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                            >  2:1534218/1‑91 (MQ=255)
                               aGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTAc                                                            <  1:1071285/91‑1 (MQ=255)
                                  aTGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAt                                                         >  1:2512506/1‑91 (MQ=255)
                                   tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                        <  2:1822485/91‑1 (MQ=255)
                                   tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATc                                                        <  2:2921897/91‑1 (MQ=255)
                                   tGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACCCCTTTTTGATTGAATTAGCATTTTACGATc                                                        <  2:1559250/91‑1 (MQ=255)
                                    ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                       >  2:1668204/1‑91 (MQ=255)
                                    ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                       >  1:951838/1‑91 (MQ=255)
                                    ggggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCa                                                       >  2:1047619/1‑91 (MQ=255)
                                     gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                      <  2:657186/91‑1 (MQ=255)
                                     gggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCAt                                                      <  2:1393197/91‑1 (MQ=255)
                                      ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                     <  2:2148640/91‑1 (MQ=255)
                                      ggCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATa                                                     <  1:1104513/91‑1 (MQ=255)
                                        cAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATaaa                                                   <  2:51549/91‑1 (MQ=255)
                                             ttAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                              <  2:3711233/91‑1 (MQ=255)
                                             ttAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTc                                              >  2:3589922/1‑91 (MQ=255)
                                             ttAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGAACATAAACAGTc                                              >  1:3510308/1‑91 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTTTTTTTATTAGCATTTTACGATCATAAACAGTCAt                                            >  2:909854/1‑91 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                            <  2:1260239/91‑1 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                            <  2:2410400/91‑1 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                            >  2:409189/1‑91 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                            >  2:677235/1‑91 (MQ=255)
                                               aGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCAt                                            >  1:2365373/1‑91 (MQ=255)
                                                 gAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATta                                          >  1:1454316/1‑91 (MQ=255)
                                                  aGTGTGGCAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGTATTAGCATTTTACGATCATAAACAGTCATtat                                         >  2:2293872/1‑91 (MQ=14)
                                                  aGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATtat                                         <  2:1732009/91‑1 (MQ=255)
                                                    tgtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATaa                                       <  2:2456269/91‑1 (MQ=255)
                                                     gtggcAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                      >  1:2888936/5‑91 (MQ=14)
                                                     gtgACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAAt                                      <  2:3402453/91‑1 (MQ=255)
                                                       ggcAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                                    >  2:2321562/3‑91 (MQ=17)
                                                       gACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTg                                    <  2:369003/91‑1 (MQ=255)
                                                         cAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                                  <  1:861489/91‑1 (MQ=255)
                                                         cAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGGCATTATAATTGAg                                  >  1:1757239/1‑91 (MQ=255)
                                                         cAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                                  >  2:604640/1‑91 (MQ=18)
                                                         cAGAAATGAATTTTATAAAAAT‑‑TTCATCTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAg                                  >  2:3037331/1‑91 (MQ=21)
                                                          agagATGAAATTTTTATAAAA‑‑TTCATTTCTGCCCCACACCTTTTTGCTTGAATTAGCATCTTACGATCATAAACAGTCATTATAATTGAGt                                 >  1:3093321/1‑91 (MQ=255)
                                                          agagATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                                 <  1:1298986/91‑1 (MQ=255)
                                                          agaaaTGAATTTTATAAAAAT‑‑TTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGt                                 >  2:3535043/5‑91 (MQ=21)
                                                            agaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAt                               >  2:3308908/1‑91 (MQ=255)
                                                             gaTGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTAtt                              >  1:1182461/1‑91 (MQ=255)
                                                               tGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTg                            >  2:36355/1‑91 (MQ=255)
                                                                 aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                          >  2:80216/1‑91 (MQ=255)
                                                                 aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                          <  2:2862882/91‑1 (MQ=255)
                                                                 aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                          <  1:470648/91‑1 (MQ=255)
                                                                 aaaTTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                          >  1:2562273/1‑91 (MQ=255)
                                                                   aattttATAAAAATTTCATCTCTGTCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGaa                          >  1:2786620/3‑91 (MQ=25)
                                                                    tttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt                       >  1:144697/1‑91 (MQ=255)
                                                                    tttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACAt                       <  1:2668342/91‑1 (MQ=255)
                                                                     ttttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATa                      >  2:1079356/1‑91 (MQ=255)
                                                                       ttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa                    <  1:346043/91‑1 (MQ=255)
                                                                       ttATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaa                    <  1:1737600/91‑1 (MQ=255)
                                                                         ataAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATaaaaa                  >  1:247187/1‑91 (MQ=255)
                                                                          taAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAAt                 <  1:208215/91‑1 (MQ=255)
                                                                           aaaa‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATg                <  2:1119650/91‑1 (MQ=255)
                                                                              a‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTaa             >  1:1274769/1‑91 (MQ=255)
                                                                              a‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTaa             >  2:2083937/1‑91 (MQ=255)
                                                                                 ttCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAt            <  2:950016/91‑1 (MQ=255)
                                                                                  tCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAtt           <  1:569473/91‑1 (MQ=255)
                                                                                  tCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAAtt           >  1:2921497/1‑91 (MQ=255)
                                                                                    aTTTCTGCCACACTCCTTTTTTTTTGTTTTTTCATTTTACGATCATAAACATTCATTATAATTGAGTATTTGAACATAAAAATGTAAtttt         >  2:278875/1‑91 (MQ=255)
                                                                                      ttCTGCCACACTCCTTTTTGCTTGAATTAGCATTTTACGATCATAAACAGTCATAATAATTGAGTATTTGAACATAAAAATGTAATTTTac       >  2:3835593/1‑90 (MQ=255)
                                                                                      ttCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTAt       >  1:164572/1‑91 (MQ=255)
                                                                                       tCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATc      >  2:3515779/1‑91 (MQ=255)
                                                                                       tCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATc      <  2:1462310/91‑1 (MQ=255)
                                                                                        cTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCg     <  2:2562273/91‑1 (MQ=255)
                                                                                          gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc   >  2:3035904/1‑91 (MQ=255)
                                                                                          gCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTc   <  1:2344669/91‑1 (MQ=255)
                                                                                           |                                                                                          
AGTTATCCATGAATGATTAATAACAACAAAAAGTATGGGGCAAGATTAGGAGTGTGACAGAGATGAAATTTTTATAAAA‑‑TTCATTTCTGCCACACTCCTTTTTGATTGAATTAGCATTTTACGATCATAAACAGTCATTATAATTGAGTATTTGAACATAAAAATGTAATTTTATCGTCA  >  USA300TCH1516_ALE/770799‑770978

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: