BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	418,825	A→T	7.2%	N53I (AAT→ATT)	rpsF →	30S ribosomal protein S6

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	418,825	0	A	T	7.2%	772.8 / 33.2	228	N53I (AAT→ATT)	rpsF	30S ribosomal protein S6
Reads supporting (aligned to +/- strand): ref base A (110/100); new base T (7/10); total (118/110)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

TTGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAG  >  USA300TCH1516_ALE/418735‑418915

ttGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCaa                                                                                            >  1:952550/1‑91 (MQ=255)

 tGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAAt                                                                                           <  1:2848548/91‑1 (MQ=255)

 tGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAAt                                                                                           <  1:2750967/91‑1 (MQ=255)

   aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGa                                                                                         <  2:2509222/91‑1 (MQ=255)

   aaCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGa                                                                                         <  1:1638081/91‑1 (MQ=255)

    aCGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAt                                                                                        >  1:301365/1‑91 (MQ=255)

     cGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCACTGAtt                                                                                       >  1:272862/1‑91 (MQ=255)

     cGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAtt                                                                                       >  1:4761/1‑91 (MQ=255)

     acgttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGAtt                                                                                       >  2:114442/4‑91 (MQ=255)

       tttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCATTGATTTc                                                                                     <  2:2552571/91‑1 (MQ=255)

        ttCAACGGTATCTTATCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAGACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                    >  1:3377083/1‑91 (MQ=255)

        ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                    >  2:927850/1‑91 (MQ=255)

        ttCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCa                                                                                    >  2:901359/1‑91 (MQ=255)

         tCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCaa                                                                                   <  2:2861169/91‑1 (MQ=255)

         tCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCaa                                                                                   <  2:1922631/91‑1 (MQ=255)

          cAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCaaa                                                                                  >  1:3174222/1‑91 (MQ=255)

          cAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCaaa                                                                                  >  2:3053924/1‑91 (MQ=255)

           aaCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAg                                                                                 >  1:980035/1‑91 (MQ=255)

           aaCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAg                                                                                 >  1:1783145/1‑91 (MQ=255)

            aCGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGa                                                                                <  1:1501326/91‑1 (MQ=255)

               gTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATgg                                                                             >  1:2620365/1‑91 (MQ=255)

                tATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGc                                                                            <  2:2241164/91‑1 (MQ=255)

                  tCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCtt                                                                          >  1:3865986/1‑91 (MQ=255)

                    ttAGCTGCTGAAGGTGCAGAAGTTTTAGAAGGAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCt                                                                        <  2:1391619/91‑1 (MQ=255)

                    ttAGCTACTGAAGGTGCAGACGTTTTAGACGCAACAGCCTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGGTGGCTTCt                                                                        >  2:3440802/1‑91 (MQ=255)

                    ttAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCt                                                                        <  1:2639643/91‑1 (MQ=255)

                    ttAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCt                                                                        >  1:3263470/1‑91 (MQ=255)

                     tAGCTACTGAAGTTGCAGAAGTTTTAGAAGCAAAAGACGGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTa                                                                       <  2:2356393/91‑1 (MQ=255)

                     tAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTa                                                                       <  2:2030128/91‑1 (MQ=255)

                     tAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTa                                                                       <  1:755034/91‑1 (MQ=255)

                     tAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTa                                                                       <  2:1111817/91‑1 (MQ=255)

                      aGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTac                                                                      >  2:2579725/1‑91 (MQ=255)

                      aGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTac                                                                      <  1:2113426/91‑1 (MQ=255)

                      aGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTac                                                                      >  1:2333446/1‑91 (MQ=255)

                       gCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTaca                                                                     >  1:1599863/1‑91 (MQ=255)

                       gCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTaca                                                                     <  2:2643363/91‑1 (MQ=255)

                       gCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTaca                                                                     <  2:2194309/91‑1 (MQ=255)

                            tGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATc                                                                >  1:2174132/1‑91 (MQ=255)

                            tGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATc                                                                >  1:225460/1‑91 (MQ=255)

                             gAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCg                                                               >  2:3600596/1‑91 (MQ=255)

                             gAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCg                                                               >  1:905503/1‑91 (MQ=255)

                                ggTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTAc                                                            >  2:2780648/1‑91 (MQ=255)

                                ggTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTAc                                                            >  1:2413556/1‑91 (MQ=255)

                                ggTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTAc                                                            >  1:3527571/1‑91 (MQ=255)

                                 gTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACg                                                           <  2:212295/91‑1 (MQ=255)

                                  tGCAGAAGTCTTAGAAGCAACAGACTGGGGTAAACGTCGCCTCGCTTATGAAATCAATGCTTTCAACGATGGCTTCTACAACATCGTACgt                                                          >  2:1769386/1‑91 (MQ=255)

                                   gCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtg                                                         >  1:1853571/1‑91 (MQ=255)

                                    cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                        <  1:1421601/91‑1 (MQ=255)

                                    cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                        >  1:1129369/1‑91 (MQ=255)

                                    cAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACgtgt                                                        <  2:3942425/91‑1 (MQ=255)

                                       aaGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTaa                                                     >  1:2622325/1‑91 (MQ=255)

                                          ttttAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATc                                                  >  1:1689471/1‑91 (MQ=255)

                                           tttAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCt                                                 <  2:3483077/91‑1 (MQ=255)

                                            ttAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTg                                                <  2:2423001/91‑1 (MQ=255)

                                             tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTgg                                               >  2:3794034/1‑90 (MQ=255)

                                             tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                               >  1:890778/1‑91 (MQ=255)

                                             tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                               <  1:180611/91‑1 (MQ=255)

                                             tAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGa                                               >  1:2289986/1‑91 (MQ=255)

                                              aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAACGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  2:3900515/1‑91 (MQ=255)

                                              aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  2:1181237/1‑91 (MQ=255)

                                              aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  2:1445546/1‑91 (MQ=255)

                                              aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  2:1538424/1‑91 (MQ=255)

                                              aGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGAt                                              >  1:723575/1‑91 (MQ=255)

                                               gAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATa                                             <  1:3692570/91‑1 (MQ=255)

                                                 aGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATaac                                           <  1:3208418/91‑1 (MQ=255)

                                                     aaaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAacaaca                                       >  1:491073/1‑91 (MQ=255)

                                                      aaGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAAcaacaa                                      <  1:1837671/91‑1 (MQ=255)

                                                       aGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACaaa                                     <  1:757475/91‑1 (MQ=255)

                                                           tGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTa                                 >  1:882736/1‑91 (MQ=255)

                                                             gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                               >  1:2523739/1‑91 (MQ=255)

                                                             gggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACt                               <  2:1821390/91‑1 (MQ=255)

                                                              ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                              >  1:1650782/1‑91 (MQ=255)

                                                              ggTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTg                              >  2:1718438/1‑91 (MQ=255)

                                                                 aaaCGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACg                           <  2:3712016/91‑1 (MQ=255)

                                                                 aaaCGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACg                           <  2:3077133/91‑1 (MQ=255)

                                                                  aaCGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGa                          >  1:3804167/1‑91 (MQ=255)

                                                                   aCGTCGCCTAGCTTATGAAATCAATGATTTAAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGaa                         <  2:3144294/91‑1 (MQ=255)

                                                                      tCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                      >  1:1873725/1‑91 (MQ=255)

                                                                      tCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTc                      >  1:2937824/1‑91 (MQ=255)

                                                                       cGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTcc                     >  1:2480628/1‑91 (MQ=255)

                                                                        gCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCa                    >  2:3417607/1‑91 (MQ=255)

                                                                        gCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCa                    >  2:654557/1‑91 (MQ=255)

                                                                        gCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCa                    >  2:1884673/1‑91 (MQ=255)

                                                                         ccTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCaa                   >  2:782422/1‑91 (MQ=255)

                                                                            aGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGt                >  1:1014145/1‑91 (MQ=255)

                                                                              cTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCt              <  1:1294385/91‑1 (MQ=255)

                                                                               ttATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTa             <  1:947115/91‑1 (MQ=255)

                                                                               ttATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTa             >  2:3634220/1‑91 (MQ=255)

                                                                                tATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAg            <  1:508061/91‑1 (MQ=255)

                                                                                 aTGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGc           <  1:2858598/91‑1 (MQ=255)

                                                                                  tGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCt          <  1:2010668/91‑1 (MQ=255)

                                                                                   gAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTa         >  1:1077651/1‑91 (MQ=255)

                                                                                   gAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTa         <  1:3275695/91‑1 (MQ=255)

                                                                                    aaaTCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaa        <  2:743383/91‑1 (MQ=255)

                                                                                      aTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaaa      >  2:3570821/1‑91 (MQ=255)

                                                                                      aTCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTaaaa      >  1:3036604/1‑91 (MQ=255)

                                                                                       tCATTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAACGCTACTGACGAATTCCAACGTCTAGCTacact     >  2:582067/1‑87 (MQ=255)

                                                                                         aaTGATTTCAAAGATGGCGTCTACAACATCGTACGTGTTAAGTCTGATAACAGCAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCa   <  2:1311596/91‑1 (MQ=255)

                                                                                          aTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAg  <  2:3755120/91‑1 (MQ=255)

                                                                                          aTGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAg  >  2:1590416/1‑91 (MQ=255)

TTGAACGTTTCAACGGTATCTTAGCTACTGAAGGTGCAGAAGTTTTAGAAGCAAAAGACTGGGGTAAACGTCGCCTAGCTTATGAAATCAATGATTTCAAAGATGGCTTCTACAACATCGTACGTGTTAAATCTGATAACAACAAAGCTACTGACGAATTCCAACGTCTAGCTAAAATCAG  >  USA300TCH1516_ALE/418735‑418915

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`