Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,126 T→G 9.7% 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,1260TG 9.7% 854.2 / 78.6 278T352T (ACA→ACCpatA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (145/105);  new base G (22/5);  total (168/110)
Fisher's exact test for biased strand distribution p-value = 2.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 282 total aligned reads displayed.

ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAATA  >  USA300TCH1516_ALE/1055036‑1055216
                                                                                          |                                                                                          
aTTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCt                                                                                            <  1:345216/91‑1 (MQ=255)
  ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGt                                                                                          >  1:69750/1‑91 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTa                                                                                         <  1:1877658/91‑1 (MQ=255)
    aaaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaa                                                                                        >  1:403986/1‑91 (MQ=255)
    aaaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaa                                                                                        >  1:1561192/1‑91 (MQ=255)
      atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTAAAg                                                                                      <  2:3659143/91‑1 (MQ=255)
      atagTTAGCGAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGTAAAg                                                                                      <  1:3549810/91‑1 (MQ=255)
       tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                                     <  2:3405386/91‑1 (MQ=255)
       tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                                     <  1:659699/91‑1 (MQ=255)
       tatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                                     <  2:680442/91‑1 (MQ=255)
        atTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGa                                                                                    <  2:2706706/91‑1 (MQ=255)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAAGCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                  <  2:2230883/91‑1 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAGACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              <  1:3694254/91‑1 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              >  1:1066409/1‑91 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              >  1:1067187/1‑91 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                              <  1:3811955/91‑1 (MQ=255)
               aaGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCt                                                                             >  1:322780/1‑91 (MQ=255)
                aGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTg                                                                            <  1:1669010/91‑1 (MQ=255)
                 gTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTgg                                                                           <  1:2012965/91‑1 (MQ=255)
                    tctTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAAc                                                                        <  2:2564081/91‑1 (MQ=255)
                      ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                                      >  2:1962996/1‑91 (MQ=255)
                        cATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAtt                                                                    <  2:159052/91‑1 (MQ=255)
                         aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACACATCCTATACCGAATTCTGTAAAGGATGAGCCTGGAACGAttt                                                                   >  2:3915384/1‑90 (MQ=255)
                         aTACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTg                                                                   <  2:474493/91‑1 (MQ=255)
                          tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATGCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                                  <  2:31013/91‑1 (MQ=255)
                          tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                                  <  1:1290220/91‑1 (MQ=255)
                          tACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGc                                                                  <  2:2389802/91‑1 (MQ=255)
                          gaCCTTCTTTTAACACGTCCATGTCAGGAGCAGAGGAAATCCTTACAAATCCTTTTCCGAATTCTGTAAAGGAGGAGCCTGGAACGATTGc                                                                  <  2:1994244/90‑1 (MQ=255)
                           aCCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCt                                                                 <  1:2663593/91‑1 (MQ=255)
                            ccTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTa                                                                >  2:2804016/1‑91 (MQ=255)
                               tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                             >  1:2775685/1‑91 (MQ=255)
                               tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                             <  2:2876939/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                            <  1:2607168/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                            <  1:1431175/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                            <  2:1845139/91‑1 (MQ=255)
                                 ttttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGTATGAGCCTGGACCGATTGCTAAAtgt                                                           >  2:786041/1‑91 (MQ=255)
                                  tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                          <  2:3747593/91‑1 (MQ=255)
                                  tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCAAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                          <  2:125022/91‑1 (MQ=255)
                                   ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                         <  2:3805458/91‑1 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGTCTGGAACGATTACTAAATGTGtt                                                        >  2:3011386/1‑91 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                        >  1:306495/1‑91 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                        <  2:1657749/91‑1 (MQ=255)
                                      aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                                      >  1:3519242/1‑91 (MQ=255)
                                      aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATCAGCCTGGAATGATTGCTAAATGTGTTGa                                                      >  2:195020/1‑91 (MQ=255)
                                       aaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAt                                                     >  2:1687481/1‑91 (MQ=255)
                                        aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                                    >  2:1143606/1‑91 (MQ=255)
                                        aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                                    >  2:3589588/1‑91 (MQ=255)
                                        aCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAtt                                                    <  1:2462095/91‑1 (MQ=255)
                                         cgtCCATTTCGTAAGCATAGGAAAGCGTTACAAATCGTTTACCGAATTCTGGAAAGGATGAGCCTGGAACGATTGCTAAATGTGTCGAttt                                                   <  2:1893339/89‑2 (MQ=255)
                                         cATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTc                                                   <  2:820202/91‑1 (MQ=255)
                                           tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                                 >  1:151203/1‑91 (MQ=255)
                                           tCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCaa                                                 >  1:3302954/1‑91 (MQ=255)
                                            ccATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAg                                                >  1:1534453/1‑91 (MQ=255)
                                             cATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGt                                               >  1:1273019/1‑91 (MQ=255)
                                             cATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGGGTTGATTCAAGt                                               <  1:3056030/91‑1 (MQ=255)
                                              aTTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTa                                              >  1:794895/1‑91 (MQ=255)
                                                   aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                         >  2:1209798/1‑91 (MQ=255)
                                                   aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                         >  2:1536497/1‑91 (MQ=255)
                                                   aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                         >  2:3304029/1‑91 (MQ=255)
                                                   aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGCAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                         >  1:3426191/1‑91 (MQ=255)
                                                     aaGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGac                                       >  2:2724645/1‑91 (MQ=255)
                                                      aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                                      <  1:2370440/91‑1 (MQ=255)
                                                       gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                                     >  1:2063230/1‑91 (MQ=255)
                                                          tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                                  <  2:3352383/91‑1 (MQ=255)
                                                          tAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaa                                  >  2:1353378/1‑91 (MQ=255)
                                                           aGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaaaa                                 <  2:2789259/91‑1 (MQ=255)
                                                            ggAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAt                                <  2:3204736/91‑1 (MQ=255)
                                                             gAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAAtt                               >  1:855651/1‑91 (MQ=255)
                                                                 tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                           >  1:3633596/1‑91 (MQ=255)
                                                                  ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaca                          >  2:1488373/1‑89 (MQ=255)
                                                                  ccTTACAAAGCCTTTACCGAGTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                          <  1:1879259/91‑1 (MQ=255)
                                                                   cTTACAAATCCTTTACCGTATTCTGTAAAGGATGAGCCTGGAGCGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAt                         <  2:26117/91‑1 (MQ=255)
                                                                    ttACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATc                        >  2:2519169/1‑91 (MQ=255)
                                                                    ttACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATc                        <  2:3826858/91‑1 (MQ=255)
                                                                     tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg                       >  1:1627159/1‑91 (MQ=255)
                                                                      aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGt                      >  1:3845023/1‑91 (MQ=255)
                                                                      aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGATAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAtaga                      >  1:603007/1‑88 (MQ=255)
                                                                       cAAGTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAGATGTGTTGATTCAAGTAAATCGACACAAATTTCGACATCGTc                     <  2:2630449/91‑1 (MQ=255)
                                                                       cAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAGTTCGAAATCGTc                     <  1:3041354/91‑1 (MQ=255)
                                                                         aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt                   <  2:1171985/91‑1 (MQ=255)
                                                                         aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt                   <  2:6180/91‑1 (MQ=255)
                                                                          aTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAACGACATc                  >  1:490617/1‑91 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                <  2:1383229/91‑1 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg                >  1:3902716/1‑91 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               <  1:2871953/91‑1 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               >  1:967010/1‑91 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               >  1:552274/1‑91 (MQ=255)
                                                                             cTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGt               >  1:2682520/1‑91 (MQ=255)
                                                                               ttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa             >  1:1770783/1‑91 (MQ=255)
                                                                                tacagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACACAATTCGAAATCGTCATCGGTGat            >  2:3883935/5‑91 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGACATCGTCATCGGTGata           >  2:3780868/4‑91 (MQ=255)
                                                                                 aCCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata           >  1:1808603/1‑91 (MQ=255)
                                                                                  cagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGatat          >  1:768862/3‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg         >  1:2116891/2‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg         >  1:1732945/2‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg         >  2:2307562/2‑91 (MQ=255)
                                                                                    gAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt        >  1:2228736/1‑91 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt      >  2:1717634/1‑91 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAACCCGAAATCGTGATCGGTGATATGttt      >  2:3421623/1‑91 (MQ=255)
                                                                                          tGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAATa  >  2:1104444/1‑91 (MQ=255)
                                                                                          |                                                                                          
ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAATA  >  USA300TCH1516_ALE/1055036‑1055216

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: