BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,726,545	(C)_6→5	noncoding (2640/2904 nt)	rrlG ←	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,726,545	0	C	.	76.8%	627.5 / 282.7	233	noncoding (2640/2904 nt)	rrlG	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base C (24/30); new base . (81/98); total (105/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 233 total aligned reads displayed.

TCTTCCGCATTTAGCTACCGGGCAGTGCCATTGGCATGACAACCCGAACACCAGTGATGCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCC  >  NC_000913/2726449‑2726642

tccgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGt                                                                                              <  2:4346444‑M1/10‑1 (MQ=255)

tccgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGt                                                                                              <  1:6521970‑M1/10‑1 (MQ=255)

  cgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTtc                                                                                            <  1:4346448‑M1/12‑1 (MQ=255)

  cgcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTtc                                                                                            <  1:5393679‑M1/12‑1 (MQ=255)

   gcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTtct                                                                                           <  1:6581852‑M1/13‑1 (MQ=255)

   gcatttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTtct                                                                                           <  1:5393682‑M1/13‑1 (MQ=255)

    catttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTtctc                                                                                          >  2:5393380‑M1/88‑101 (MQ=255)

     atttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTcc                                                                                         >  2:4346215‑M1/87‑101 (MQ=255)

      tttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCa                                                                                        >  2:5393432‑M1/86‑101 (MQ=255)

      tttagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCa                                                                                        <  2:6329045‑M1/16‑1 (MQ=255)

        tagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAgc                                                                                      <  1:6521823‑M1/18‑1 (MQ=255)

        tagctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAgc                                                                                      <  2:5393689‑M1/18‑1 (MQ=255)

         agctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAgcg                                                                                     <  2:6581929‑M1/19‑1 (MQ=255)

           ctaccgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGcc                                                                                   <  2:335315‑M1/21‑1 (MQ=255)

             accgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCa                                                                                 <  2:5393695‑M1/23‑1 (MQ=255)

             accgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCa                                                                                 >  2:6522265‑M1/79‑101 (MQ=255)

             accgggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCa                                                                                 <  1:335094‑M1/23‑1 (MQ=255)

                gggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACgg                                                                              <  1:6328816‑M1/26‑1 (MQ=255)

                gggcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACgg                                                                              >  1:6522260‑M1/76‑101 (MQ=255)

                  gcagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCa                                                                            <  2:5393702‑M1/28‑1 (MQ=255)

                   cagtgccattggcatgacaaccctaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAg                                                                           >  1:4346316‑M1/73‑101 (MQ=255)

                   cagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAg                                                                           >  1:5393573‑M1/73‑101 (MQ=255)

                   cagtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAg                                                                           <  2:6193115‑M1/29‑1 (MQ=255)

                    agtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAgcg                                                                                     <  1:4346472‑M1/19‑1 (MQ=255)

                    agtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGa                                                                          <  2:6521797‑M1/30‑1 (MQ=255)

                    agtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGa                                                                          <  2:6581770‑M1/30‑1 (MQ=255)

                     gtgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGAt                                                                         <  2:6192947‑M1/31‑1 (MQ=255)

                      tgccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATa                                                                        <  2:6521794‑M1/32‑1 (MQ=255)

                       gccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAg                                                                       <  1:6582020‑M1/33‑1 (MQ=255)

                       gccattggcatgacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAg                                                                       >  1:4346105‑M1/69‑101 (MQ=255)

                                  gacaacccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTg                                                            >  2:4346100‑M1/58‑101 (MQ=255)

                                      acccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtct                                                         <  1:6521998‑M1/48‑1 (MQ=255)

                                      acccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtctc                                                        >  1:5393428‑M1/54‑101 (MQ=255)

                                      acccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtctc                                                        >  2:5393466‑M1/54‑101 (MQ=255)

                                      acccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtctc                                                        <  2:6193071‑M1/48‑1 (MQ=255)

                                       cccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCa                                                       >  2:5393530‑M1/53‑101 (MQ=255)

                                       cccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCa                                                       >  2:4346252‑M1/53‑101 (MQ=255)

                                        ccgaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCac                                                      <  2:6582161‑M1/50‑1 (MQ=255)

                                          gaacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCacga                                                    >  2:4345904‑M1/50‑101 (MQ=255)

                                           aacaccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCacgac                                                   <  1:5393738‑M1/53‑1 (MQ=255)

                                             caccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGt                                                 <  2:4346502‑M1/55‑1 (MQ=255)

                                              accagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGtt                                                >  1:4345985‑M1/46‑101 (MQ=255)

                                               ccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtc                                                          <  2:335476‑M1/46‑1 (MQ=255)

                                               ccagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTc                                               >  2:4346163‑M1/45‑101 (MQ=255)

                                                cagtgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCt                                              <  1:5393744‑M1/58‑1 (MQ=255)

                                                   tgatgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTaaa                                           <  1:335384‑M1/61‑1 (MQ=255)

                                                     atgcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAAcc                                         <  2:4346508‑M1/63‑1 (MQ=255)

                                                       gcgtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCa                                       >  1:4346211‑M1/37‑101 (MQ=255)

                                                         gtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGc                                     <  1:335180‑M1/67‑1 (MQ=255)

                                                         gtccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGc                                     >  1:6522212‑M1/35‑101 (MQ=255)

                                                          tccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCt                                    <  1:4346516‑M1/68‑1 (MQ=255)

                                                           ccactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTc                                   >  2:335506‑M1/33‑101 (MQ=255)

                                                            cactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCGCGGCAGATAGGGACCGAACTGTCTCACGACGTTCt                                              >  2:4346134‑M1/32‑89 (MQ=255)

                                                            cactccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTcg                                  >  2:5393495‑M1/32‑101 (MQ=255)

                                                              ctccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGAcc                                                                  >  2:5393487‑M1/30‑67 (MQ=255)

                                                              ccccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCt                                              >  2:4345988‑M1/30‑87 (MQ=255)

                                                               tccggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGt                               >  1:5393442‑M1/29‑101 (MQ=255)

                                                                  ggtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTAcc                            <  2:6193060‑M1/76‑1 (MQ=255)

                                                                   gtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTcg                                  <  2:6192906‑M1/70‑1 (MQ=255)

                                                                   gtcctctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTAcc                            <  1:334836‑M1/76‑1 (MQ=255)

                                                                    tcctctcgtactaggagcagcccCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCa                           <  2:6521964‑M1/78‑1 (MQ=255)

                                                                       tctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACttt                       <  1:6193087‑M1/81‑1 (MQ=255)

                                                                        ctcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTa                      >  1:5393398‑M1/20‑101 (MQ=255)

                                                                         tcgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTaa                     >  2:6193342‑M1/19‑101 (MQ=255)

                                                                          cgtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTaaa                    >  2:4346125‑M1/18‑101 (MQ=255)

                                                                           gtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAAt                   >  2:6522194‑M1/17‑101 (MQ=255)

                                                                           gtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAAt                   <  1:4346536‑M1/85‑1 (MQ=255)

                                                                           gtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAAt                   >  1:6522195‑M1/17‑101 (MQ=255)

                                                                           gtactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAAt                   <  1:335020‑M1/85‑1 (MQ=255)

                                                                            tactaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATg                  >  1:5393438‑M1/16‑101 (MQ=255)

                                                                             actaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATgg                 <  1:6328763‑M1/87‑1 (MQ=255)

                                                                             actaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATgg                 <  2:4346541‑M1/87‑1 (MQ=255)

                                                                              ctaggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGc                >  1:6522187‑M1/14‑101 (MQ=255)

                                                                               taggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCg               <  1:6521871‑M1/89‑1 (MQ=255)

                                                                                aggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGa              <  1:6328743‑M1/90‑1 (MQ=255)

                                                                                 ggagcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAAcc                                         >  2:5393440‑M1/11‑73 (MQ=255)

                                                                                   agcagcccCCCCC‑TCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACa           >  2:6522180‑M1/9‑101 (MQ=255)

                                                                                     gcagccCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAg          >  1:335476‑M1/7‑101 (MQ=255)

                                                                                        gccCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCa       >  2:4346290‑M1/4‑101 (MQ=255)

                                                                                          cCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGtc                                                          >  1:6193397‑M1/2‑48 (MQ=255)

                                                                                           CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAc    <  2:6328736‑M1/101‑1 (MQ=255)

                                                                                           CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAc    >  1:6522168‑M1/1‑101 (MQ=255)

                                                                                           CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAc    <  1:6328776‑M1/101‑1 (MQ=255)

                                                                                           CCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAc    <  2:6521899‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTa                                             >  1:4346188‑M1/1‑59 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAAt                   <  2:6521585‑M1/85‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATgg                 >  1:6329189‑M1/1‑87 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  1:6328607‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  2:6328841‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  2:6328759‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   >  2:6522166‑M1/1‑101 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  2:4346557‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  2:6582037‑M1/101‑1 (MQ=255)

                                                                                            CCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAcc   <  1:6581816‑M1/101‑1 (MQ=255)

                                                                                             CCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAccc  >  2:6329172‑M1/1‑101 (MQ=255)

                                                                                             CCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAccc  <  2:6193020‑M1/101‑1 (MQ=255)

                                                                                             CCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAccc  >  1:6329173‑M1/1‑101 (MQ=255)

                                                                                             CCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAccc  <  1:6192920‑M1/101‑1 (MQ=255)

                                                                                             CCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGCTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATAccc  <  2:6581629‑M1/101‑1 (MQ=255)

                                                                                             CCCCTCAGTGCTCCAGCGCCCACGGCAGATAGGGACCGAACTTTCTCACGACGTTCTAAACCCa                                       >  2:6522217‑M1/1‑64 (MQ=255)

TCTTCCGCATTTAGCTACCGGGCAGTGCCATTGGCATGACAACCCGAACACCAGTGATGCGTCCACTCCGGTCCTCTCGTACTAGGAGCAGCCCCCCTCAGTTCTCCAGCGCCCACGGCAGATAGGGACCGAACTGTCTCACGACGTTCTAAACCCAGCTCGCGTACCACTTTAAATGGCGAACAGCCATACCC  >  NC_000913/2726449‑2726642

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