Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 593,433 T→A 10.4% 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,4330TA10.4% 441.1 / 35.5 135intergenic (+25/‑112)rpoC/rplGBDNA‑directed RNA polymerase subunit beta'/Ribosome‑associated protein L7Ae‑like protein
Reads supporting (aligned to +/- strand):  ref base T (68/53);  new base A (6/8);  total (74/61)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 136 total aligned reads displayed.

GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAG  >  USA300TCH1516_ALE/593343‑593522
                                                                                           |                                                                                         
gTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCt                                                                                           <  1:3132164/91‑1 (MQ=255)
    aaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAg                                                                                       >  1:1289177/1‑91 (MQ=255)
    aaTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAg                                                                                       >  1:1289162/1‑91 (MQ=255)
     aTACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGc                                                                                      >  1:733756/1‑91 (MQ=255)
      tACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATGTAACAAGAGGCTAAAGCATTAGc                                                                                      <  1:482209/91‑1 (MQ=255)
       aCGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCct                                                                                    >  1:1120780/1‑91 (MQ=255)
        cGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctc                                                                                   >  1:1696921/1‑91 (MQ=255)
         gAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctct                                                                                  >  1:2695530/1‑91 (MQ=255)
         gAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCctct                                                                                  <  1:3848078/91‑1 (MQ=255)
          aaaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCtt                                                                                 <  1:755273/91‑1 (MQ=255)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctg                                                                                 <  1:242972/91‑2 (MQ=38)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctg                                                                                 <  1:1943222/91‑2 (MQ=38)
           aaaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTg                                                                                >  1:2210919/1‑91 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                               <  1:2575485/91‑1 (MQ=255)
            aaaaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGt                                                                               <  1:1807420/91‑1 (MQ=255)
              aaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAAGAGGCTAAAGCATTAGCctctgtta                                                                              <  1:1615148/91‑5 (MQ=37)
              aaCAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTa                                                                             >  1:1571845/1‑91 (MQ=255)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGATTAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                         <  1:3284141/91‑1 (MQ=255)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGACTAACAAGTAAATAACAAGAGGCTAAAGCATTAGCCTCTGTTAtata                                                                         <  1:3300752/91‑4 (MQ=18)
                  gCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAtttt                                                                         <  1:3371187/91‑1 (MQ=255)
                   cTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTAttttt                                                                        >  1:1100556/1‑91 (MQ=255)
                       aCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGt                                                                    <  1:3607486/91‑1 (MQ=255)
                       aCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGt                                                                    <  1:471594/91‑1 (MQ=255)
                       aCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGt                                                                    <  1:471512/91‑1 (MQ=255)
                       aCCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGt                                                                    <  1:1503708/91‑1 (MQ=255)
                         cAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTaa                                                                  <  1:4106345/91‑1 (MQ=255)
                          aGTTGCAGAAGTTGAATCTCAAACGGAAGTAGCGGAATAACAAGTGTATAACAGGGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTaaa                                                                 <  1:659724/91‑1 (MQ=255)
                             tGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatta                                                              >  1:1745287/1‑91 (MQ=255)
                             tGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatta                                                              >  1:1746053/1‑91 (MQ=255)
                             tGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAatta                                                              >  1:335948/1‑91 (MQ=255)
                              gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                             >  1:3536125/1‑91 (MQ=255)
                              gCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAattat                                                             >  1:879655/1‑91 (MQ=255)
                                actaGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTAttt                                                           >  1:2583738/4‑91 (MQ=255)
                                 gAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTg                                                          >  1:1067722/1‑91 (MQ=255)
                                  aaGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGACGCTAATGCTTTAGCCTCTTGTTATTTTTACGTAAATTATTTGa                                                         >  1:4144285/1‑91 (MQ=255)
                                    gTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAtt                                                       <  1:3386035/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:2182813/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:223959/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:125328/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:2182812/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:865035/91‑1 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      >  1:974707/1‑91 (MQ=255)
                                     ttGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGAttt                                                      <  1:1259007/91‑1 (MQ=255)
                                         aTCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATg                                                  >  1:1366271/1‑91 (MQ=255)
                                           ctcAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGtt                                                >  1:3802977/1‑91 (MQ=255)
                                           ctcAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGtt                                                >  1:3705044/1‑91 (MQ=255)
                                            tcAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTg                                               <  1:1136239/91‑1 (MQ=255)
                                             cAAACTGAAGTAACGGGATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                              <  1:421925/91‑1 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                              <  1:420854/91‑1 (MQ=255)
                                             cAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGa                                              <  1:847043/91‑1 (MQ=255)
                                               aaCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACg                                            <  1:3409545/91‑1 (MQ=255)
                                               aaCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACg                                            <  1:3594133/91‑1 (MQ=255)
                                                aCTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGa                                           >  1:1705290/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATtctc                                     >  1:692059/1‑91 (MQ=255)
                                                      gTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTCGATTTAATGTTGACGAATtctc                                     >  1:1396716/1‑91 (MQ=255)
                                                       tAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTctct                                    >  1:1697572/1‑91 (MQ=255)
                                                       tAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTctct                                    >  1:1697555/1‑91 (MQ=255)
                                                        aaCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCtt                                   <  1:2427710/91‑1 (MQ=255)
                                                         aCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTg                                  >  1:3317543/1‑91 (MQ=255)
                                                         aCGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTg                                  >  1:3088047/1‑91 (MQ=255)
                                                          cGGAATAGCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGt                                 <  1:2064299/91‑1 (MQ=255)
                                                          cGGAATAGCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGt                                 <  1:2064298/91‑1 (MQ=255)
                                                               tAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAAt                            >  1:3458666/1‑91 (MQ=255)
                                                               tAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAAt                            >  1:2026251/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                           <  1:54562/91‑1 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                           <  1:54552/91‑1 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                           >  1:1219988/1‑91 (MQ=255)
                                                                aaCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATg                           >  1:1531869/1‑91 (MQ=255)
                                                                 aCAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGt                          >  1:3099240/1‑91 (MQ=255)
                                                                  cAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGtt                         <  1:181232/91‑1 (MQ=255)
                                                                     gTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAat                      >  1:3787221/1‑91 (MQ=255)
                                                                      tatataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAata                     >  1:1062646/1‑91 (MQ=255)
                                                                      tatataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAata                     >  1:1062645/1‑91 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatat                    >  1:3709148/1‑91 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatat                    <  1:2951769/91‑1 (MQ=255)
                                                                       atataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatat                    >  1:614144/1‑91 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                  <  1:372944/91‑1 (MQ=255)
                                                                         ataACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAatatat                  >  1:1702548/1‑91 (MQ=255)
                                                                          taACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATAtt                 >  1:925354/1‑91 (MQ=255)
                                                                           aaCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTa                >  1:4034158/1‑91 (MQ=255)
                                                                            acaagagGCTAAAGCATTAGCCTCTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa               >  1:4057146/4‑91 (MQ=37)
                                                                            acaagagGCTAAAGCATTAGCCTCTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa               >  1:950536/4‑91 (MQ=37)
                                                                            acaagagGCTAAAGCATTAGCCTCTGTTATTTTTATGTAAATTATTTGATTTAATGTCGACGAATTCTCTTGTTCAATGTTAATATATTaa               >  1:3401462/4‑91 (MQ=37)
                                                                            aCAGAGGCTAATGCTTTAGTCTCTTGTTCTTTTTATGTAAATTACTTGACTTAATGTTGCCGAATTCTCTTGTTCAATGTTAATATATTaa               >  1:1200762/1‑91 (MQ=255)
                                                                            aCAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaa               <  1:2733665/91‑1 (MQ=255)
                                                                             cAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTaaa              <  1:988574/91‑1 (MQ=255)
                                                                                gagGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAgg            <  1:1444905/91‑1 (MQ=255)
                                                                                gagGCTAAAGCATTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAgg            >  1:2168516/1‑91 (MQ=255)
                                                                                 agGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGt           >  1:3278765/1‑91 (MQ=255)
                                                                                 agGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGt           <  1:87921/91‑1 (MQ=255)
                                                                                 agGCAAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGt           >  1:3278048/1‑91 (MQ=255)
                                                                                  ggCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGtt          <  1:396263/91‑1 (MQ=255)
                                                                                   gCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGCTTg         >  1:3343065/1‑91 (MQ=255)
                                                                                    cTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGa        <  1:2535498/91‑1 (MQ=255)
                                                                                    cTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGa        >  1:891289/1‑91 (MQ=255)
                                                                                     tAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGAt       >  1:3456925/1‑91 (MQ=255)
                                                                                     tAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGAt       >  1:1254553/1‑91 (MQ=255)
                                                                                       aTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGc     >  1:3228664/1‑91 (MQ=255)
                                                                                       aTGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGc     >  1:1262856/1‑91 (MQ=255)
                                                                                          cTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAg  >  1:3976693/1‑91 (MQ=255)
                                                                                           |                                                                                         
GTAAAATACGAAAAAACAGCTAAACCAGTTGCAGAAGTTGAATCTCAAACTGAAGTAACGGAATAACAAGTATATAACAGAGGCTAATGCTTTAGCCTCTTGTTATTTTTATGTAAATTATTTGATTTAATGTTGACGAATTCTCTTGTTCAATGTTAATATATTAAAGGTTGATGCAAG  >  USA300TCH1516_ALE/593343‑593522

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: