Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 979,347 T→C 5.2% intergenic (+17/+42) clpB_1 → / ← ampR Chaperone protein ClpB/HTH‑type transcriptional activator AmpR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE979,3470TC 5.2% 690.0 / 17.8 211intergenic (+17/+42)clpB_1/ampRChaperone protein ClpB/HTH‑type transcriptional activator AmpR
Reads supporting (aligned to +/- strand):  ref base T (92/108);  new base C (7/4);  total (99/112)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 234 total aligned reads displayed.

AGGTACAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACC  >  USA300TCH1516_ALE/979258‑979437
                                                                                         |                                                                                          
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  2:1885925/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  1:2262245/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       <  1:3642838/91‑1 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  1:1933327/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  2:1379135/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  1:406013/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  2:3731474/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGtgatga                                                                                       >  2:1122230/1‑91 (MQ=255)
    aCAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGcgatta                                                                                       >  1:544746/1‑89 (MQ=255)
     cAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTgatgat                                                                                      >  2:3553281/1‑91 (MQ=255)
     cAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTgatgat                                                                                      >  1:2159918/1‑91 (MQ=255)
             aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATcatc                                                                              <  2:2159918/91‑4 (MQ=255)
             aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCa                                                                              >  2:1202822/1‑91 (MQ=255)
             aaaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCa                                                                              >  2:3879053/1‑91 (MQ=255)
              aaGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCAt                                                                             >  1:261629/1‑91 (MQ=255)
               aGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCAtt                                                                            <  2:298200/91‑1 (MQ=255)
                gTTACTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                           <  1:3237513/91‑1 (MQ=255)
                gTTAATGTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                           <  2:349936/91‑1 (MQ=255)
                  tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                         >  1:3457542/1‑91 (MQ=255)
                  tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCtt                                                                         >  1:3458093/1‑91 (MQ=255)
                  tAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGCGATTAATCATCACATTCtt                                                                         <  1:2404622/91‑1 (MQ=25)
                  tAATTTAAATTCAGACAATAAATTAACTATTAATAATAAAAAAATTCATGAATAAAATTGAACCAAGAATGCGATTAATCATCACATTCtt                                                                         <  2:468679/91‑1 (MQ=16)
                   aaTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTg                                                                        >  1:2872509/1‑91 (MQ=255)
                    aTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTgg                                                                       >  2:2650894/1‑91 (MQ=255)
                    aTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTgg                                                                       >  2:132557/1‑91 (MQ=255)
                    aTTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTgg                                                                       <  1:1445487/91‑1 (MQ=255)
                       tAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                           >  2:2851442/1‑84 (MQ=255)
                       tAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTc                                                                           <  1:2851442/84‑1 (MQ=255)
                       tAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGttt                                                                    <  2:353541/91‑1 (MQ=255)
                             cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                              <  1:2572430/91‑1 (MQ=255)
                             cAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGa                                                              <  1:769046/91‑1 (MQ=255)
                               gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAtt                                                            <  2:1691339/91‑1 (MQ=255)
                               gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAtt                                                            >  1:955146/1‑91 (MQ=255)
                               gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAtt                                                            >  2:2023466/1‑91 (MQ=255)
                               gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAtt                                                            <  2:2885987/91‑1 (MQ=255)
                               gACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATCCTTGGTTTTTTAGAtt                                                            >  2:816937/1‑91 (MQ=255)
                                aCAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAttt                                                           >  2:795071/1‑91 (MQ=255)
                                aCAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGAttt                                                           <  2:1815759/91‑1 (MQ=255)
                                 cAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTg                                                          >  1:3027695/1‑91 (MQ=255)
                                    tAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGaaa                                                       >  2:1983815/1‑91 (MQ=255)
                                    tAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGaaa                                                       >  1:1552110/1‑91 (MQ=255)
                                    tAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGaaa                                                       >  2:2273331/1‑91 (MQ=255)
                                      aCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAtt                                                     >  1:487525/1‑91 (MQ=255)
                                      aCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAtt                                                     >  1:1301944/1‑91 (MQ=255)
                                      aCTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAtt                                                     <  2:211840/91‑1 (MQ=255)
                                       cTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAttt                                                    >  1:1606189/1‑91 (MQ=255)
                                       cTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAAttt                                                    >  2:1087373/1‑91 (MQ=255)
                                        ttAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATtta                                                   >  2:26214/1‑91 (MQ=255)
                                        ttAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATtta                                                   <  2:3019249/91‑1 (MQ=255)
                                            cGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTtattat                                               >  2:1423210/1‑91 (MQ=255)
                                             gTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATtata                                              <  2:2392760/91‑1 (MQ=255)
                                              tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                             <  2:2454446/91‑1 (MQ=255)
                                              tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                             <  2:3203790/91‑1 (MQ=255)
                                              tttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATaa                                             <  1:1113596/91‑1 (MQ=255)
                                               ttAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAAt                                            <  1:1853590/91‑1 (MQ=255)
                                                tAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAAtt                                           <  1:1202131/91‑1 (MQ=255)
                                                  atgatgAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATGTATTATAATTGt                                         <  1:3000004/91‑1 (MQ=255)
                                                  aTGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGt                                         <  1:1870890/91‑1 (MQ=255)
                                                  aTGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGt                                         <  2:578542/91‑1 (MQ=255)
                                                   tGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGtt                                        >  2:3818804/1‑91 (MQ=255)
                                                   tGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGtt                                        <  2:2521921/91‑1 (MQ=255)
                                                   tGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGtt                                        >  2:3819062/1‑91 (MQ=255)
                                                    gTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTc                                       <  2:3389346/91‑1 (MQ=255)
                                                     ttGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCa                                      >  1:3247916/1‑91 (MQ=255)
                                                      tGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCAt                                     >  1:1650954/1‑91 (MQ=255)
                                                       gAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCAtt                                    <  2:2127858/91‑1 (MQ=255)
                                                        aaaaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTg                                   >  1:2066779/1‑91 (MQ=255)
                                                         aaaaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTgg                                  <  1:729804/91‑1 (MQ=255)
                                                           aaaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCa                                <  1:2865668/91‑1 (MQ=255)
                                                            aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                               <  1:2440058/91‑1 (MQ=255)
                                                            aaTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATCCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCaa                               <  1:3814546/91‑1 (MQ=255)
                                                             aTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAat                              <  1:1358994/91‑1 (MQ=255)
                                                             aTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAat                              <  1:1894287/91‑1 (MQ=255)
                                                             aTTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAat                              <  1:2023466/91‑1 (MQ=255)
                                                              ttCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAata                             >  1:208880/1‑91 (MQ=255)
                                                               tCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAatat                            <  1:535546/91‑1 (MQ=255)
                                                               tCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAatat                            >  2:2347508/1‑91 (MQ=255)
                                                               tCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAatat                            <  1:2485544/91‑1 (MQ=255)
                                                                cATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAAtata                           <  1:2295933/91‑1 (MQ=255)
                                                                cATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAAtata                           <  1:3553890/91‑1 (MQ=255)
                                                                  tgcaTAAAATTGAACCAAGCATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATCGGCAATCTACt                         >  2:3630774/4‑91 (MQ=255)
                                                                  tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                         <  1:3188511/91‑1 (MQ=255)
                                                                  tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                         >  1:659673/1‑91 (MQ=255)
                                                                  tGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACt                         >  1:2918468/1‑91 (MQ=255)
                                                                   gAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtt                        >  1:2006055/1‑91 (MQ=255)
                                                                     aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                      <  1:889553/91‑1 (MQ=255)
                                                                     aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                      <  1:170902/91‑1 (MQ=255)
                                                                     aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                      <  1:2775404/91‑1 (MQ=255)
                                                                     aTAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACtttt                      <  1:3699943/91‑1 (MQ=255)
                                                                          aTTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAAt                 <  1:2240657/91‑1 (MQ=255)
                                                                            tGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGa               <  2:33318/91‑1 (MQ=255)
                                                                             gAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAt              >  1:3333419/1‑91 (MQ=255)
                                                                               aCCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt            >  1:935765/1‑91 (MQ=255)
                                                                               aCCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt            <  1:1047477/91‑1 (MQ=255)
                                                                               aCCAAGAATGCGATTAATCATCACATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGAttt            >  2:1910702/1‑91 (MQ=25)
                                                                                 cAAGAATGCGATTAATCATCACATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaa          >  1:1251606/1‑91 (MQ=25)
                                                                                   aGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaata        <  2:3301156/91‑1 (MQ=255)
                                                                                    gAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTtaataa       >  1:1495974/1‑91 (MQ=255)
                                                                                     aaTGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTAATATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATaaa      >  2:3477801/1‑91 (MQ=255)
                                                                                       tgtgATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAaca    <  2:2754126/91‑1 (MQ=255)
                                                                                         |                                                                                          
AGGTACAACGATTAAAGTTAATTTAAATTCAGACAATAACTTAACGTTTAATGTTGAAAAAATTCATGAATAAAATTGAACCAAGAATGTGATGATTAATCGCATTCTTGGTTTTTTAGATTTGAAATTTATTATAATTGTTCATTGGCAATATACTTTTTAAATGATTTAATAAACACC  >  USA300TCH1516_ALE/979258‑979437

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: