BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,189,117	A→G	6.0%	D86G (GAC→GGC)	PP_4564 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,189,117	0	A	G	6.0%	1210.1 / 41.2	379	D86G (GAC→GGC)	PP_4564	lipoprotein
Reads supporting (aligned to +/- strand): ref base A (173/183); new base G (6/17); total (179/200)
Fisher's exact test for biased strand distribution p-value = 5.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 392 total aligned reads displayed.

CCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCC  >  NC_002947/5188969‑5189266

ccGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCcgac                                                                                                                                                      <  1:6559868/150‑1 (MQ=255)

  gTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTcggcga                                                                                                                                                    <  2:6590688/150‑1 (MQ=255)

       cGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAgg                                                                                                                                               <  2:2392678/150‑1 (MQ=255)

        gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGt                                                                                                                                              <  1:7441787/150‑1 (MQ=255)

        gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGt                                                                                                                                              <  2:2286259/150‑1 (MQ=255)

        gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                              >  1:870168/1‑150 (MQ=255)

        gCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGt                                                                                                                                              >  1:7785704/1‑150 (MQ=255)

         cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTc                                                                                                                                             <  2:5859380/150‑1 (MQ=255)

         cTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGTCGGCGACCATGTc                                                                                                                                             <  1:8181140/150‑1 (MQ=255)

          tGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                                            >  2:4329124/1‑150 (MQ=255)

          tGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCa                                                                                                                                            <  2:1805713/150‑1 (MQ=255)

               gACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACCTGGTCGTCGGCGACCATGTCAcgtag                                                                                                                                       <  2:2868199/150‑5 (MQ=255)

                 cTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCa                                                                                                                                     <  2:5076956/150‑1 (MQ=255)

                  tACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAg                                                                                                                                    >  1:8574696/1‑150 (MQ=255)

                  tACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAg                                                                                                                                    <  2:7619931/150‑1 (MQ=255)

                   aCCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAgg                                                                                                                                   <  1:3588541/150‑1 (MQ=255)

                    ccTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGc                                                                                                                                  >  2:2036435/1‑150 (MQ=255)

                    ccTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGc                                                                                                                                  >  1:6029125/1‑150 (MQ=255)

                     cTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGcc                                                                                                                                 >  1:7316902/1‑150 (MQ=255)

                      tGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGccc                                                                                                                                >  1:6888018/1‑150 (MQ=255)

                         aCATGAACCGGGTGCAAAGGGGGCGGCGGGATGATGGGTTGCGCGAGCACCTGGCCTGGCTGTGCAAGACTCAACGG‑CCCGGCTACCCCGACTGGCTGGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgc                                                                                                                             <  1:154663/150‑1 (MQ=255)

                         aCATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgc                                                                                                                             >  1:5298122/1‑150 (MQ=255)

                            tGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgcagc                                                                                                                          <  1:547897/150‑1 (MQ=255)

                             gAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTgcagca                                                                                                                         >  1:8481213/1‑150 (MQ=255)

                               aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                       <  1:1435971/150‑1 (MQ=255)

                               aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                       <  2:8007618/150‑1 (MQ=255)

                               aCCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCaaa                                                                                                                       <  1:4861610/150‑1 (MQ=255)

                                ccGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAg                                                                                                                      >  1:4288415/1‑150 (MQ=255)

                                  gcgAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCa                                                                                                                    >  1:2720423/1‑150 (MQ=255)

                                     aGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCAccc                                                                                                                 <  1:562266/150‑1 (MQ=255)

                                     aGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCAccc                                                                                                                 <  1:6498834/150‑1 (MQ=255)

                                     aGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCAccc                                                                                                                 <  1:3093577/150‑1 (MQ=255)

                                      gCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCcg                                                                                                                >  1:3349709/1‑150 (MQ=255)

                                        aaaGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCcgcg                                                                                                              <  1:6455553/150‑1 (MQ=255)

                                         aaGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGa                                                                                                             <  1:6996629/150‑1 (MQ=255)

                                             gcgcTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTc                                                                                                         <  2:8335982/150‑1 (MQ=255)

                                                  gctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                                    <  2:5449183/150‑1 (MQ=255)

                                                  gctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                                    <  2:3361685/150‑1 (MQ=255)

                                                  gctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGa                                                                                                    <  2:6731393/150‑1 (MQ=255)

                                                   ctgGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAg                                                                                                   <  1:2055175/150‑1 (MQ=255)

                                                     ggATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGc                                                                                                 >  2:4755980/1‑150 (MQ=255)

                                                         gatCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAAtt                                                                                             <  2:3888836/150‑1 (MQ=255)

                                                           tCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTgg                                                                                           >  2:8156006/1‑150 (MQ=255)

                                                             ggTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                         <  1:8377802/150‑1 (MQ=255)

                                                             ggTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTc                                                                                         <  2:5830719/150‑1 (MQ=255)

                                                                tGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAt                                                                                      <  2:7340186/150‑1 (MQ=255)

                                                                tGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAt                                                                                      <  2:6938330/150‑1 (MQ=255)

                                                                 gcgcGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGAtt                                                                                     <  2:5774377/150‑1 (MQ=255)

                                                                  cgcgAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTg                                                                                    >  2:4251840/1‑150 (MQ=255)

                                                                   gcgAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTgg                                                                                   >  2:8756544/1‑150 (MQ=255)

                                                                    cgAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGc                                                                                  <  1:7966672/150‑1 (MQ=255)

                                                                         aCCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                             <  2:1194928/150‑1 (MQ=255)

                                                                         aCCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                             >  1:3745075/1‑150 (MQ=255)

                                                                         aCCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                             >  2:2655721/1‑150 (MQ=255)

                                                                         aCCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGCTCGCCGACGACCAGCTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAgg                                                                             >  1:4966545/1‑150 (MQ=255)

                                                                            tggcctggcATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGaac                                                                          >  2:6384682/1‑149 (MQ=255)

                                                                                ctggcATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGTCTCGACCGGCTACGTGACATGGTAGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                                      >  1:188613/1‑150 (MQ=255)

                                                                                ctggcATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAc                                                                      >  1:5700894/1‑150 (MQ=255)

                                                                                 tggcATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAcc                                                                     <  1:2645681/150‑1 (MQ=255)

                                                                                 tggcATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCAcc                                                                     <  2:3224343/150‑1 (MQ=255)

                                                                                      ttGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTc                                                                <  1:7613207/150‑1 (MQ=255)

                                                                                      ttGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTc                                                                <  1:4164508/150‑1 (MQ=255)

                                                                                         cAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGc                                                             >  1:6706961/1‑150 (MQ=255)

                                                                                          aaGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGcc                                                            >  1:946146/1‑150 (MQ=255)

                                                                                          aaGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGcc                                                            >  1:1162958/1‑150 (MQ=255)

                                                                                           aGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCa                                                           >  2:7084111/1‑150 (MQ=255)

                                                                                           aGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCa                                                           <  2:5028717/150‑1 (MQ=255)

                                                                                            gACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAc                                                          <  1:2593441/150‑1 (MQ=255)

                                                                                             atgcAACTGACCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAc                                                          >  1:4048674/4‑150 (MQ=255)

                                                                                             aCTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAcc                                                         >  1:1214544/1‑150 (MQ=255)

                                                                                               tCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGa                                                       >  2:8238933/1‑150 (MQ=255)

                                                                                               tCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGa                                                       >  1:5182602/1‑150 (MQ=255)

                                                                                               tCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGa                                                       >  1:1051441/1‑150 (MQ=255)

                                                                                                 aaCTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAAc                                                     >  1:3445283/1‑150 (MQ=255)

                                                                                                    tG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAActgc                                                  >  2:1594540/1‑150 (MQ=255)

                                                                                                        ccGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACtgctgc                                               <  2:5273226/150‑1 (MQ=255)

                                                                                                               ccTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCAt                                        <  2:5651842/150‑1 (MQ=255)

                                                                                                                 tCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGa                                      >  1:1676724/1‑150 (MQ=255)

                                                                                                                  cGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAg                                     <  1:6120725/150‑1 (MQ=255)

                                                                                                                   gACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGc                                    >  2:2630010/1‑150 (MQ=255)

                                                                                                                   gACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGc                                    <  2:6900548/150‑1 (MQ=255)

                                                                                                                   gACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGc                                    >  2:148927/1‑150 (MQ=255)

                                                                                                                       ggCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACa                                >  1:808024/1‑150 (MQ=255)

                                                                                                                        gCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACaa                               >  1:6207593/1‑150 (MQ=255)

                                                                                                                               cGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTa                        <  2:6286170/150‑1 (MQ=255)

                                                                                                                                ggCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAg                       <  2:7430068/150‑1 (MQ=255)

                                                                                                                                 gCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgc                      <  2:292903/150‑1 (MQ=255)

                                                                                                                                   tACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTaacag                    >  1:2176055/1‑146 (MQ=255)

                                                                                                                                   tACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAgcgc                    >  2:4052235/1‑150 (MQ=255)

                                                                                                                                    aCGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCa                   <  1:6620834/150‑1 (MQ=255)

                                                                                                                                      gTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGa                 <  1:2265926/150‑1 (MQ=255)

                                                                                                                                      gTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGa                 <  2:8574696/150‑1 (MQ=255)

                                                                                                                                        gACCTGGTCGTCGGCGACCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATg               >  2:1439039/1‑150 (MQ=255)

                                                                                                                                        gACCTGGTCGTCGGCGACCATGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATg               >  2:7211699/1‑150 (MQ=255)

                                                                                                                                          cATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCg             >  1:4751619/1‑150 (MQ=255)

                                                                                                                                          cATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCg             >  2:354559/1‑150 (MQ=255)

                                                                                                                                               tCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCAc                                                          <  2:6505740/100‑1 (MQ=255)

                                                                                                                                               tCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAgg        <  2:1148779/150‑1 (MQ=255)

                                                                                                                                               tCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAgg        >  2:5997912/1‑150 (MQ=255)

                                                                                                                                                 gCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGcc      <  2:3387745/150‑1 (MQ=255)

CCGTGGTCGCTGAGCGACTACCTGGACATGAACCGCGAGCAAAGGGCGCTGCTGGATGATCGGTTGCGCGAGCACCTGGCCTGGCATTGCAAGACTCAACTG‑CCCGGCTACCTCGACTGGCTCGACCGGCTACGTGACATGGTCGCCGACGACCAGGTCACCGACCAGGCCCTGCAGCAAAGCACCCGCGAGGTCCGTGAGGCAATTGGTCGATTGGCCGAGGAAATCACCCCTTCGGCCACCGAACTGCTGCGCGGCATGAGCGACAACCAGGTAGCGCAGATGCGTGAGGCCTTCC  >  NC_002947/5188969‑5189266

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