BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,134,010	G→C	5.2%	W348C (TGG→TGC) ‡	tolC →	agglutination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,134,010	0	G	C	5.2%	358.9 / 9.3	115	W348C (TGG→TGC) ‡	tolC	agglutination protein
Reads supporting (aligned to +/- strand): ref base G (63/46); new base C (3/3); total (66/49)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 116 total aligned reads displayed.

TGGCAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTAC  >  NC_002947/5133861‑5134157

tGGCAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGCCATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTgg                                                                                                                                                       >  1:321773/1‑150 (MQ=255)

   cAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCCAGAC‑TTCGCAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAAc                                                                                                                                                    <  1:42296/150‑1 (MQ=255)

         ctcgttctcactgtcgaggacgtccagcagggtgcgctggcccaggccgaactggtcctggtaggccgcacgcacacgcttggtggtttcggcgtattcacgcgcggtcggcagctgcttgctggcgttgttcaTGGCGTGGAACGCCAt                                                                                                                                               >  1:1271915‑M1/135‑150 (MQ=255)

              ggAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGaacaa                                                                                                                                         <  1:1586120/150‑1 (MQ=255)

                agaTGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑g                                                                                                                                     >  2:1420128/1‑150 (MQ=255)

                   tGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGCGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCa                                                                                                                                  >  1:337483/1‑150 (MQ=255)

                    gAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCa                                                                                                                                                >  1:617372/1‑137 (MQ=255)

                      aCTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagca                                                                                                                               >  1:138865/1‑150 (MQ=255)

                      aCTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagca                                                                                                                               >  1:140517/1‑150 (MQ=255)

                        tACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGAACGAAAACCTGAACCTGGCGTGGAACGCCATAAACAA‑C‑GCCagcaaa                                                                                                                             >  2:144143/1‑149 (MQ=255)

                        tACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAAAGCCCCGCGTAAATTGCCCGAGAACCGCTGCCGGGCGTGGAACCCCATGAACAA‑C‑GCCagcaag                                                                                                                             >  2:2365068/1‑150 (MQ=255)

                             ccTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCt                                                                                                                        >  1:2299769/1‑150 (MQ=255)

                              cTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTg                                                                                                                       <  1:504163/150‑1 (MQ=255)

                               tGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGc                                                                                                                      <  1:2339439/150‑1 (MQ=255)

                                gTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGcc                                                                                                                     <  2:652321/150‑1 (MQ=255)

                                             aGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAATATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCa                                                                                                                                                <  2:1380746/112‑1 (MQ=255)

                                             aGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCa                                                                                                                                                >  1:1380746/1‑112 (MQ=255)

                                              gCGACAAGGCTCGCATACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAAt                                                                                                       <  2:1133549/150‑1 (MQ=255)

                                                 aCAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACg                                                                                                    >  1:1640617/1‑150 (MQ=255)

                                                  cAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCACGCGTGAATACGc                                                                                                   <  1:172829/150‑1 (MQ=255)

                                                    aGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCg                                                                                                 <  2:1640617/150‑1 (MQ=255)

                                                     ggCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGa                                                                                                >  1:499415/1‑150 (MQ=255)

                                                      gCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaa                                                                                               >  2:1375603/1‑150 (MQ=255)

                                                       cTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                              <  2:138865/150‑1 (MQ=255)

                                                        tCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACCACGTACTGCGTTAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATGCGCCGAAac                                                                                             >  2:919694/1‑150 (MQ=255)

                                                        tCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAac                                                                                             <  2:1467985/150‑1 (MQ=255)

                                                           ccTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATAc                                                                                                     >  2:913423/1‑139 (MQ=255)

                                                            cTACAGTCC‑GACGCGCACAAGATCCACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAaccacc                                                                                         >  2:1622817/1‑150 (MQ=255)

                                                             tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGTAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAAccacca                                                                                        >  2:1484252/1‑150 (MQ=255)

                                                             tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCTTGAATACGCCGAAAccacca                                                                                        >  1:1062626/1‑150 (MQ=255)

                                                             tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAAccacca                                                                                        >  2:1875767/1‑150 (MQ=255)

                                                               cAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                              <  1:944307/142‑1 (MQ=255)

                                                               cAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                              >  2:944307/1‑142 (MQ=255)

                                                                aGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGc                                                                                     <  1:1622817/150‑1 (MQ=255)

                                                                    c‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCTCGTGAATACGCCGAAACCACCAAGCgtgt                                                                                 >  1:2192340/1‑150 (MQ=255)

                                                                      gACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCGTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGtgtg                                                                                <  2:1655536/150‑1 (MQ=255)

                                                                      gACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGtgtg                                                                                >  2:1488823/1‑150 (MQ=255)

                                                                         gcgcACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTgtgcgt                                                                             >  1:1371799/1‑150 (MQ=255)

                                                                           gcACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTgtgcgtgc                                                                           <  1:144143/150‑1 (MQ=255)

                                                                             acaAGATCAACCAGGCACTGGACATTCGGAACAACGCCCTGCGTGAGTTGACCGAGGACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                         >  2:1255870/1‑150 (MQ=255)

                                                                             acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                         >  2:491432/1‑150 (MQ=255)

                                                                             acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                         >  2:1931310/1‑150 (MQ=255)

                                                                             acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                         >  2:2113566/1‑150 (MQ=255)

                                                                             acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                         >  2:1978671/1‑150 (MQ=255)

                                                                                 gATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTa                                                                     >  1:1622746/1‑150 (MQ=255)

                                                                                 gATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTa                                                                     >  1:1030462/1‑150 (MQ=255)

                                                                                   tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGCCGAATCTTACTGTATGCAGGTG‑‑GACCACGCGTTAAGACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTAcc                                                                   >  1:280203/1‑150 (MQ=255)

                                                                                   tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaa                                                                                                                              <  1:1684237/91‑1 (MQ=255)

                                                                                   tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaa                                                                                                                              >  2:1684237/1‑91 (MQ=255)

                                                                                     aaCCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCCTGTGCGTGCGGCATACCAg                                                                 >  2:228241/1‑150 (MQ=255)

                                                                                       ccAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAggc                                                               >  2:1010572/1‑149 (MQ=255)

                                                                                        cAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTAGCACCACCTGTTTCC‑T‑GCTAGCAAGCAGCTGCTTTCCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                              <  1:1484252/150‑1 (MQ=255)

                                                                                             aCTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGtt                                                         >  2:630527/1‑150 (MQ=255)

                                                                                              cTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTc                                                        >  1:996431/1‑150 (MQ=255)

                                                                                               tGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCg                                                       >  2:79721/1‑150 (MQ=255)

                                                                                               tGGACATTCGCAACAACGCCCGGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCg                                                       <  2:227995/150‑1 (MQ=255)

                                                                                                 gACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGc                                                     <  2:499415/150‑1 (MQ=255)

                                                                                                    aTTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTg                                                  <  2:51288/150‑1 (MQ=255)

                                                                                                    aTTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTg                                                  >  2:2306714/1‑150 (MQ=255)

                                                                                                     ttCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGTTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                 >  2:1582682/1‑150 (MQ=255)

                                                                                                     ttCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                 <  1:79721/150‑1 (MQ=255)

                                                                                                     ttCGCAACAACACCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                 <  1:1073434/150‑1 (MQ=255)

                                                                                                      tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaagca                                                                                                                           >  1:336468/1‑75 (MQ=255)

                                                                                                      tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaagca                                                                                                                           <  2:336468/75‑1 (MQ=255)

                                                                                                      tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTggg                                                >  2:913793/1‑150 (MQ=255)

                                                                                                         caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCa                                                                  >  1:652528/1‑129 (MQ=255)

                                                                                                         caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCa                                                                  <  2:652528/129‑1 (MQ=255)

                                                                                                         caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCa                                             <  1:1875767/150‑1 (MQ=255)

                                                                                                          aacaacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAg                                            >  1:80729/1‑150 (MQ=255)

                                                                                                           acaacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCGGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAgc                                           <  1:1978671/150‑1 (MQ=255)

                                                                                                             aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                              <  2:847181/129‑1 (MQ=255)

                                                                                                             aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                              >  1:847181/1‑129 (MQ=255)

                                                                                                             aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAgcgc                                         >  1:948563/1‑150 (MQ=255)

                                                                                                               cGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCAc                                       >  2:1900660/1‑150 (MQ=255)

                                                                                                                gCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCAcc                                      >  2:1129779/1‑150 (MQ=255)

                                                                                                                gCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCAcc                                      <  1:630527/150‑1 (MQ=255)

                                                                                                                    tGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGACCAC‑A‑GTCCGTGGTCAGCCGCTCACCCCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCctgc                                  <  1:2113566/150‑1 (MQ=255)

                                                                                                                        tGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGa                              <  1:1900660/150‑1 (MQ=255)

                                                                                                                              gACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCt                        <  1:2306714/150‑1 (MQ=255)

                                                                                                                              gACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCt                        >  1:2348062/1‑150 (MQ=255)

                                                                                                                               ttgctggcgttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGa                                                               >  2:2200573‑M1/17‑110 (MQ=255)

                                                                                                                               ttgctggcgttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGa                                                               <  1:2200573‑M1/94‑1 (MQ=255)

                                                                                                                               aCCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                 >  1:964302/1‑124 (MQ=255)

                                                                                                                                 cGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAcc                                                             >  2:2199416/1‑110 (MQ=255)

                                                                                                                                   agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                   >  1:2224784/1‑150 (MQ=255)

                                                                                                                                   agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                   <  2:2224784/150‑1 (MQ=255)

                                                                                                                                   agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                   <  2:80729/150‑1 (MQ=255)

                                                                                                                                       gttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGt                                                                                                           >  2:1522919‑M1/9‑58 (MQ=255)

                                                                                                                                       gttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGt                                                                                                           <  1:1522919‑M1/50‑1 (MQ=255)

                                                                                                                                       ccTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTga               >  2:2379812/1‑150 (MQ=255)

                                                                                                                                        cTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTgag              >  2:1175378/1‑150 (MQ=255)

                                                                                                                                        cTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTgag              >  2:1632117/1‑150 (MQ=255)

                                                                                                                                           ttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGt                 >  1:195413‑M1/5‑144 (MQ=255)

                                                                                                                                           ttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGt                 <  2:195413‑M1/140‑1 (MQ=255)

                                                                                                                                             ccTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCTTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGa         <  2:1062626/150‑1 (MQ=255)

                                                                                                                                              cTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAg        <  2:948563/150‑1 (MQ=255)

                                                                                                                                               TGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGc       <  2:639085‑M1/150‑1 (MQ=255)

                                                                                                                                                GGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCt      >  1:1255870‑M1/1‑150 (MQ=255)

                                                                                                                                                   GTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATATGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTa   >  2:1765782‑M1/1‑150 (MQ=255)

                                                                                                                                                    TGGACCGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTAc  >  1:111261‑M1/1‑150 (MQ=255)

TGGCAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTAC  >  NC_002947/5133861‑5134157

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