BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,777,043	C→G	5.5%	A169P (GCA→CCA)	PP_2429 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,777,043	0	C	G	5.5%	343.3 / 9.1	109	A169P (GCA→CCA)	PP_2429	membrane protein
Reads supporting (aligned to +/- strand): ref base C (58/44); new base G (4/2); total (63/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 110 total aligned reads displayed.

CGGCGCTTTCGCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCA  >  NC_002947/2776894‑2777187

            gAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGNTCAANTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGGTGACCACAAAACCCCCGCCAACACTGGGTTGGAATCATACTACCTCCCCGGTGCCTGCACCCCCC‑g                                                                                                                                      >  1:641/1‑150 (MQ=255)

             gggCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                     <  2:1318236/149‑1 (MQ=255)

             aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                     <  1:192962/150‑1 (MQ=255)

             aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg                                                                                                                                     <  1:454253/150‑1 (MQ=255)

                 gTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGgcgc                                                                                                                                 >  2:1822878/1‑150 (MQ=255)

                  ttGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCt                                                                                                                                >  1:1602575/1‑150 (MQ=255)

                  ttGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCt                                                                                                                                >  2:243678/1‑150 (MQ=255)

                    gACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGcccc                                                                                                                              >  2:33665/1‑150 (MQ=255)

                     acacGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTccc                                                                                                                             <  1:243678/150‑1 (MQ=255)

                                    ttcGCTCAGGCCGCTTGTTGGTCAAGTATTTTTTCAGTCCAATACTGAAGTCTATGCCGATATCAACAAAACCCCCGCCAACAGTTGGTTGAATATTCGCTACCTCTCCTTTTCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGc                                                                                                              >  2:451171/3‑150 (MQ=255)

                                      cGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGccc                                                                                                            >  2:1686543/1‑150 (MQ=255)

                                       gatCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGTAGCCCATGCCGAGATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCa                                                                                                           >  2:1756595/3‑150 (MQ=255)

                                        tgcAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCa                                                                                                                                              >  2:517327/3‑115 (MQ=255)

                                        tgcAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCa                                                                                                                                              <  1:517327/113‑1 (MQ=255)

                                             gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGGCTGCAGGCACC‑ggg                                                                                                                                    <  1:574278/119‑1 (MQ=255)

                                              ccccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACCCCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAactc                                                                                                                 <  2:71603/137‑3 (MQ=255)

                                               cccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCa                                                                                                   >  2:1831942/1‑150 (MQ=255)

                                               cccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCa                                                                                                   <  1:1233222/150‑1 (MQ=255)

                                                ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                  >  2:1495817/1‑150 (MQ=255)

                                                ccTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt                                                                                                  >  1:114689/1‑150 (MQ=255)

                                                        gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGcc                                                                                                             <  1:93493/131‑1 (MQ=255)

                                                        gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGcc                                                                                                             >  2:93493/1‑131 (MQ=255)

                                                        gTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCa                                                                                          <  1:121512/150‑1 (MQ=255)

                                                        gTCAAGTAGAGGCGCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCa                                                                                          <  1:322578/150‑1 (MQ=255)

                                                         tCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCag                                                                                         >  2:1036884/1‑150 (MQ=255)

                                                           aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTcccccccccc                                                                                                                      <  1:482478/119‑1 (MQ=255)

                                                           aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag                                                                                       >  1:1050163/1‑150 (MQ=255)

                                                           aaGTAGAGGCGCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCAccccc                                                                                                                      >  2:482478/1‑119 (MQ=255)

                                                             gTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                     >  1:622509/1‑150 (MQ=255)

                                                             gTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                     >  1:980744/1‑150 (MQ=255)

                                                             gTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCg                                                                                     >  2:547349/1‑150 (MQ=255)

                                                                  ggCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAg                                                                                >  1:1794241/1‑150 (MQ=255)

                                                                   gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                 >  1:1155615/1‑148 (MQ=255)

                                                                   gCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa                                                                                 <  2:1155615/148‑1 (MQ=255)

                                                                    cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCAATGCTCGGCAGAGCGATCAGc                                                                               >  1:1913806/1‑149 (MQ=255)

                                                                    cTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCAATGCTCGGCAGAGCGATCAGc                                                                               <  2:1913806/149‑1 (MQ=255)

                                                                     tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTc                                                                                                                               >  1:516472/1‑100 (MQ=255)

                                                                     tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTc                                                                                                                               <  2:516472/100‑1 (MQ=255)

                                                                     tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGCGCGATCCGCaa                                                                             >  2:856287/1‑150 (MQ=255)

                                                                     tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                             <  2:1794241/150‑1 (MQ=255)

                                                                      cAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAg                                                                            <  1:1495112/150‑1 (MQ=255)

                                                                        ggCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGaa                                                                          <  1:1036884/150‑1 (MQ=255)

                                                                           cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAacac                                                                       <  2:114689/150‑1 (MQ=255)

                                                                           cAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTAAATGAAGTCACC‑GAGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAacac                                                                       >  1:1875487/1‑150 (MQ=255)

                                                                               aCTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgt                                                                   >  1:1897686/1‑149 (MQ=255)

                                                                               aCTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgt                                                                   >  1:722311/1‑149 (MQ=255)

                                                                                 tGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgcgc                                                                 >  2:1258538/1‑150 (MQ=255)

                                                                                      cccATGCCGATATCAACAAAACCCCCGCCAACCCTTGGTTGAACACCCGCACCCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt                                                            <  2:384492/150‑1 (MQ=255)

                                                                                      cacATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt                                                            >  2:472027/3‑150 (MQ=255)

                                                                                       ccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                             >  2:1639918/1‑132 (MQ=255)

                                                                                       ccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa                                                                             <  1:1639918/132‑1 (MQ=255)

                                                                                        cATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGc                                                          <  2:199479/150‑1 (MQ=255)

                                                                                        cATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGc                                                          >  2:619604/1‑150 (MQ=255)

                                                                                        cATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGc                                                          <  1:472027/150‑1 (MQ=255)

                                                                                         aTGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCa                                                                                                                                              >  2:153428/1‑66 (MQ=255)

                                                                                         aTGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCa                                                                                                                                              <  1:153428/66‑1 (MQ=255)

                                                                                          tGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCca                                                        <  1:1201567/150‑1 (MQ=255)

                                                                                               atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                   <  2:822782/150‑1 (MQ=255)

                                                                                               atatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTg                                                   <  1:1258538/150‑1 (MQ=255)

                                                                                                tatCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTgg                                                  <  1:619604/150‑1 (MQ=255)

                                                                                                 atCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                 >  1:1428600/1‑150 (MQ=255)

                                                                                                   caacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCaa                                               <  1:1831942/150‑1 (MQ=255)

                                                                                                    aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCca                                                        >  2:1433267/1‑140 (MQ=255)

                                                                                                    aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              >  2:1179532/1‑150 (MQ=255)

                                                                                                    aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              >  2:171563/1‑150 (MQ=255)

                                                                                                     acaaAACCCCCGTCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTCCCACACTGTCAacg                                             >  2:1384640/1‑150 (MQ=255)

                                                                                                     acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg                                                                                                      >  1:1823654/1‑93 (MQ=255)

                                                                                                     acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg                                                                                                      <  2:1823654/93‑1 (MQ=255)

                                                                                                     acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAacac                                                                       <  1:618875/124‑1 (MQ=255)

                                                                                                     acaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAacac                                                                       >  2:618875/1‑124 (MQ=255)

                                                                                                        aaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgacg                                          <  2:1875487/150‑1 (MQ=255)

                                                                                                         aaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGGCCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACgg                                         >  1:904167/1‑150 (MQ=255)

                                                                                                            ccccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcgc                                      <  2:1602575/150‑1 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              <  2:939467/141‑1 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              <  2:938631/141‑1 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              >  1:938631/1‑141 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAac                                              >  1:939467/1‑141 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGcgcg                                     >  2:196404/1‑150 (MQ=255)

                                                                                                             cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGcgcg                                     >  1:920710/1‑150 (MQ=255)

                                                                                                               cGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc                                   <  2:1413347/150‑1 (MQ=255)

                                                                                                                   aaCACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAAt                               <  1:1686543/150‑1 (MQ=255)

                                                                                                                            ttGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                 >  1:238622/1‑123 (MQ=255)

                                                                                                                            ttGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGc                                                 <  2:238622/123‑1 (MQ=255)

                                                                                                                              gAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCCAGCACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGaccac                    >  2:222661/1‑150 (MQ=255)

                                                                                                                                acacCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCg                  >  1:217049/1‑150 (MQ=255)

                                                                                                                                   cccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGccc               >  2:367576/1‑149 (MQ=255)

                                                                                                                                   cccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCa               >  2:433915/1‑150 (MQ=255)

                                                                                                                                    ccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc              >  1:950806/1‑150 (MQ=255)

                                                                                                                                    ccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc              >  1:1387969/1‑150 (MQ=255)

                                                                                                                                    ccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc              >  2:8944/1‑150 (MQ=255)

                                                                                                                                    ccGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCAc              >  2:1619567/1‑150 (MQ=255)

                                                                                                                                           cGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATc       >  2:835317/1‑150 (MQ=255)

                                                                                                                                            cgcccGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcg                                       <  2:1777555/116‑1 (MQ=255)

                                                                                                                                            cgcccGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGgcg                                       >  1:1777555/3‑118 (MQ=255)

                                                                                                                                              cccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc    >  2:1672385/1‑150 (MQ=255)

                                                                                                                                              cccGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc    <  2:1387969/150‑1 (MQ=255)

                                                                                                                                               ccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAcc   >  1:757071/1‑150 (MQ=255)

                                                                                                                                                cGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCa  <  1:1495817/150‑1 (MQ=255)

                                                                                                                                                cGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCa  <  1:171563/150‑1 (MQ=255)

                                                                                                                                                cGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCCGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCa  <  1:1756595/150‑1 (MQ=255)

CGGCGCTTTCGCGAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAACCA  >  NC_002947/2776894‑2777187

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