BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	549,745	C→T	14.3%	D139D (GAC→GAT)	tufB →	elongation factor Tu‑B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	549,745	0	C	T	14.3%	494.2 / 60.0	188	D139D (GAC→GAT)	tufB	elongation factor Tu‑B
Reads supporting (aligned to +/- strand): ref base C (73/88); new base T (8/19); total (81/107)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 189 total aligned reads displayed.

AGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGG  >  NC_002947/549597‑549891

agaaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGAcc                                                                                                                                                   <  1:2687798/150‑1 (MQ=37)

  aaCATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTg                                                                                                                                                 >  2:146997/1‑150 (MQ=37)

    cATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                                                                                               >  1:786427/1‑150 (MQ=37)

     aTGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTa                                                                                                                                              <  2:3050406/150‑1 (MQ=37)

        aTCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgac                                                                                                                                           <  1:2398082/150‑1 (MQ=37)

         tCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgacg                                                                                                                                          >  1:4907852/1‑150 (MQ=37)

              ggTGCTGCCCAGATGGACGGCGCGAACCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGACTTCCTGAACAAGTCCGACCTGGTAGACGACGCt                                                                                                                                     >  2:2539921/1‑150 (MQ=255)

              ggTGCTGCACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCt                                                                                                                                     >  1:3063394/1‑150 (MQ=2)

                   tgcACAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAact                                                                                                                                >  1:570552/1‑150 (MQ=2)

                     cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                                                                                              >  2:4097632/1‑150 (MQ=37)

                     cccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgc                                                                                                                              <  1:995926/150‑1 (MQ=37)

                     cacAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttc                                                                                                                              >  2:480945/3‑148 (MQ=2)

                      ccAGATGGACGGCGCTATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                                                                                             >  1:3158865/1‑150 (MQ=37)

                      ccAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAgctgct                                                                                                                             >  1:747200/1‑150 (MQ=37)

                       cAGATGGACGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAacttctc                                                                                                                            <  1:2741506/150‑2 (MQ=2)

                        aGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTgg                                                                                                                           <  1:2561002/150‑1 (MQ=37)

                          aTGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGaa                                                                                                                         <  1:2237055/150‑1 (MQ=38)

                          aTGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGaa                                                                                                                         <  2:1133813/150‑1 (MQ=38)

                               cGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGt                                                                                                                    >  2:3686878/1‑150 (MQ=39)

                               cGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGt                                                                                                                    >  1:5249573/1‑150 (MQ=39)

                               cGGCGCGATCCTGGTTGTTTCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGt                                                                                                                    <  1:2930023/150‑1 (MQ=255)

                                ggCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGCCGCTGAGCTGCTGGAACTGGTc                                                                                                                   <  1:2335151/150‑1 (MQ=255)

                                ggCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGACCTggcc                                                                                                                   <  2:2360943/150‑3 (MQ=37)

                                ggCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTc                                                                                                                   <  1:163001/150‑1 (MQ=39)

                                  cgcgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCga                                                                                                                 <  1:1348759/150‑1 (MQ=39)

                                    cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                               <  2:3980692/150‑1 (MQ=39)

                                    cgATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCgaga                                                                                                               <  2:5103261/150‑1 (MQ=39)

                                     gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                                                                                              <  2:1052631/150‑1 (MQ=39)

                                     gATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGAt                                                                                                              <  2:2230726/150‑1 (MQ=39)

                                       tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                                                                                            <  2:1998407/150‑1 (MQ=39)

                                       tCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATgg                                                                                                            <  1:3319624/150‑1 (MQ=39)

                                         cTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGaa                                                                                                          <  2:354782/150‑1 (MQ=39)

                                            gTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGtt                                                                                                       <  1:3229547/150‑1 (MQ=39)

                                             tttGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTTGTCGAGATGGAAGTTc                                                                                                      >  2:2737200/1‑150 (MQ=255)

                                             tttGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTc                                                                                                      <  1:1767351/150‑1 (MQ=39)

                                               tGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                    <  2:5383447/150‑1 (MQ=39)

                                               gtttCGGCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTcgc                                                                                                    <  2:1239721/147‑1 (MQ=255)

                                                gCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgcg                                                                                                   <  2:3253137/150‑1 (MQ=39)

                                                 cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                                                                                                  <  2:2085160/150‑1 (MQ=39)

                                                 cTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGa                                                                                                  >  2:1282437/1‑150 (MQ=39)

                                                   cGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGAcc                                                                                                >  1:1626241/1‑150 (MQ=39)

                                                    ggCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACct                                                                                               >  1:1666801/1‑150 (MQ=39)

                                                    ggCAGCTGACGGCCCGATGCCACAGACCCGTGAGCACATCCTGCTCGCTCGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTgacgatg                                                                                                                                       <  2:2965963/110‑3 (MQ=255)

                                                          cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                         >  1:2884663/1‑150 (MQ=39)

                                                          cgATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTc                                                                                         >  1:194483/1‑150 (MQ=39)

                                                            aTGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCAGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCa                                                                                       <  1:591917/150‑1 (MQ=255)

                                                             tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                      <  2:385246/150‑1 (MQ=39)

                                                             tGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCAc                                                                                      <  2:350656/150‑1 (MQ=39)

                                                               gTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCt                                                                                    >  1:322605/1‑150 (MQ=39)

                                                                  cGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGt                                                                                                                                               >  1:4953369/1‑88 (MQ=255)

                                                                   gATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGa                                                                                <  2:2301326/150‑1 (MQ=39)

                                                                        cGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTcc                                                                           >  1:1487984/1‑150 (MQ=38)

                                                                         gCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTccc                                                                          <  2:2006377/150‑1 (MQ=38)

                                                                               ccGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGcga                                                                    >  1:585500/1‑150 (MQ=37)

                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACctgct                                                                                            >  2:78147/1‑120 (MQ=39)

                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGacac                                                              >  1:1100095/1‑150 (MQ=37)

                                                                                     gCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGacac                                                              >  2:4007078/1‑150 (MQ=37)

                                                                                      cacaTCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACt                                                             >  1:432520/1‑150 (MQ=37)

                                                                                          tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                         <  2:4314122/150‑1 (MQ=37)

                                                                                          tCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                         >  2:3089821/1‑150 (MQ=37)

                                                                                                 gTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcat                                                  <  1:2672153/150‑1 (MQ=34)

                                                                                                 cgctcGTCAGGTTGGCGTTCCTGCGATCGTCGTGTTCCTGAACAAGTGCGATCAGGTTGACGATGCTGAACTTCTCGAACTCGTAGAACTCGAAGTTCGCGAACTGCTTTCGAAGTACGACTTCCCAGGCGACGAAGTTCCAatcatcaa                                                  <  1:4039689/146‑2 (MQ=2)

                                                                                                  tCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGatcatcatc                                                 >  2:501842/1‑150 (MQ=34)

                                                                                                    ccGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCgg                                               >  2:4554927/1‑150 (MQ=34)

                                                                                                            tAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAgac                                                                                                                                           >  2:3121370/1‑50 (MQ=38)

                                                                                                              ggCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGt                                     >  2:1510989/1‑150 (MQ=34)

                                                                                                                 gTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATg                                  <  1:5053137/150‑1 (MQ=34)

                                                                                                                   tCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGc                                >  2:2653622/1‑150 (MQ=34)

                                                                                                                      gTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                             <  1:3960075/150‑1 (MQ=37)

                                                                                                                      gTACATCGTGGTCGTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCt                             <  2:2180445/150‑1 (MQ=37)

                                                                                                                       tACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGATCCGCTCGTATGGCCCTg                            >  1:2793184/1‑150 (MQ=37)

                                                                                                                         cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                          <  1:1097917/150‑1 (MQ=37)

                                                                                                                         cATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGa                          <  1:5299957/150‑1 (MQ=37)

                                                                                                                          aTCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGaa                         <  1:722588/150‑1 (MQ=37)

                                                                                                                           tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                        <  2:2142348/150‑1 (MQ=37)

                                                                                                                           tCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                        <  1:210207/150‑1 (MQ=37)

                                                                                                                               ggTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCaa                    <  2:1188918/150‑1 (MQ=37)

                                                                                                                                 tCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCCTCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAg                  >  1:4487655/1‑150 (MQ=37)

                                                                                                                                 tCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAg                  >  1:885793/1‑150 (MQ=37)

                                                                                                                                  cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                 <  1:2562115/150‑1 (MQ=37)

                                                                                                                                  cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                 >  2:1825209/1‑150 (MQ=37)

                                                                                                                                  cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                 <  1:3607060/150‑1 (MQ=37)

                                                                                                                                  cTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAga                 <  2:3158865/150‑1 (MQ=37)

                                                                                                                                   ttCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGa                                                         >  1:1241140/1‑109 (MQ=35)

                                                                                                                                      cTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAgacgac             <  2:1337406/150‑1 (MQ=37)

                                                                                                                                        gatcAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                    >  2:4277122/4‑61 (MQ=255)

                                                                                                                                        gAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTcgc                                                                                                    <  1:4277122/61‑1 (MQ=255)

                                                                                                                                         aaCAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAAc          >  1:4101659/1‑150 (MQ=37)

                                                                                                                                           cAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGa        >  1:4601671/1‑150 (MQ=37)

                                                                                                                                            aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                  <  2:2509401/75‑1 (MQ=39)

                                                                                                                                            aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTAc                                                                                  >  1:2509401/1‑75 (MQ=39)

                                                                                                                                            aaGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGaa       >  2:3884077/1‑150 (MQ=37)

                                                                                                                                              ggCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAAt     <  1:4802643/150‑1 (MQ=37)

                                                                                                                                                ccGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATgg   <  1:1206205/150‑1 (MQ=35)

                                                                                                                                                 cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg  >  1:1632586/1‑150 (MQ=35)

                                                                                                                                                 cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg  >  2:1451477/1‑150 (MQ=35)

                                                                                                                                                 cGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATggg  <  2:3577833/150‑1 (MQ=35)

                                                                                                                                                   aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                        <  1:3421201/126‑1 (MQ=35)

                                                                                                                                                   aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                        >  2:3421201/1‑126 (MQ=35)

                                                                                                                                                   aCCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAg                        >  2:3421160/1‑126 (MQ=35)

AGAACATGATCACCGGTGCTGCCCAGATGGACGGCGCGATCCTGGTTTGCTCGGCCGCCGATGGTCCGATGCCGCAAACCCGTGAGCACATCCTGCTGTCCCGTCAGGTAGGCGTTCCGTACATCGTGGTCTTCCTGAACAAGGCCGACCTGGTAGACGACGCTGAGCTGCTGGAACTGGTCGAGATGGAAGTTCGCGACCTGCTGTCCACCTACGACTTCCCAGGCGACGACACTCCGATCATCATCGGTTCCGCTCGTATGGCCCTGGAAGGCAAAGACGACAACGAAATGGG  >  NC_002947/549597‑549891

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