BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,730,407	A→C	7.7%	L62R (CTG→CGG)	PP_1519 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,730,407	0	A	C	7.7%	896.9 / 17.3	322	L62R (CTG→CGG)	PP_1519	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (125/163); new base C (32/2); total (157/165)
Fisher's exact test for biased strand distribution p-value = 4.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

GGTGTTCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAG  >  NC_002947/1730258‑1730555

ggtgTTCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCa                                                                                                                                                      <  2:6340758/150‑1 (MQ=255)

  tgtTCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAgg                                                                                                                                                    <  2:4781099/150‑1 (MQ=255)

  tgtTCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAgg                                                                                                                                                    <  1:2010413/150‑1 (MQ=255)

    ttCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGt                                                                                                                                                  >  2:6843160/1‑150 (MQ=255)

    ttCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGt                                                                                                                                                  >  2:3614540/1‑150 (MQ=255)

       ggTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCAc                                                                                                                                               >  1:4371055/1‑150 (MQ=255)

       ggTATTGATGAACAGCAACTTGCCAGTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCAc                                                                                                                                               <  1:3230450/150‑1 (MQ=255)

        gTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACg                                                                                                                                              <  2:6424042/150‑1 (MQ=255)

         tATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACgg                                                                                                                                             >  2:6601984/1‑150 (MQ=255)

         tATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACgg                                                                                                                                             <  2:392053/150‑1 (MQ=255)

            tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGTGACCCTGGCGCGCCAggccggccc                                                                                                                                             <  1:4482624/150‑8 (MQ=255)

            tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAg                                                                                                                                          <  1:1678492/150‑1 (MQ=255)

             gatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCGGGTCACGGCAgg                                                                                                                                         >  1:1600772/1‑150 (MQ=255)

             gatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAgg                                                                                                                                         <  1:3675272/150‑1 (MQ=255)

             gatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAgg                                                                                                                                         <  1:997642/150‑1 (MQ=255)

             gatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAgg                                                                                                                                         <  1:4858357/150‑1 (MQ=255)

              atgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAggg                                                                                                                                        <  1:6823176/150‑1 (MQ=255)

               tgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAgggg                                                                                                                                       >  2:6178698/1‑150 (MQ=255)

                   cAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCGGGTCACGGCAGGGGCGTc                                                                                                                                   >  2:2422907/1‑150 (MQ=255)

                   cAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTc                                                                                                                                   >  2:5686611/1‑150 (MQ=255)

                   cAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCACGGTCACGGCAGGGGCggc                                                                                                                                   >  1:639880/1‑148 (MQ=255)

                       aaCTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCGGGCGCGCCAGGGGCACGGCAGGGGCGTCCCTc                                                                                                                               >  1:1966717/1‑150 (MQ=255)

                        aCTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCa                                                                                                                              <  2:2164993/150‑1 (MQ=255)

                          ttGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGc                                                                                                                            <  1:6341000/150‑1 (MQ=255)

                           tGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCg                                                                                                                           >  1:5004675/1‑150 (MQ=255)

                              cATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACCGCAGGGGCGTCCAACAGCGaca                                                                                                                        >  1:5445033/1‑150 (MQ=255)

                              cATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCGGGCGCGCCACGGTCACGGCAGGGGCGGCCATCACCGaca                                                                                                                        >  1:1243907/1‑150 (MQ=255)

                               aTTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCAACAGCGacac                                                                                                                       >  1:1021527/1‑150 (MQ=255)

                               aTTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                       >  2:6451601/1‑150 (MQ=255)

                               aTTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                       >  2:4436644/1‑150 (MQ=255)

                               aTTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                       >  2:1771305/1‑150 (MQ=255)

                               aTTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                       >  2:2940937/1‑150 (MQ=255)

                               aTTCAGCTCAATCTTCGCACTTACCGCATAGCGGTGGCCAGGATAGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGacac                                                                                                                       >  2:2329252/1‑150 (MQ=255)

                                  cAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                                                    <  2:4958935/150‑1 (MQ=255)

                                   aGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATcc                                                                                                                   <  2:333675/150‑1 (MQ=255)

                                    tcTCAATGAGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCGGGGGCGTCCATCAGCGACACAAccg                                                                                                                  >  2:4374441/2‑149 (MQ=255)

                                     cTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCCCGGCAGGGGCGGCCATCAGCGACACAACCTg                                                                                                                 >  1:767462/1‑150 (MQ=255)

                                      tCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGc                                                                                                                <  2:1302434/150‑1 (MQ=255)

                                      tCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGc                                                                                                                >  2:6087657/1‑150 (MQ=255)

                                                  cTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTc                                                                                                    <  2:144459/150‑1 (MQ=255)

                                                     aCCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACAACCTGCAGTTCCACCCTCAgg                                                                                                 >  1:6047114/1‑150 (MQ=255)

                                                       cGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCggccgg                                                                                                                                                <  2:3823262/101‑5 (MQ=255)

                                                         cATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTgg                                                                                             <  2:3270631/150‑1 (MQ=255)

                                                           tAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                                                           <  2:1203580/150‑1 (MQ=255)

                                                               ggTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcag                                                                                       <  2:1673556/150‑1 (MQ=255)

                                                                gTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCGGGCGCGCCCGGGTCAGGGCAGGGGCGGCCCACAGCGACACCTCCTGCCGTGCCACGCGCCGGGTggcgggggg                                                                                      >  1:298753/1‑146 (MQ=255)

                                                                    ccAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑GCGGGCGCGCCCGGGGGACGGCCGGGGCGTCCAACACCGACACAAACAGCACATACACGCGCCGGGGggcggccggcggc                                                                                  >  1:572618/1‑150 (MQ=255)

                                                                     cAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCGGGTCACGGCCGGGGCGTCCATCAGCGACACAACCTGCAGTTCCACGCGCCGGGTggcggcaggcggca                                                                                 >  1:4079851/1‑150 (MQ=255)

                                                                     cAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcaggcggca                                                                                 <  2:3395724/150‑1 (MQ=255)

                                                                     cAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcaggcggca                                                                                 >  2:6541155/1‑150 (MQ=255)

                                                                     cAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcaggcggca                                                                                 <  2:5988620/150‑1 (MQ=255)

                                                                      aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGgcggcaggcggcag                                                                                <  1:4266067/150‑1 (MQ=255)

                                                                      aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGgcggcaggcggcag                                                                                <  2:1531366/150‑1 (MQ=255)

                                                                      aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGgcggcaggcggcag                                                                                <  1:5857649/150‑1 (MQ=255)

                                                                      aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCAGCAGTTCCACGCTCAGGGTGgcggcaggcggcag                                                                                >  2:188258/1‑150 (MQ=255)

                                                                        gCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACATTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAggg                                                                              <  1:6238421/150‑1 (MQ=255)

                                                                          ttGACCTGCGCGGGGTCATAGGTGAGATGGAATGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCAGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACAAACAGCAGTTCCACGCACAGGGTGGCGGCAGGCGGCAGGGCa                                                                            >  2:2256255/1‑150 (MQ=255)

                                                                          ttGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑GCTGGCGCGCCAGGGTCACCGCAGGGGCGTCCATCAGCGACACAACCAGCAGTTACACGCTCAGGGTGGGGGCAGGCGGCAGGGCa                                                                            >  1:853711/1‑150 (MQ=255)

                                                                          ttGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCCGGTCACGGCAGGGGCGGCCATAAGCGACACAAACAGCAGTTACACGCACAGGGTGGCGGCAGGCGGCCGGGCa                                                                            >  1:3420740/1‑150 (MQ=255)

                                                                               cTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCAGGCGCGCCAGGGGCACGGCCGGGGCGTGCCTCAGCGACACAACCAGCAGTTCAACCCACAGGGTGGCGGCCGGCGGCCGGGCAAAgcg                                                                       >  1:4650919/1‑150 (MQ=255)

                                                                                tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                                                      >  2:5755228/1‑150 (MQ=255)

                                                                                   gcgGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGc                                                                   <  2:5288511/150‑1 (MQ=255)

                                                                                    cgGGGTCATAGGTGAGATGGAACGACAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCa                                                                  <  2:6776853/150‑1 (MQ=255)

                                                                                      gggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAgg                                                                <  1:5952257/150‑1 (MQ=255)

                                                                                          cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTaga                                                            <  1:3614540/150‑1 (MQ=255)

                                                                                            tAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCACGGTCACGGCAGGGGCGTCCATCAGCGACACAAACAGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTagaag                                                          >  1:2128383/1‑150 (MQ=255)

                                                                                                tGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCACCGACACAACCAGCAGTTCAACGCTCAGGGTGGCGGCACGCGGCAGGGCAATGCGCTGCAGGTAGAAGACtt                                                      >  1:4126967/1‑150 (MQ=255)

                                                                                                  agaTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCa                                                    <  1:1438621/150‑1 (MQ=255)

                                                                                                   gaTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCAc                                                   <  1:4561503/150‑1 (MQ=255)

                                                                                                    aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCAcc                                                  <  1:2420664/150‑1 (MQ=255)

                                                                                                    aTGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCGGGCGCGCCCGGGGCACGGCAGGGGCGGCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACGTCAcc                                                  >  1:1090480/1‑150 (MQ=255)

                                                                                                        aaCGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTcc                                              <  1:2505688/150‑1 (MQ=255)

                                                                                                            gCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCCGGGTCACGGCAGGGGCGTCCAACAGCGACACATCCAGCAGTTACACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAATACTTCACCGTCCAGGc                                          >  2:2863422/1‑150 (MQ=255)

                                                                                                              aGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTg                                        <  2:3196128/150‑1 (MQ=255)

                                                                                                                      gTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCCAGCCTggac                                <  2:6240605/150‑3 (MQ=255)

                                                                                                                            tttGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCGGGgtgggtgg                          >  2:5577349/1‑150 (MQ=255)

                                                                                                                            tttGACCGGGCCG‑G‑‑CCTGGCGCGCCACGGTCACGGCAGGGGCGTCCATCAACGACACAAACAGCAGTTACACGCACAGGGTGGCGGCAGGCGGCAGGGCAAATAGCTGCAGGTATAATAATTAAACCTACACGCTGGCCGGGTtgggtgg                          >  1:2587607/1‑150 (MQ=255)

                                                                                                                             ttGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtggg                         <  2:3578794/150‑1 (MQ=255)

                                                                                                                             ttGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtggg                         <  1:1318013/150‑1 (MQ=255)

                                                                                                                             ttGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtgggtggg                         <  1:3596349/150‑1 (MQ=255)

                                                                                                                             ttGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACGTCACCGTCCAGGCTGGCCTGGTtgggtggg                         <  2:2761292/150‑1 (MQ=255)

                                                                                                                              tGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTgggg                        <  1:969837/150‑1 (MQ=255)

                                                                                                                                  cGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCAACAGCGACACATCCTGCAGTTACACGCTCAGGGTGGCGGCAGGCGGCAGGGCAAAGCGCTGCAGGTAGAAGACTTAACCGTCCAGGCTGGCCGGGTTGGGTGGGGTTTg                    >  2:5244952/1‑150 (MQ=255)

                                                                                                                                    ggccg‑g‑‑ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGc                  <  2:2656310/150‑1 (MQ=255)

                                                                                                                                      ccg‑g‑‑ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCTGGGCAATGCGCTGCAGGTAGCAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTg                <  1:2329252/150‑1 (MQ=255)

                                                                                                                                       cgtgacccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGc                  >  1:6801911/7‑150 (MQ=255)

                                                                                                                                       cg‑g‑‑ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTgg               >  2:4041735/1‑150 (MQ=255)

                                                                                                                                       cg‑g‑‑ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTgg               >  1:5177085/1‑150 (MQ=255)

                                                                                                                                        g‑g‑‑ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCAACAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGa              <  1:3551810/150‑1 (MQ=255)

                                                                                                                                             ccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGc            >  2:4065366/1‑150 (MQ=255)

                                                                                                                                              cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCa           <  2:772842/150‑1 (MQ=255)

                                                                                                                                              cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCa           <  1:5053244/150‑1 (MQ=255)

                                                                                                                                              cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCa           <  1:6601984/150‑1 (MQ=255)

                                                                                                                                              cTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCa           <  2:4275583/150‑1 (MQ=255)

                                                                                                                                                ggCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAgg         >  1:6510985/1‑150 (MQ=255)

                                                                                                                                                 gcgcgcAAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCTGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGTCCTGGTTGGGTGGGGTTTGGCTGGAGCAgga        >  1:3853617/1‑149 (MQ=255)

                                                                                                                                                   gcgcCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTg      >  1:3076402/1‑150 (MQ=255)

                                                                                                                                                   gcgcCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTg      >  2:196642/1‑150 (MQ=255)

                                                                                                                                                     gcCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCa    >  2:2228405/1‑150 (MQ=255)

                                                                                                                                                       cAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAg  >  2:6094384/1‑150 (MQ=255)

GGTGTTCGGTATTGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCG‑G‑‑CCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGGGTGGGGTTTGGCTGGAGCAGGCTGCAAG  >  NC_002947/1730258‑1730555

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