Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,126 T→G 13.2% T352T (ACA→ACC 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,1260TG13.2% 397.9 / 48.2 132T352T (ACA→ACCpatA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (56/57);  new base G (12/6);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAATA  >  USA300TCH1516_ALE/1055036‑1055216
                                                                                          |                                                                                          
aTTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCt                                                                                            <  1:745677/91‑1 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                         <  1:3303354/91‑1 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                         <  1:1073899/91‑1 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTa                                                                                         <  1:1836094/91‑1 (MQ=255)
     aaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaaa                                                                                       >  1:2724008/1‑91 (MQ=255)
      atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                                      <  1:2961462/91‑1 (MQ=255)
      atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                                      <  1:2401980/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)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                  <  1:2970802/91‑1 (MQ=255)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                  <  1:2971378/91‑1 (MQ=255)
           tAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGa                                                                                 >  1:107030/1‑91 (MQ=255)
             gCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGc                                                                               <  1:3066846/91‑1 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              <  1:28182/91‑1 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              <  1:3916689/91‑1 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              <  1:1127531/91‑1 (MQ=255)
               aaGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                             <  1:2294575/91‑1 (MQ=255)
                aGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTg                                                                            <  1:305865/91‑1 (MQ=255)
                aGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTg                                                                            >  1:2140173/1‑91 (MQ=255)
                 gTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTgg                                                                           >  1:2946407/1‑91 (MQ=255)
                  tctcTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGa                                                                          >  1:3352467/1‑91 (MQ=255)
                     ctTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACg                                                                       <  1:2157667/91‑1 (MQ=255)
                      ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                                      >  1:1620282/1‑91 (MQ=255)
                         aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTg                                                                   >  1:1591743/1‑91 (MQ=255)
                          tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                                  <  1:454934/91‑1 (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:895575/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                            <  1:3167426/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                            <  1:402021/91‑1 (MQ=255)
                                  tttAAAACATCCATTTCATACGCATAGGAAATCCTCACAAATCCGTTACCGAATTCTGTAGAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                          <  1:2338350/91‑1 (MQ=255)
                                      aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                                      >  1:593569/1‑91 (MQ=255)
                                        aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                                    <  1:2170961/91‑1 (MQ=255)
                                        aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                                    <  1:1668479/91‑1 (MQ=255)
                                          aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                                  >  1:2418734/1‑91 (MQ=255)
                                           tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                                 <  1:2282406/91‑1 (MQ=255)
                                           tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                                 >  1:2991705/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:1065727/91‑1 (MQ=255)
                                                      aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                                      >  1:2332654/1‑91 (MQ=255)
                                                      aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                                      >  1:569941/1‑91 (MQ=255)
                                                       gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                                     >  1:1995743/1‑91 (MQ=255)
                                                         aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                                   >  1:3880341/1‑91 (MQ=255)
                                                         aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                                   <  1:540514/91‑1 (MQ=255)
                                                           aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                                 <  1:3848853/91‑1 (MQ=255)
                                                           aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                                 >  1:1034862/1‑91 (MQ=255)
                                                            ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                                >  1:2218783/1‑91 (MQ=255)
                                                            ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                                <  1:2705701/91‑1 (MQ=255)
                                                              aaaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTc                              >  1:1608149/1‑91 (MQ=255)
                                                               aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                             >  1:2767107/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:4034052/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:2850339/1‑91 (MQ=255)
                                                                       cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTc                     <  1:3710571/91‑1 (MQ=255)
                                                                        aaaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCa                    <  1:3580864/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)
                                                                           tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCg                 <  1:3625642/91‑1 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                <  1:1046718/91‑1 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                <  1:3949986/91‑1 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                >  1:3450620/1‑91 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                <  1:3126530/91‑1 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               >  1:2836405/1‑91 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               <  1:4109588/91‑1 (MQ=255)
                                                                              tttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg              >  1:573727/1‑91 (MQ=255)
                                                                               ttACCGAATTCTGTAAAGGATGAGCCTGGTACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTgg             >  1:3077852/1‑90 (MQ=255)
                                                                               ttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa             >  1:2050794/1‑91 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata           >  1:3651624/4‑91 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata           >  1:396847/4‑91 (MQ=255)
                                                                                 aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata           <  1:3118666/91‑1 (MQ=255)
                                                                                    gAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        <  1:1966067/91‑1 (MQ=255)
                                                                                    gAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        >  1:2301704/1‑91 (MQ=255)
                                                                                    gAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        >  1:2330106/1‑91 (MQ=255)
                                                                                    gAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        >  1:84832/1‑91 (MQ=255)
                                                                                    gAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        >  1:260595/1‑91 (MQ=255)
                                                                                     aaTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGtt       <  1:651438/91‑1 (MQ=255)
                                                                                     aaTTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGtt       >  1:1839738/1‑91 (MQ=255)
                                                                                     aaTTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGtt       >  1:2059417/1‑91 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa     >  1:2204335/1‑91 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa     >  1:761556/1‑91 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa     >  1:1925747/1‑91 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa     >  1:1927484/1‑91 (MQ=255)
                                                                                        tCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTaa    <  1:1448238/91‑1 (MQ=255)
                                                                                         cTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAAt   >  1:1354943/1‑91 (MQ=255)
                                                                                          tGTAAAGGATGAGCCTGGAACGATTGCTAAATGCGTTAATTCAAGTCAATCGACACAAACTGCGAAATCGTCATCGGTGATATGTTTAATa  >  1:3316297/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAATA  >  USA300TCH1516_ALE/1055036‑1055216

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: