Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 593,433 T→A 8.2% intergenic (+25/‑112) rpoC → / → rplGB DNA‑directed RNA polymerase subunit beta'/Ribosome‑associated protein L7Ae‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE593,4330TA 8.2% 803.3 / 42.5 251intergenic (+25/‑112)rpoC/rplGBDNA‑directed RNA polymerase subunit beta'/Ribosome‑associated protein L7Ae‑like protein
Reads supporting (aligned to +/- strand):  ref base T (137/93);  new base A (10/11);  total (147/104)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 263 total aligned reads displayed.

GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGC  >  USA300TCH1516_ALE/593343‑593523
                                                                                           |                                                                                          
gTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCt                                                                                            <  2:397407/91‑1 (MQ=255)
 tAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCtt                                                                                           <  2:2092471/91‑1 (MQ=255)
 aaaaaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCtt                                                                                           <  2:1155043/90‑1 (MQ=255)
  aaaaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCttt                                                                                          <  2:685976/91‑1 (MQ=255)
   acaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAAGCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGCGGCTAAAGCatt                                                                                          <  2:2023733/89‑1 (MQ=255)
    aaaaCGAAAAAACAGCTACACCGGTTGCAGAAGGTGAATCCCAAACTGACGTAACGGAATGACAAGTATATAACAGAGGCTAATGCTTTAg                                                                                        <  2:2991850/88‑1 (MQ=255)
    aaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAg                                                                                        <  2:2833365/91‑1 (MQ=255)
     aaaCGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAg                                                                                        <  2:1544413/89‑1 (MQ=255)
      tACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACACGTATATAACAGAGGCTAATGCTTTAgac                                                                                      >  1:1902308/1‑89 (MQ=255)
        cGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctc                                                                                    <  2:2403002/91‑1 (MQ=255)
        aaaaaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctc                                                                                    <  2:3487996/89‑1 (MQ=255)
         gAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctct                                                                                   <  1:1086468/91‑1 (MQ=255)
          aaaaaaCAGCTAGACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctct                                                                                   <  2:637451/91‑1 (MQ=38)
          aaaaaaCAGCTAAGCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctct                                                                                   <  2:2725174/91‑1 (MQ=38)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                 <  1:2928257/91‑1 (MQ=255)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                 >  1:689016/1‑91 (MQ=255)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                 >  1:2284525/1‑91 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                <  1:1644468/91‑1 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                >  2:2821240/1‑91 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                <  1:2453941/91‑1 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                                >  2:2153831/1‑91 (MQ=255)
             aaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGtt                                                                               >  1:2130643/1‑91 (MQ=255)
               ccggCTAAACCAGGTGCAGAAGTGGAATCTCAAACTGAAGTAGGGGAATAACAAGTATAGAACAAGGGGCTAAAGCATTAGCCTCTGTTAt                                                                             <  2:1615081/88‑1 (MQ=16)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                          >  1:743954/1‑91 (MQ=255)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                          >  2:2418912/1‑91 (MQ=255)
                   cTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAttttt                                                                         <  2:38279/91‑1 (MQ=255)
                       aCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGt                                                                     >  1:48341/1‑91 (MQ=255)
                         cAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTaa                                                                   >  2:3422862/1‑91 (MQ=255)
                           gTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAat                                                                 <  2:3281135/91‑1 (MQ=255)
                            ttGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatt                                                                >  1:2218667/1‑91 (MQ=255)
                             tGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatta                                                               >  2:2126871/1‑91 (MQ=255)
                              gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                              >  1:2108766/1‑91 (MQ=255)
                              gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                              <  1:3593621/91‑1 (MQ=255)
                              gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                              >  1:1236469/1‑91 (MQ=255)
                               cAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattatt                                                             >  2:365839/1‑91 (MQ=255)
                                aGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTAttt                                                            >  1:294868/1‑91 (MQ=255)
                                 gAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTg                                                           >  2:495216/1‑91 (MQ=255)
                                  aaGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGa                                                          <  2:3390709/91‑1 (MQ=255)
                                   aGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAt                                                         >  2:1744602/1‑91 (MQ=255)
                                    gTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAtt                                                        >  1:2994139/1‑91 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                       <  1:1539020/91‑1 (MQ=255)
                                       gAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTaa                                                     >  1:2026066/1‑91 (MQ=255)
                                        aaTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAAt                                                    <  2:476102/91‑1 (MQ=255)
                                         aTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATg                                                   >  1:1214274/1‑91 (MQ=255)
                                          tctcAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGt                                                  <  1:2829324/91‑1 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGCTTTAATGTTGa                                               >  2:3768540/1‑91 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                               >  2:3727792/1‑91 (MQ=255)
                                              aaaCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGAc                                              >  2:2793107/1‑91 (MQ=255)
                                               aaCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTACGTAAATTATTTGATTTAATGTTGACg                                             <  1:3241006/91‑1 (MQ=255)
                                                aCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGGCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGa                                            >  1:3949555/1‑91 (MQ=255)
                                                aCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGa                                            >  1:1773121/1‑91 (MQ=255)
                                                 cTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGaa                                           <  1:2362852/91‑1 (MQ=255)
                                                 cTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGaa                                           <  1:343813/91‑1 (MQ=255)
                                                    aaGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtc                                        >  2:2704197/1‑91 (MQ=255)
                                                     aGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAACTTATTTGATTTAATGTCGACGAATtct                                       >  2:1887073/1‑91 (MQ=255)
                                                     aGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtct                                       >  2:3227507/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  1:3013606/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      <  1:1435399/91‑1 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  1:2924980/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  1:525457/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                      >  2:1004487/1‑91 (MQ=255)
                                                       tAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTctct                                     >  2:1548773/1‑91 (MQ=255)
                                                       tAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTctct                                     <  1:1547560/91‑1 (MQ=255)
                                                        aaCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCtt                                    <  1:592838/91‑1 (MQ=255)
                                                         aCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTg                                   >  2:1146679/1‑91 (MQ=255)
                                                             aaTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCa                               >  2:1958213/1‑91 (MQ=255)
                                                             aaTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCa                               >  2:1294189/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                            >  1:912137/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                            >  2:1836244/1‑91 (MQ=255)
                                                                 aCAAGGATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGt                           <  1:1316425/91‑1 (MQ=255)
                                                                  cAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGtt                          >  1:1471941/1‑91 (MQ=255)
                                                                  cAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGtt                          >  2:3392429/1‑91 (MQ=255)
                                                                      tatataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATTTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTActa                      >  2:919675/1‑91 (MQ=255)
                                                                      aaaataACAAGAGGCTAAAGCATTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAat                       >  1:1842836/4‑91 (MQ=37)
                                                                        tataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatata                    <  2:2607354/91‑1 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                   <  1:431568/91‑1 (MQ=255)
                                                                          taACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATAtt                  <  2:3861826/91‑1 (MQ=255)
                                                                          taACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATGATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATAtt                  <  1:3155538/91‑1 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGCTTTAACGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 >  2:3871789/1‑91 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCATTGTTTATATATTa                 >  2:2554189/1‑91 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                 >  1:3829919/1‑91 (MQ=255)
                                                                            acaagagGCTAAAGCATTAGCCTCTGTTATTTTAATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                >  1:1193233/4‑91 (MQ=37)
                                                                            aCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                <  2:3576401/91‑1 (MQ=255)
                                                                            aCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                <  2:1580746/91‑1 (MQ=255)
                                                                             caagagGCTAAAGCATTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa                >  1:2662145/3‑91 (MQ=39)
                                                                             cAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaaa               <  1:1612768/91‑1 (MQ=255)
                                                                                gagGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAgg             >  2:227577/1‑91 (MQ=255)
                                                                                gagGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTAACGAATTCTCTTGTTCAATGTTACTATATTAAAgg             >  2:2447798/1‑91 (MQ=255)
                                                                                  ggCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGtt           >  2:778725/1‑91 (MQ=255)
                                                                                   gCTAAAGCATTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTg          >  2:23069/1‑91 (MQ=255)
                                                                                   gCTAAAGCATTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTg          >  1:3259575/1‑91 (MQ=255)
                                                                                    cTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGa         >  2:270593/1‑91 (MQ=255)
                                                                                    cTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGa         <  2:1481329/91‑1 (MQ=255)
                                                                                     tAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGAt        >  2:1345444/1‑91 (MQ=255)
                                                                                      aaTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATg       >  1:2505327/1‑91 (MQ=255)
                                                                                      aaTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATg       >  1:1425386/1‑91 (MQ=255)
                                                                                       aTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGc      >  2:2226170/1‑91 (MQ=255)
                                                                                        tGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCa     >  2:2972587/1‑91 (MQ=255)
                                                                                         gCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCaa    >  1:3650860/1‑91 (MQ=255)
                                                                                         gCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCaa    >  1:1702751/1‑91 (MQ=255)
                                                                                           |                                                                                          
GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAGC  >  USA300TCH1516_ALE/593343‑593523

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