Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,055,126 T→G 11.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,1260TG11.2% 805.0 / 88.4 258T352T (ACA→ACCpatA_2Putative N‑acetyl‑LL‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (127/101);  new base G (26/3);  total (154/104)
Fisher's exact test for biased strand distribution p-value = 2.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 262 total aligned reads displayed.

ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAA  >  USA300TCH1516_ALE/1055036‑1055214
                                                                                          |                                                                                        
  ttAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACAGAATTCGGt                                                                                        <  2:1503597/91‑1 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                       >  1:2388280/1‑91 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTa                                                                                       >  1:1619903/1‑91 (MQ=255)
   tAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTATGTa                                                                                       >  1:3237658/1‑91 (MQ=255)
     aaTATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTaaa                                                                                     >  2:1680132/1‑91 (MQ=255)
      atatTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAg                                                                                    <  1:1010195/91‑1 (MQ=255)
        atTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGa                                                                                  >  1:1063751/1‑91 (MQ=255)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                <  1:3384104/91‑1 (MQ=255)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                <  1:3246327/91‑1 (MQ=255)
          ttAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATg                                                                                >  2:2730468/1‑91 (MQ=255)
            aGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAg                                                                              >  2:321210/1‑91 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAgg                                                                                   >  2:2130043/1‑84 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                            >  1:1066567/1‑91 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                            >  2:24911/1‑91 (MQ=255)
              cAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGcc                                                                            >  2:24154/1‑91 (MQ=255)
                      ttCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGa                                                                    <  2:1082648/91‑1 (MQ=255)
                        cATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAtt                                                                  >  2:3050303/1‑91 (MQ=255)
                        cATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGAtt                                                                  >  2:529576/1‑91 (MQ=255)
                            ccTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTa                                                              >  2:2964760/1‑91 (MQ=255)
                             cttcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTaa                                                             <  1:2308568/91‑1 (MQ=255)
                               tcttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAt                                                           >  2:2860323/1‑91 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACGAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                          <  1:1349436/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                          <  2:546082/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                          <  2:3359265/91‑1 (MQ=255)
                                cttTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtg                                                          <  1:115395/91‑1 (MQ=255)
                                  tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                        >  1:2166961/1‑91 (MQ=255)
                                  tttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAAtgtg                                                        <  1:3185825/91‑1 (MQ=255)
                                   ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                       <  1:2489549/91‑1 (MQ=255)
                                   ttAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                       <  1:489441/91‑1 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                      <  1:2542074/91‑1 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                      <  1:2699600/91‑1 (MQ=255)
                                    tAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGtt                                                      >  1:732400/1‑91 (MQ=255)
                                      aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                                    <  2:3594846/91‑1 (MQ=255)
                                      aaaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGa                                                    >  1:1753874/1‑91 (MQ=255)
                                       aaCATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGAt                                                   >  2:2373427/1‑91 (MQ=255)
                                         cATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTc                                                 <  1:1543050/91‑1 (MQ=255)
                                          aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                                >  1:851520/1‑91 (MQ=255)
                                          aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                                <  2:13914/91‑1 (MQ=255)
                                          aTCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCa                                                >  2:745925/1‑91 (MQ=255)
                                            ccATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAg                                              >  2:333017/1‑91 (MQ=255)
                                              aTTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTa                                            >  2:3085048/1‑91 (MQ=255)
                                               tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                       <  1:1421152/79‑1 (MQ=255)
                                               tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATgtgt                                                       >  2:1421152/1‑79 (MQ=255)
                                               tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTaa                                           >  1:3414686/1‑91 (MQ=255)
                                               tttCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATCCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTaa                                           >  2:124899/1‑91 (MQ=255)
                                                   aTAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCg                                       >  1:2897882/1‑91 (MQ=255)
                                                      aGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGaca                                    <  1:2121083/91‑1 (MQ=255)
                                                       gCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGacac                                   <  2:2089465/91‑1 (MQ=255)
                                                         aTAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACaa                                 <  2:2326852/91‑1 (MQ=255)
                                                               aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                           >  1:528643/1‑91 (MQ=255)
                                                               aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                           >  2:328987/1‑91 (MQ=255)
                                                               aaTCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCg                           <  2:2778116/91‑1 (MQ=255)
                                                                 tCCTTACGAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTAGaa                         <  2:3437245/91‑1 (MQ=255)
                                                                 tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                         <  1:2865800/91‑1 (MQ=255)
                                                                 tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                         <  1:2361388/91‑1 (MQ=255)
                                                                 tCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaa                         <  1:329050/91‑1 (MQ=255)
                                                                  ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                        <  1:1680132/91‑1 (MQ=255)
                                                                  ccTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGaaa                        <  1:1009551/91‑1 (MQ=255)
                                                                   cTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAAt                       <  1:1079542/91‑1 (MQ=255)
                                                                    ttACAAATCCTTTACCGAATTCTGTAAAGGATGTGCCAGGAACTATTACTCATTGTGTTGAGTCAAGTAAAACGACACAATGTTCGAACTc                      >  2:2600405/1‑91 (MQ=255)
                                                                     tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg                     >  2:1664332/1‑91 (MQ=255)
                                                                     tACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCg                     <  2:437196/91‑1 (MQ=255)
                                                                      aCAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGt                    >  2:3693456/1‑91 (MQ=255)
                                                                        tcaTCCTTTACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCa                  >  1:3427192/3‑91 (MQ=38)
                                                                         aaTCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCAt                 <  1:3853009/91‑1 (MQ=255)
                                                                          aTCCTTTACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATc                >  2:2691160/1‑91 (MQ=255)
                                                                           tCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATGGTCATCg               >  1:1947933/1‑91 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg              >  1:2383286/1‑91 (MQ=255)
                                                                            ccTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCgg              <  1:3334857/91‑1 (MQ=255)
                                                                              tttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg            >  2:2268864/1‑91 (MQ=255)
                                                                              tttACAGAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg            >  2:2011488/1‑91 (MQ=255)
                                                                               ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa           >  1:1507690/1‑91 (MQ=255)
                                                                               ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa           <  1:3685492/91‑1 (MQ=255)
                                                                               ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa           <  1:2674537/91‑1 (MQ=255)
                                                                               ttACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGa           <  2:3801199/91‑1 (MQ=255)
                                                                                tacagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGat          >  1:792782/5‑91 (MQ=255)
                                                                                tacagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGat          >  1:3220025/5‑91 (MQ=255)
                                                                                tacagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGat          >  2:654850/5‑91 (MQ=255)
                                                                                tACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTg            <  1:614501/89‑1 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata         >  2:3695463/4‑91 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata         >  2:2051187/4‑91 (MQ=255)
                                                                                 acagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGata         >  1:3556713/4‑91 (MQ=255)
                                                                                  ccGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGatat        <  2:1575707/91‑1 (MQ=255)
                                                                                  ccGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGatat        <  2:1175344/91‑1 (MQ=255)
                                                                                  cagAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGatat        >  1:3076798/3‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg       >  1:2034658/2‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg       >  2:1086444/2‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg       >  2:1612858/2‑91 (MQ=255)
                                                                                   agAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATg       >  1:2174227/2‑91 (MQ=255)
                                                                                    gAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt      >  1:2329435/1‑91 (MQ=255)
                                                                                    gAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt      >  1:1252129/1‑91 (MQ=255)
                                                                                    gAATTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGt      >  2:764836/1‑91 (MQ=255)
                                                                                     aaTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGtt     <  1:28149/91‑1 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt    <  1:3829167/91‑1 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt    <  2:19736/91‑1 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt    >  1:1214548/1‑91 (MQ=255)
                                                                                      aTTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt    >  1:386500/1‑91 (MQ=255)
                                                                                      aTTCGGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGttt    <  1:654850/91‑1 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa   <  2:3769343/91‑1 (MQ=255)
                                                                                       ttCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTa   <  2:80421/91‑1 (MQ=255)
                                                                                          |                                                                                        
ATTTAAATATTTAGCAAGTCTCTTCATACCTTCTTTTAAAACATCCATTTCATAAGCATAGGAAATCCTTACAAATCCTTTACCGAATTCTGTAAAGGATGAGCCTGGAACGATTGCTAAATGTGTTGATTCAAGTAAATCGACACAAAATTCGAAATCGTCATCGGTGATATGTTTAA  >  USA300TCH1516_ALE/1055036‑1055214

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: