Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,119 C→A 7.8% G355C (GGT→TGT)  patA_2 ← Putative N‑acetyl‑LL‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,055,1190CA 7.8% 439.3 / 22.9 128G355C (GGT→TGT) patA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (54/64);  new base A (2/8);  total (56/72)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 132 total aligned reads displayed.

ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATA  >  USA300TCH1516_ALE/1055029‑1055207
                                                                                          |                                                                                        
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAGACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                          >  1:3293247/1‑91 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                          >  1:526143/1‑91 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                          <  1:2242177/91‑1 (MQ=255)
atttaGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTAcc                                                                                          >  1:1562519/1‑91 (MQ=255)
   taGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGaa                                                                                       <  1:2552424/91‑1 (MQ=255)
   taGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGaa                                                                                       <  1:2552399/91‑1 (MQ=255)
    tatttttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAAt                                                                                      >  1:1928449/5‑91 (MQ=255)
    tatttttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAAt                                                                                      >  1:3702730/5‑91 (MQ=255)
    aGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAAt                                                                                      <  1:457801/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTc                                                                                    <  1:311339/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTc                                                                                    <  1:2677348/91‑1 (MQ=255)
      tATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTc                                                                                    <  1:2677297/91‑1 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                <  1:1073899/91‑1 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                <  1:3303354/91‑1 (MQ=255)
          tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTa                                                                                <  1:1836094/91‑1 (MQ=255)
            aaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaaa                                                                              >  1:2724008/1‑91 (MQ=255)
            aaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaaa                                                                              >  1:2723724/1‑91 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                             <  1:2401980/91‑1 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                             <  1:2961462/91‑1 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                             <  1:2961510/91‑1 (MQ=255)
             atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACAGCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTAAAg                                                                             <  1:644235/91‑1 (MQ=255)
             atatTTAGCAAGTCGCTTCAGACCTTCTTATAAATCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTAAAg                                                                             <  1:1209/91‑1 (MQ=255)
              tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                            <  1:2006346/91‑1 (MQ=255)
              tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                            >  1:2701149/1‑91 (MQ=255)
              tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTAAAgg                                                                            <  1:818098/91‑1 (MQ=255)
               atTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGa                                                                           <  1:372461/91‑1 (MQ=255)
                tttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGAt                                                                          >  1:2168245/1‑91 (MQ=255)
                tttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGAt                                                                          >  1:2168267/1‑91 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                         <  1:2970802/91‑1 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                         <  1:2971378/91‑1 (MQ=255)
                 ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                         >  1:1935635/1‑91 (MQ=255)
                  tAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGa                                                                        >  1:107030/1‑91 (MQ=255)
                    gCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGc                                                                      <  1:3066846/91‑1 (MQ=255)
                    gCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGc                                                                      >  1:844437/1‑91 (MQ=255)
                     cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                     <  1:28182/91‑1 (MQ=255)
                      aaGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                    <  1:2294575/91‑1 (MQ=255)
                      aaGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                    <  1:2372777/91‑1 (MQ=255)
                        gTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTgg                                                                  >  1:2946407/1‑91 (MQ=255)
                         tctcTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGa                                                                 >  1:2779675/1‑91 (MQ=255)
                            ctTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACg                                                              <  1:2157667/91‑1 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                             >  1:1620282/1‑91 (MQ=255)
                             ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                             >  1:1621716/1‑91 (MQ=255)
                                 tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                         >  1:398915/1‑91 (MQ=255)
                                   ccTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTa                                                       <  1:4044962/91‑1 (MQ=255)
                                    cttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTaa                                                      <  1:1933673/91‑1 (MQ=255)
                                     ttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCCGTAAAGGATGAGCCTGGAACGATTGCTaaa                                                     <  1:2373951/91‑1 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                   <  1:402021/91‑1 (MQ=255)
                                       cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                   <  1:895575/91‑1 (MQ=255)
                                           tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                               <  1:3248815/91‑1 (MQ=255)
                                             aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                             >  1:2696557/1‑91 (MQ=255)
                                             aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                             >  1:593569/1‑91 (MQ=255)
                                               aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                           <  1:2170961/91‑1 (MQ=255)
                                                cATCCATTTCATAAGCATAGGAAATCCTNACAANTCCTTTACCGANTTCTGTAAAGGATGNGCCTGGAACGATTGCTAAATGTGTTGATTc                                          <  1:1569963/91‑1 (MQ=255)
                                                  tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                        >  1:2991705/1‑91 (MQ=255)
                                                       ttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTaaa                                   >  1:2117672/1‑91 (MQ=255)
                                                        tCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAAt                                  >  1:1005061/1‑91 (MQ=255)
                                                        tCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAAt                                  <  1:2429857/91‑1 (MQ=255)
                                                           tAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGa                               <  1:370777/91‑1 (MQ=255)
                                                            aaGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGac                              <  1:2420598/91‑1 (MQ=255)
                                                            aaGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGac                              <  1:1065620/91‑1 (MQ=255)
                                                             aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                             >  1:569941/1‑91 (MQ=255)
                                                             aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                             >  1:2332654/1‑91 (MQ=255)
                                                              gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                            >  1:1995743/1‑91 (MQ=255)
                                                              gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                            >  1:3423406/1‑91 (MQ=255)
                                                                aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                          >  1:3880341/1‑91 (MQ=255)
                                                                aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                          <  1:540514/91‑1 (MQ=255)
                                                                 tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                         >  1:668972/1‑91 (MQ=255)
                                                                  aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                        >  1:1034862/1‑91 (MQ=255)
                                                                  aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                        <  1:3848853/91‑1 (MQ=255)
                                                                  aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                        <  1:3649795/91‑1 (MQ=255)
                                                                   ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                       >  1:2218783/1‑91 (MQ=255)
                                                                   ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                       >  1:2317721/1‑91 (MQ=255)
                                                                     aaaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTc                     >  1:1608149/1‑91 (MQ=255)
                                                                      aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                    >  1:623435/1‑91 (MQ=255)
                                                                       aTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGa                   >  1:24482/1‑91 (MQ=255)
                                                                        tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                  >  1:1258146/1‑91 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                 >  1:2804419/1‑91 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCTGGAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                 <  1:2847246/91‑1 (MQ=255)
                                                                         ccTTACAAATCCTTTACCGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                 <  1:1535249/91‑1 (MQ=255)
                                                                           ttACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATc               >  1:1125758/1‑91 (MQ=255)
                                                                           ttACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATc               >  1:2850339/1‑91 (MQ=255)
                                                                              cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTc            <  1:3710571/91‑1 (MQ=255)
                                                                               aaaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCa           <  1:3580864/91‑1 (MQ=255)
                                                                               aaaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCa           <  1:127217/91‑1 (MQ=255)
                                                                                aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt          >  1:285973/1‑91 (MQ=255)
                                                                                  tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCg        >  1:2719848/1‑91 (MQ=255)
                                                                                  tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCg        >  1:2066981/1‑91 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg       >  1:3450620/1‑91 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg       <  1:3126530/91‑1 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg       <  1:1046718/91‑1 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg       <  1:2250061/91‑1 (MQ=255)
                                                                                   ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg       <  1:3949986/91‑1 (MQ=255)
                                                                                    cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt      <  1:4109588/91‑1 (MQ=255)
                                                                                    cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt      >  1:2836405/1‑91 (MQ=255)
                                                                                     tttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg     >  1:573727/1‑91 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGTACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTgg    >  1:3077852/1‑90 (MQ=255)
                                                                                      ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa    >  1:3077886/1‑91 (MQ=255)
                                                                                      ttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa    >  1:2050794/1‑91 (MQ=255)
                                                                                        aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata  >  1:4069698/1‑91 (MQ=255)
                                                                                        aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata  <  1:3118666/91‑1 (MQ=255)
                                                                                          |                                                                                        
ATTTAGTATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATA  >  USA300TCH1516_ALE/1055029‑1055207

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: