BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,260	A→T	7.6%	L16Q (CTG→CAG)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,260	0	A	T	7.6%	357.8 / 12.7	144	L16Q (CTG→CAG)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (58/69); new base T (4/8); total (65/79)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 148 total aligned reads displayed.

ACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCA  >  NC_002947/1383113‑1383409

acacGTCGCTGCGGAAGCTCAATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGAGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGt                                                                                                                                                     >  2:4633536/1‑150 (MQ=255)

   cGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcg                                                                                                                                                  <  1:209102/150‑1 (MQ=255)

    gTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcgg                                                                                                                                                 <  2:954543/150‑1 (MQ=255)

     tCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTgcggc                                                                                                                                                <  1:3481652/150‑1 (MQ=255)

        cTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                             <  2:516309/150‑1 (MQ=255)

        cTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGcc                                                                                                                                             <  1:5439648/150‑1 (MQ=255)

         tGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGccc                                                                                                                                            >  1:4484657/1‑150 (MQ=255)

         tGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGccc                                                                                                                                            >  2:4276409/1‑150 (MQ=255)

          gCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTGTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGccct                                                                                                                                           <  1:833962/150‑2 (MQ=255)

             gAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCTTGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGCACAGTGCGGCGCCCCCaa                                                                                                                                        <  1:4262852/150‑1 (MQ=255)

              aaGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCAGTGCGGCGCCCCCacg                                                                                                                                       >  1:2541568/1‑148 (MQ=255)

               aGCTCTATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAgg                                                                                                                                      <  2:2584386/150‑1 (MQ=255)

               aGCTCCATCACCTTGGGCTCGCCGACATCGATCTGCTGGGTTTGCTCGAGCTGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGTACATCTGGTTGGCGCTGCCGCTGACGAGCAGGGCGCCGCACAGTGGGGCGCCCCCAAgg                                                                                                                                      <  1:5240530/150‑1 (MQ=255)

                gCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAggg                                                                                                                                     <  2:5335182/150‑1 (MQ=255)

                    cATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTaaa                                                                                                                                 >  1:4996640/1‑150 (MQ=255)

                    cATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTaaa                                                                                                                                 >  2:3786611/1‑150 (MQ=255)

                       cACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCGAGGGTAAAggg                                                                                                                              >  1:2833998/1‑149 (MQ=255)

                         ccTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTAc                                                                                                                            >  2:4114977/1‑150 (MQ=255)

                            tCGGCGCGCCGACATCGATGTGCTGGGTTTGCTCGAGCATTTTGCGCTCGACGCGCTCCTCTTGGTCGGTGCGGTGGGACATGGGGTTGGGGCCGCCGCTGACGTGCTGGGCGCCGAACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTc                                                                                                                         <  1:2955239/150‑1 (MQ=255)

                             cGGCTCGCCGACATCGATCTGCAGGGTATGCTCGCGCAGTGTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCg                                                                                                                        <  2:2800908/150‑1 (MQ=255)

                                cTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCtt                                                                                                                     <  2:2713079/150‑1 (MQ=255)

                                cTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCtt                                                                                                                     <  2:2713072/150‑1 (MQ=255)

                                 tCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTg                                                                                                                    >  2:3593399/1‑150 (MQ=255)

                                  cGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGt                                                                                                                   <  1:3168803/150‑1 (MQ=255)

                                  cGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGt                                                                                                                   <  2:2874549/150‑1 (MQ=255)

                                   gCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTa                                                                                                                  >  2:4957591/1‑150 (MQ=255)

                                    ccGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTAc                                                                                                                 <  2:2628921/150‑1 (MQ=255)

                                     cGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAATGTACTTAGCGTGTACa                                                                                                                >  1:2525007/1‑150 (MQ=255)

                                      gACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCCGGGTAAAGGTACTTCGCTTGTACAt                                                                                                               >  1:5141881/1‑150 (MQ=255)

                                             aTCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCAcg                                                                                                        >  1:873677/1‑149 (MQ=255)

                                                tGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGt                                                                                                     <  2:1089926/150‑1 (MQ=255)

                                                    gggTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCAGGGTAAAGGTAATTCGCTTGTACATGATAACTCGTGCAt                                                                                                 >  1:185925/1‑150 (MQ=255)

                                                       tATGCTCGAGCATTTTGCGCTCGACGCGCCCCTCTTGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTgg                                                                                              <  1:2633986/150‑1 (MQ=255)

                                                       tATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTgg                                                                                              <  2:2632284/150‑1 (MQ=255)

                                                       tATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTgg                                                                                              <  2:1928379/150‑1 (MQ=255)

                                                        aTGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGt                                                                                             <  1:5303152/150‑1 (MQ=255)

                                                          gCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTa                                                                                           >  2:2250844/1‑150 (MQ=255)

                                                          gCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTa                                                                                           >  2:2250658/1‑150 (MQ=255)

                                                              gAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAg                                                                                       >  1:2018562/1‑150 (MQ=255)

                                                                   gTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTgcg                                                                                                                                                  >  2:938622/1‑86 (MQ=255)

                                                                     tggcgcACGACGCGCCCCGCATGGTGGCTGCGCTGGCACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGGGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGaa                                                                                <  1:2798505/148‑1 (MQ=255)

                                                                         gcTCGACGCGCTCCTCATGGTCGCTGCGCTGGGCCATCTGGTTGGCGCAGCCGCTGACTAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAAcgcg                                                                            <  1:863739/150‑1 (MQ=255)

                                                                         gcTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCCCCCACTGCGGCGCCCCCCACGGGAAACGTTCATGGCGTGTTCAAGGTCACACGTGGAATGGTTATGATCGGCGAAcgcg                                                                            >  1:246669/1‑150 (MQ=255)

                                                                          cTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCAGTGCGGCGCCCCCATGGGGAGAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                           >  2:3373985/1‑150 (MQ=255)

                                                                          cTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCGTGTACATGATCACTCGTGCAGTGGTTATCAGCGGCGAACGCgg                                                                           >  2:3373937/1‑150 (MQ=255)

                                                                          cTCGACGCGCTCCCCATGGTCGCTGCGCTGGGACTTCTGTTTGGCGCAGCCGCTCACGAGCAGGGCGCCGCACCGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCgg                                                                           <  1:1112956/150‑1 (MQ=255)

                                                                            cGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGc                                                                         <  2:3725168/150‑1 (MQ=255)

                                                                            cGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGc                                                                         <  2:3909262/150‑1 (MQ=255)

                                                                             gACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACACTGCGGCGCCCCCAAGGGTAAAGGTAATTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGcc                                                                        >  1:1159426/1‑150 (MQ=255)

                                                                               cgcgCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTCAAGGTACTTCGCTTGTACATGATCACTCGTGCCTTGGTTATCAGCGGCGAACGCGGGCCTg                                                                      >  1:4349685/1‑150 (MQ=255)

                                                                               cgcgCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTg                                                                      <  1:3593399/150‑1 (MQ=255)

                                                                                gcgcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCCGCCGCTGACGAGCTGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGc                                                                     <  1:4314433/150‑1 (MQ=255)

                                                                                  gcTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAg                                                                   >  1:3050063/1‑150 (MQ=255)

                                                                                      ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAg                                                               >  1:4377437/1‑150 (MQ=255)

                                                                                      ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAg                                                               <  1:1716646/150‑1 (MQ=255)

                                                                                      ctcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACGTCGCTTGTACATGATCACTCGTGCATTGGTTATAAGCGGCGAACGCGGGCCTGCAGGAAg                                                               >  1:1095480/1‑150 (MQ=255)

                                                                                       tcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCCGTGCGGCGCCCCCCAGGGGTAAGGTACTTCGCTTGTTCATGATCACTCGTGCAGTGGTTTTCAGCGGCGCACGCGGGCCGGCAGGAAgg                                                              >  1:3349429/1‑150 (MQ=255)

                                                                                       tcATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCCGTGCGGCGCCCCCAAGGGTAAAGGGACGTCGCTTGTACATGATCACAAGGGCATTGGTTATCAGCGGCGAACGCCGGCCGGCAGGAAgg                                                              >  1:3349669/1‑150 (MQ=255)

                                                                                          tGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCa                                                           <  1:3521560/150‑1 (MQ=255)

                                                                                          tGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCa                                                           <  2:3044460/150‑1 (MQ=255)

                                                                                          tGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCa                                                           <  1:369858/150‑1 (MQ=255)

                                                                                           ggTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCCCACTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACACCTGCAATGGTTATAAGCCGCGAACCCCGGCCTGCAGGAAGGtacg                                                          >  1:502845/1‑147 (MQ=255)

                                                                                               gctgcgctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCCAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACt                                                      >  2:2415229/1‑150 (MQ=255)

                                                                                               gctgcgctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACt                                                      >  2:1203451/1‑150 (MQ=255)

                                                                                                ctgcgctgGTACATCTGGTTGGCGCAGCCGCTGACGAGCTGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACtt                                                     <  1:1386885/150‑1 (MQ=255)

                                                                                                ctgcgctgGGACATCTGGTTGGCGCAGCCGCTGACTAGCAGGGCGCCGCACAGAGCGGCGCCCCCAAGGTTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACtt                                                     <  1:4198005/150‑1 (MQ=255)

                                                                                                  gcggtgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCa                                                   <  1:3915352/150‑1 (MQ=255)

                                                                                                  gcgctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCa                                                   <  2:2008404/150‑1 (MQ=255)

                                                                                                   cgctgGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCAATCGTGCATTGGTTATAAGCGGCGAACGCGGGCCTGCAGGAAGGTAAGCAcataat                                                  >  1:1492003/1‑145 (MQ=255)

                                                                                                        ggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGccct                                                                                                                                           <  1:5357566/56‑2 (MQ=255)

                                                                                                        ggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACTGTGCGGCGccct                                                                                                                                           >  2:5357566/1‑55 (MQ=255)

                                                                                                        ggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGc                                             <  2:5041373/150‑1 (MQ=255)

                                                                                                        ggACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGc                                             <  2:1236455/150‑1 (MQ=255)

                                                                                                        ggACATCTGGTGGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGc                                             <  2:5041377/150‑1 (MQ=255)

                                                                                                         gACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCg                                            <  2:3800490/150‑1 (MQ=255)

                                                                                                               tGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGGGGCGCCCCCAAGGGTAAAGGTACTTGGCTTGTACATGATGACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACg                                      >  2:4662240/1‑150 (MQ=255)

                                                                                                                 gTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACCGTGCGGCGCCCCCCAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAATGAGCGGCAACGcc                                    >  2:5148022/1‑150 (MQ=255)

                                                                                                                 gTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCAGTGAGCGGCAACGcc                                    >  2:4072516/1‑150 (MQ=255)

                                                                                                                    ggCGCCGCCGCGGCCGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 <  1:711279/150‑1 (MQ=255)

                                                                                                                      cgcAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGc                               <  2:1559014/150‑1 (MQ=255)

                                                                                                                         aGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTc                            <  1:3675996/150‑1 (MQ=255)

                                                                                                                                  cGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggt                   <  2:748304/150‑1 (MQ=255)

                                                                                                                                     gCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcg                >  1:2786576/1‑150 (MQ=255)

                                                                                                                                      cAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcgg               >  2:3011033/1‑150 (MQ=255)

                                                                                                                                        gggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 <  1:4956086/130‑1 (MQ=255)

                                                                                                                                        gggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTgg                                 >  2:4956086/1‑130 (MQ=255)

                                                                                                                                        gggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTcggtgcggtg             >  2:5400435/1‑150 (MQ=255)

                                                                                                                                         ggCGCCGCACAGTGCGGCGCCCCCCAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATGGGTTATCAGCGGCGAACGCGGGCCTGCAGGACGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            >  2:1703843/1‑150 (MQ=255)

                                                                                                                                         ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            >  2:2117824/1‑150 (MQ=255)

                                                                                                                                         ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            >  2:1632713/1‑150 (MQ=255)

                                                                                                                                         ggCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCggtgcggtgc            >  2:4725295/1‑150 (MQ=255)

                                                                                                                                          gcgcCGCACTGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCt           >  2:4797470/1‑150 (MQ=255)

                                                                                                                                           cgccgcCCAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          <  1:3101638/150‑1 (MQ=255)

                                                                                                                                           cgccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          <  2:3991308/150‑1 (MQ=255)

                                                                                                                                           cgccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCtt          <  1:4276409/150‑1 (MQ=255)

                                                                                                                                             ccgcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGt        >  2:2117976/1‑150 (MQ=255)

                                                                                                                                               gcACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTAc      <  2:1424633/150‑1 (MQ=255)

                                                                                                                                                cacaGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACt     >  1:2398007/1‑150 (MQ=255)

                                                                                                                                                 acaGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTc    <  1:3952294/150‑1 (MQ=255)

                                                                                                                                                   aGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCa  >  1:815720/1‑150 (MQ=255)

ACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACATCTGGTTGGCGCAGCCGCTGACGAGCAGGGCGCCGCACAGTGCGGCGCCCCCAAGGGTAAAGGTACTTCGCTTGTACATGATCACTCGTGCATTGGTTATCAGCGGCGAACGCGGGCCTGCAGGAAGGTCAGCACTTCATTGAGCGGCAACGCCTGGGCGTCCTGCTCGGTGCGGTGCTTGTACTCCA  >  NC_002947/1383113‑1383409

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