BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	USA300TCH1516_ALE	1,089,908	T→G	5.5%	intergenic (+61/‑364)	USA300HOU_RS05380 → / → rlmI	hypothetical protein/Ribosomal RNA large subunit methyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	1,089,908	0	T	G	5.5%	743.8 / 28.7	221	intergenic (+61/‑364)	USA300HOU_RS05380/rlmI	hypothetical protein/Ribosomal RNA large subunit methyltransferase I
Reads supporting (aligned to +/- strand): ref base T (111/95); new base G (8/4); total (122/99)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 229 total aligned reads displayed.

ATTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAA  >  USA300TCH1516_ALE/1089819‑1089997

aTTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGata                                                                                          >  1:3818086/1‑91 (MQ=255)

aTTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGata                                                                                          >  1:1782712/1‑91 (MQ=255)

  tCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGatata                                                                                        <  1:3558373/91‑1 (MQ=255)

     ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGagatatat                                                                                     <  2:2687781/91‑3 (MQ=255)

     ttACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtct                                                                                     >  1:2417130/1‑91 (MQ=255)

      tACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATAtctc                                                                                    <  2:3751554/91‑1 (MQ=255)

         aaTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGGTGTGTACGCTATAAGTTAGATATATCTCTaa                                                                                 >  2:270898/1‑91 (MQ=255)

         aaTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTaa                                                                                 <  2:3230027/91‑1 (MQ=255)

         aaTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAGATATATCTaa                                                                                 <  1:42259/91‑1 (MQ=255)

          aTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAac                                                                                <  2:2912045/91‑1 (MQ=255)

          aTTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAac                                                                                >  2:823939/1‑91 (MQ=255)

            tACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAactt                                                                              <  2:1975715/91‑1 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                            <  1:1912764/91‑1 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                            >  1:3339799/1‑91 (MQ=255)

              cGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttac                                                                            >  2:1386406/1‑91 (MQ=255)

                agagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttactt                                                                          >  2:2106404/1‑91 (MQ=255)

                agagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAacttactt                                                                          >  2:674475/1‑91 (MQ=255)

                agagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAGATATATCTAacttatag                                                                          <  2:1608656/91‑4 (MQ=39)

                agagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGAGATATATCTAacttatag                                                                          <  1:759792/91‑4 (MQ=38)

                 gagCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAActtactta                                                                         <  2:1368038/91‑1 (MQ=255)

                  agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTACCTTACATAg                                                                        >  1:3488347/1‑91 (MQ=255)

                  agCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAg                                                                        >  2:3513225/1‑91 (MQ=255)

                   gCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGa                                                                       <  1:2193027/91‑1 (MQ=255)

                     gACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATa                                                                     <  2:910148/91‑1 (MQ=255)

                       ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                   >  1:1785310/1‑91 (MQ=255)

                       ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                   >  1:99876/1‑91 (MQ=255)

                       ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                   >  1:25157/1‑91 (MQ=255)

                       ccATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACa                                                                   >  2:2627264/1‑91 (MQ=255)

                         aTAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAgg                                                                 <  2:1582302/91‑1 (MQ=255)

                          tAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGt                                                                >  2:930231/1‑91 (MQ=255)

                             ttAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATAGAGGTCaa                                                             >  2:3454068/1‑91 (MQ=255)

                              tAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATTTCTAACTTACTTAGATACAGGTCAAt                                                            <  2:2212639/91‑1 (MQ=255)

                                aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGa                                                          >  1:670438/1‑91 (MQ=255)

                                aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGa                                                          >  2:1143847/1‑91 (MQ=255)

                                aaaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGa                                                          >  1:3808653/1‑91 (MQ=255)

                                 aaGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGaa                                                         <  1:574073/91‑1 (MQ=255)

                                   gCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGt                                                       <  1:2600278/91‑1 (MQ=255)

                                    cATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGtt                                                      >  1:1623628/1‑91 (MQ=255)

                                      ttATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTa                                                    >  1:1994006/1‑91 (MQ=255)

                                       tataAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTAt                                                   >  2:2754259/1‑91 (MQ=255)

                                          aaaaGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGa                                                <  2:1993492/91‑1 (MQ=255)

                                            aaGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATa                                              <  1:52040/91‑1 (MQ=255)

                                             aGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAg                                             <  2:2345476/91‑1 (MQ=255)

                                              gTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGt                                            <  1:2415546/91‑1 (MQ=255)

                                               tACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTa                                           <  2:2494861/91‑1 (MQ=255)

                                                aCTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTAc                                          <  1:316440/91‑1 (MQ=255)

                                                 ctatctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACt                                         >  1:734065/1‑91 (MQ=255)

                                                     ctatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttttt                                     >  2:876745/1‑91 (MQ=255)

                                                      tatTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACtttttt                                    <  1:328711/91‑1 (MQ=255)

                                                       atTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttttttt                                   >  1:1624810/1‑91 (MQ=255)

                                                       atTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACttttttt                                   <  2:1326065/91‑1 (MQ=255)

                                                        ttAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTg                                  >  2:3156099/1‑91 (MQ=255)

                                                          aaaCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTa                                <  2:117910/91‑1 (MQ=255)

                                                            aCATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAAc                              >  2:2362570/1‑91 (MQ=255)

                                                             cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAGCTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                             <  2:1589571/91‑1 (MQ=255)

                                                             cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                             <  1:3710276/91‑1 (MQ=255)

                                                             cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                             >  1:2032097/1‑91 (MQ=255)

                                                             cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                             <  1:2280322/91‑1 (MQ=255)

                                                             cATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACt                             >  2:2204633/1‑91 (MQ=255)

                                                               ttttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAg                           >  2:609327/1‑91 (MQ=255)

                                                               ttttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAg                           >  2:2451576/1‑91 (MQ=255)

                                                               ttttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAg                           >  1:3384942/1‑91 (MQ=255)

                                                                tttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTACGGATAGTACTTTTTTTGTAACTAGa                          >  1:3610679/1‑91 (MQ=255)

                                                                 ttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATATAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGAt                         >  1:3866254/1‑91 (MQ=255)

                                                                 ttGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAg                           <  1:2133236/89‑1 (MQ=255)

                                                                    aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                      >  2:3120621/1‑91 (MQ=255)

                                                                    aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                      >  2:418512/1‑91 (MQ=255)

                                                                    aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                      >  1:906079/1‑91 (MQ=255)

                                                                    aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                      >  1:1990121/1‑91 (MQ=255)

                                                                    aTGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttg                      >  1:229709/1‑91 (MQ=255)

                                                                     tgtgTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttga                     >  2:1774852/1‑91 (MQ=255)

                                                                     tgtgTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttga                     >  2:1278987/1‑91 (MQ=255)

                                                                      gtgtACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgat                    <  1:2723053/91‑1 (MQ=255)

                                                                         tACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgattga                 >  2:2074288/1‑91 (MQ=255)

                                                                         tACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATttgattga                 <  2:1804784/91‑1 (MQ=255)

                                                                          aCGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTtgattgat                <  1:2821790/91‑1 (MQ=255)

                                                                            gCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTc              <  1:826622/91‑1 (MQ=255)

                                                                             cTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCa             <  2:774628/91‑1 (MQ=255)

                                                                              tataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg            <  2:3457863/91‑1 (MQ=255)

                                                                              agtaAGTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg            >  1:1043447/3‑91 (MQ=39)

                                                                              agtaAGTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAg            >  1:1282028/3‑91 (MQ=39)

                                                                               ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg           <  1:61171/91‑1 (MQ=255)

                                                                               ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg           <  1:61518/91‑1 (MQ=255)

                                                                               ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg           >  1:1861337/1‑91 (MQ=255)

                                                                               ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg           <  1:61909/91‑1 (MQ=255)

                                                                               ataAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAgg           >  2:2243615/1‑91 (MQ=255)

                                                                                taAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGt          >  2:3569114/1‑91 (MQ=255)

                                                                                taAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTCTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGt          >  1:3077577/1‑91 (MQ=255)

                                                                                 aaGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTg         <  1:1386406/91‑1 (MQ=255)

                                                                                 aaGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTg         <  1:1214549/91‑1 (MQ=255)

                                                                                  aGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGa        >  2:1819041/1‑91 (MQ=255)

                                                                                   gTTAGAGATATATCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGAGAt       >  1:1153849/1‑91 (MQ=255)

                                                                                    ttAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtg      >  2:387643/1‑91 (MQ=255)

                                                                                     tAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt     >  2:2811990/1‑91 (MQ=255)

                                                                                     tAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt     <  1:967870/91‑1 (MQ=255)

                                                                                     tAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgt     <  1:216786/91‑1 (MQ=255)

                                                                                      aGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGAtgtg    >  2:825277/1‑91 (MQ=255)

                                                                                       gATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTTTAACTAGATTTGATTGATTCAGGTGATGTGa   >  2:146957/1‑91 (MQ=255)

                                                                                       gATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGa   >  1:3615868/1‑91 (MQ=255)

                                                                                        atatatCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGaa  <  1:2922078/91‑1 (MQ=255)

ATTCGTTACAATTACGAGAGCGACCATAATTAAAAGCATTATAAAAGTACTATCTATTAAACATTTTGATGTGTACGCTATAAGTTAGATATATCTCTAACTTACTTAGATACAGGTCAATGAAGTTTATGGATAGTACTTTTTTTGTAACTAGATTTGATTGATTCAGGTGATGTGAA  >  USA300TCH1516_ALE/1089819‑1089997

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: `‑`